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General scalarization method of dynamic elastic fields in transversally isotropic media and its 


new applications’ 
I. P. Miroshnichenko’, V. P. Sizov” 


"Don State Technical University, Rostov-on-Don, Russian Federation 


? Rostov Scientific Research Institute for Radiocommunication, Rostov-on-Don, Russian Federation 


OdoOmenHEI MeTO, CKasApH3allun AWHAMHUYeCKHX YIPyrHx Noel B TpancBepCcaJIbHO-H30TPONHBIX cpeyqax A 


cer 
ero HOBBbIe IIpHMeHeHHA 


M. 1. Mupomnuuenko’, B. I. Cu30B" 


‘JloncKoii rocyqapcTBennpiit TexHuueckuit yHuBepcuter, r. Pocros-Ha-Jlony, Poccutickas Denepanna 


? PocroscKuii-Ha-Jlony Hay4HO-HcceqOBaTeNbCKUi MHCTUTYT pagMocsa3u, r. Poctos-Ha-Jlony, Poccniicxas Dezepanus 


Introduction. An efficient technique of tensor field scalariza- 
tion is successfully used while investigating tensor elastic 
fields of displacements, stresses and deformations in the lay- 
ered structures of different materials, including transversally 
isotropic composites. These fields can be expressed through 
the scalar potentials corresponding to the quasi-longitudinal, 
quasi-transverse, and transverse-only waves. Such scalariza- 
tion is possible if the objects under consideration are tensors 
relating to the subgroup of general coordinate conversions, 
when the local affine basis has one invariant vector that coin- 
cides with the material symmetry axis of the material. At this, 
the known papers consider structures where this vector coin- 
cides with the normal to the boundary between layers. Howev- 
er, other cases of the mutual arrangement of the material 
symmetry axis of the material and the boundaries between 
layers are of interest on the practical side. 

Materials and Methods. The work objective is further devel- 
opment of the scalarization method application in the bounda- 
ry value problems of the dynamic elasticity theory for the 
cases of an arbitrary arrangement of the material symmetry 
axis relative to the boundary between layers. The present re- 
search and methodological apparatus are developed through 
the general technique of scalarization of the dynamic elastic 
fields of displacements, stresses and strains in the transversally 
isotropic media. 

Research Results. New design ratios for the determination of 
the displacement fields, stresses and deformations in the trans- 
versally isotropic media are obtained for the cases of an arbi- 
trary arrangement of the material symmetry axes of the layer 


materials with respect to the boundaries between layers. 





Beedenue. pu uccheqoBaHuu TeH30pHbIX yIIpyrux Moller 
TlepeMeleHHii, HalipskKeHH HW ZePpopMalHii B CIIOHCTBIX KOH- 
CTPyKUMAX V3 pa3JIMYHbIX MaTepHasiOB, BKJIKOUad TpaHCBep- 
CaJIbHO-H30TpONHbIe KOMIO3HTbI, YCIeWIHO MpHMeHsAeTCA 
oekTHBHBIM MeTO CKaJlapH3allHH TeH30pHBIX Hostel. Jlay- 
HBbIe MOI MOTT OBITS BbIPAKeHLI Yepe3 CKaJIApHble NOTeHIH- 
alIbl, COOTBCTCTBYIOWHe KBa3HIIPOMOJIbHBIM, KBa3HUOMepe4- 
HbIM HM 4YHCTO MOMepedHbIM BOJIHAM. Takad cKaslapH3alua 
BO3MO2%KHA, CCIM paCCMaTpHBaeMbIe OOBEKTHI ABJIAIOTCA TeH- 
30PaMH OTHOCHTeIbHO MOArpynibl oOmjux mpeoOpa3z0BaHHit 
KOOpAMHAT, KOra JIOKaJIbHEIM adPuHHBI Oa3uc HMeeT ODMH 
MHBapHaHTHbIM BeKTOP, KOTOPbIM COBMayaeT C OCbEIO MaTepH- 
ayIbHOHM CHMMeTpHH Matepuasia. IIpu 3TOM B H3BeCTHBIX pado- 
TaX PpaCCMaTPHBAIOTCA KOHCTPYKIUMH, Pye TOT BEKTOP COBIa- 
qaeT C HOpMasIbiO K rpaHule MexKLYy closAMH. OyHakKO, WA 
IIpakTHKH NpeACTaBIAIOT HHTepec UM Apyrve cilyyav B3aHMHO- 
TO pacHoslo2xkeHHA OCH MaTepHasIbHOH CHMMeTpH MaTepuala 
VW TpaHHbl MexK Ly COAMH. 

[lesbo ABJLIeTCA JabHeliee pa3sBuTHe IPHMeHeHHA MeTOya 
CKaJIApH3alluH B rpaHWyHbIx 3ajayax TWHaMMyuecKOH TeopHu 
yUupyrocTu Ha CJly¥aH TPOW3BOJIbHOrO pacioOxKeHHA OCH 
MaTepHasIbHOH CHMMETPHH M0 OTHOMICHHEO K rpaHulle Mex Ly 
CIOAMH. 

Memoovi uccneOo6aHua. II[peqmaraempiii Hay4dHo- 
MeTOWMYeCKHH allmapaT pa3spaboTaH Ha OCHOBe HCIOJIb30Ba- 
Hua OOOOMjeHHOrO MeTONa cCKallapw3alMH WMHaMM4yecKHx 
yiipyrux oyleli WepeMeleHuii, HalpwKeHHH u TedopmMalHii B 
TPaHCBepCaJIbHO-H30TPOMHBIX Cpeslax. 

Pesyiomamel ucciedoeanun. UomyueHbl HOBbIe pac4eTHbIe 
COOTHOMCHHA Id OlpeyeleHua nose lepemMenieHuit, 
HallpwKeHHM H TePopMalluii B TpaHcBepcasIbHO-H30TPOMHBIX 
cpeyax Ha Cily4av MpOW3BOJbHOLO pacioso2xKeHHA Oceli MaTe- 
PHaIbHOH CHMMeTPHH MaTepHasIOB CyIOeB MO OTHOLICHHIO K 


TpaHhllaM MexK Ty CJIOAMH. 


* The research is done within the frame of government task of RF Ministry of Education and Science no. 9.9770.2017/ BY and is supported in part by 
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Discussion and Conclusions. The present research and meth- 
odological apparatus are successfully used in determining the 
stress-strain state in the layered structures of transversally 
isotropic materials, and in analyzing the diagnosis results of 
the state of the plane-layered and layered cylindrical structures 


under operation. 


Keywords: scalarization method, transversally isotropic me- 
dium, acoustic waves, composite materials. 


For citation: 1.P. Miroshnichenko, V.P. Sizov. General sca- 
larization method of dynamic elastic fields in transversally 


O6cyacdenue u  3akmovenun. Tlpeqnaraempiii Hay4Ho- 
MeTOAMYeCKHH allMapaT ycnelwHO MCiOb30BaH UpH ompeye- 
JICHHM Hallps2%KeHHO-ZetOpMupoOBaHHOrO COCTOAHHA B CJIOH- 
CTBIX KOHCTPyKUMAX, BBINIOJIHCHHBIX V3 TPaHCBepcasIbHO- 
W30TPOMHBIX MaTepuasloB, H Ip aHasIM3e pe3yiIbTaTOB jHa- 
THOCTHKH COCTOAHMA MJIOCKOCOMCTBIX H CJIOMCTbIX UMJIMH- 
UpW4eckHX KOHCTPyKUMH, HaXOAAMIMXCA B OKCIUIyaTalMn. 


Kusrrouesbie C10Ba: MeTOT CKaslapH3auwHH, TpaHCBepCcasIbHo- 
M30TpOnHad Cpeyla, aKYCTHYCCKHE BOJIHbI, KOMITO3HWHOHHBIC 
MaTepHaJibi. 


O6paszeu 01a YUmupoeanua: Mupournnuenko, VU. I. O606- 
INCHHBIM MeTO] CKaslApH3allMu JMHAaMWYeCKHMX yiIpyrHux Mosel 


isotropic media and its new applications. Vestnik of DSTU, 
2018, vol. 18, no.3, pp. 258-264. 
https://doi.org/10.23947/1992-5980-2018-18-3-258-264 


B TpaHCBepCasIbHO-H30TPOMHBIX Cpeyax H erO HOBbIe MpHMe- 
Henna / V1. Il. Mupoumuuenxo, B. I. Cu30s // Becruux J[ou- 
CKOrO roc. TexH. yH-Ta. — 2018. — T.18, Ne3. — C. 258-264. 
https://doi.org/10.23947/1992-5980-2018-18-3-258-264 


Introduction. When investigating the displacement, stress, and strain tensor elastic fields in layered struc- 
tures of various materials, including transversely isotropic composites [1-3], an efficient method for scalarization 
of tensor fields proposed in [4] is successfully used. Then, the checked fields can be expressed through scalar po- 
tentials corresponding to quasi-longitudinal, quasi-transverse, and transverse-only waves, respectively. This sca- 
larization is possible if the objects under consideration are tensors relative to the subgroup of general coordinate 
transformations, when the local affine basis has one invariant vector that coincides with the material symmetry 
axis of the material. Structures where this vector coincides with the normal to the boundary between layers are 
considered in [5]. However, on the practical side, other cases of the mutual arrangement of the material symmetry 
axis of the material and the boundaries between layers are of interest. 

Materials and Methods. The work objective is the further development of the application of the scalari- 
zation method for the boundary problems of the dynamic elasticity theory to the cases of an arbitrary arrangement 
of the material symmetry axis with respect to the boundary between layers. 

At first, in the coordinate system with an admissible reference point, there are components of the dis- 
placement and stress tensors relative to this frame. Then, knowing all the components of the tensor fields in the 
given reference coordinates, we can find the normal and tangential displacement and stress components corre- 
sponding to the area lying on the boundary by passing to the coordinate system associated with the boundary be- 
tween layers. These components are used later to satisfy the boundary problem conditions [5]. 

To dig deeper into the foregoing, we will discuss it in the context of solving problems for the plane- 
layered structures. 

Let us consider the case when the principal symmetry axis of a transversally isotropic material makes an- 
gle a with the normal of the plane boundary surface between the layers (Fig. 1). At this, we assume that the field 


is independent of the coordinate y : 





ad 2 


Fig. 1 Computational pattern 
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Design ratios 
In accordance with [4]: 


(L) (L) 
u,-(9,+8/B.v, Jee] Davy, +25) w+ [e(SkV" -8°V* yy; 


& 
(L) Ze 7 (L) 
o,=(d, g, +4, 878" 486%, +d,V,V,)o+ (1) 


(T) 
+(d, g,V, +d,8! 3y Vv; a av +d, V,V,V,)w+ 


+2,/g[a,(V34V" V5) V*) +4,(8,54V" —8/5" VV, Jv. 


Components of the tensor fields of displacements U; and stresses 6, in the coordinate system associ- 





ated with the anisotropy of the material x”,x“ can be written as follows: 
() 





(L) (L) (L) 7) E? (7) 
U, =1£(D+D,b+(-D, ate)w; 
& 
(T)(T) 
CL) (L) a E 
U, =iB D,o-D, 1 (T) Ww; 
7 (2) 
(L) (L) (ZL) (L) @) (T) (T) (T) (f) (ry @) 
eee (d,+d,)é ]o+i E(d,+d,+d,-d,&)w; 
ee ea 7 1a (® 


a 36 B+ SP 2a ac d,&)w; 
“) wy) 2) @) a) @) 
Ox. =(d,+ d,B )o+i §(d,-d, Bw. 
Here, the potentials of quasi-longitudinal @ and quasi-transverse w waves should satisfy the wave equation 


(1) 
with the corresponding wave number g (g or a ), and they have the form: 


(ZL) 

















b= d,e% ae 
(T) (T) (3) 
w= We e Ex’ e B xX ; 
Where &’ +B’ =9°; (Magy? =6, (4) 
w 
The coefficients D and d entering (2) are determined in [5] and are as follows: 
(L) 
ey gs (C; +Cy) 7 
te ee D : 
ap-h Cy, -(2°-h XC, -C;,-Cy) 
Doss cae a 
pee a p-h (2C,+C,)-(g -h Cc, __ § a 
2 (Z) (L) (Ll) (L) (ZL) () 
he op he Cu (g° hy, Ci3 — Cay) gh’ 
“3 
(T) = 
g (C,+C,,) 
D,= 7) SAR (S) 
o'p—h’(2C,+C,)-(2'-#)C, 
(L) (L) (L) (L) (L) (L) (L) (L) (L) (L) (L) 


d,=-—g° a,-h’ a,(D,+D,); d, =—g° a,—N (2a, D,+a,D,+a, D,); 


(Z) (L) (LI) () (L) (L) 


d, = 2a, D,+2a,(2D,+D,); d,=2a,D,; 
(7) 
7 Mm @ @) Z20 
d, =a, g(1-D,) +a, 8-DD > 
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(T) 





(T) (T) (T) (T) (T) h? (T) 
d, =a, g(I-D,)+2a, g+a;| g-—ZD, |; 
& 
(T) 
(7) (1) (T) (1) p? | (1) | 
d, =2a, g+2a,| g-2D. 5 3d,=2a,D, =; (6) 
& & 
1 1 ; 
a, =C,, —2a, ; a =—(C p= C3): m=—-(C 2G, )+C,+C,=a3 
2 2 
1 
a, =A =Ca)* Casa, =C, +C3 —2(C); +2C,,), (7) 


where C;, are elasticity moduli of the material written down by the contracted index [6-7]. 


In the formulas (5) and (6) for these coefficients: 
cos 9 sin 8 g. 
-_ ——— > = = ain oe. > & = bg 
cos(a + 8) cos(a + 8) cos(a + 8) 
The boundary conditions of the dynamic elasticity problems include the components of displace- 


ments U,, U,, and stresses G.,, G,, written in the coordinate system related to the boundary (Fig. 1). 


The reference coordinates Z , ¥ andx’, x“, u, are interlinked through the relations[8]: 
Zz =-x* sina+x’ cosa; x =x* cosa+x’ sina (8) 
or: 
x* =xcosa—Zsina;x” =xXsina+Zcosa. (9) 
Using relations (8) and (9), from the formulas [5]: 
G2 ee e. (10) 
© ae ax! ax! oe 
where the coefficients with a bar tab refer to the coordinate system z, x (let us call them “new” coordinates), 
and without a bar tab — to the coordinate system x”, x“ (“old” coordinates). 


We can write the components of the displacement U,, U, and stressG,,, G,, fields entering the bounda- 


zz? 


ry conditions through the components (2): 
U, =cosaU, —sinaU, ; 
U, =sinaU, +cosau, ; (11) 


G., =COS’ AG, —SiN2AG , +SiN? ASK, ; 
= 1. . 
S.. = sin 2a(oy ~ xq) + (— 2sin’ O)O jx - (12) 


In the “new” coordinates, the potential functions (3) have the form: 


ad © i ) BX i 

i X¥sina+Z cosa) i Xcosa—Zsina) , 
b = hye e ? 
(PT) (T) 
it (Fsina+z cosa) Gf B (¥cosa—Z sina) 


W=We : 


and the displacements (11) are written as follows: 
(T) () 


=. “om @) _ (7) €2 () (L) L) (7) 2) 
U, =cosafi §(D,+ D,)o+(-D, B+ g)w)- sina{i € D,o+(-D, wt s)w] (13) 
& & 
(T) (T)(T) 
- | @ @ _ (1) 2 (1) OW OER _ 
U, =sina[i(D,+ D,)o+(—D, B+ g)w)+cos oli B D, o-D,-— =a Wl (14) 
& & 


Here, the wave numbers €, 8, g are determined with respect to the “old” frame (Fig. 2). 
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Fig. 2.Wave numbers “relations” diagram 


The wave numbers €, B which are the projections of the vector g onto the “old” frame x”, x“, and the 
wave numbers g,, g, that are the projections of the vector g onto the “new frame” Zz, x are interlinked 


through the following relations: 
__cos0 _ _cos® . ga. sind _ sind : 
sin(a +) = cos(a + 0) a sin(a +0) ee cos(a + 0) Be 





=— Ss = _*: _ (15) 
sin(a+8) cos(a+@) 





Substituting (15) into (14), we obtain expressions for displacements in the “new” coordinates: 





= (L) (L) (Lt) _ (fT) (T) 
U_=ig, yee Te eke a meee : 
: cos(a + 8) ° cos(a + 8) 
— © cosOsina “ ©. sinOcosa.— © sin(a+0)cosd ™ sina __ 
Ci oe ee (16) 
cos(a + @) cos(a + 8) cos(a + 8) cos(a + 8) 


where, in accordance with (13) and (15): 
@)_ ) _ )_ 


b= 6,68 e's WH Wels ee, (17) 
Similarly, using (2), (14), and (15), from (12), we obtain the relations for the stress components ente r- 


ing the boundary problem conditions: 























OO @. oO |, 1] © cosbsind | 25.8 © sin? |. 
G.. = {cos’ afd, +d,—(d,+d,) ————- 2’ ] + sin 2a(—d,+ d, ) —————- 2’ + sin’ a(d, -—d, ————- 2’ )} b+ 
ant La eae ) cos(046) =| G .  o05?(014 6) = @ “cost(ateyer? 
+icos a cos@ ere med : cos’ 0 7) nor sin® a 17 © — ¢cos’@ ya 
cos(a +0) Pee * cos*(a +0) ms cos(a +0) Pea? ‘eae 
: Q (T) (1) (1) 29 © | 
+sin 2a = ig.(d,—-d, — gw; 
cos(a + @) cos’ (a + 8) 
1. ® © © e960 © © gin? QQ © : 14 © eosOsind “”. — 
o.. = {sin 2a0[d,—(d,+ d,) ————- 2” >]-(1-2sin’ a)(—d,+d Vbt 
as 5 da (4s ) cos(a+6) = * cos(a+0) °"! ( NG . *cost(a+0)e2?* 
(T) (7) () 2Q_ cin2 Q ©) ; (T) (fT) 2 (T) 
4 gin pO? ia gg, oO dainty gg, — oye (18) 
2 cos(a + 0) cos‘(a+0) ~~ cos(a+0) 2 cos’ (a +0) 


Thus, the displacement and stress components entering the boundary conditions (the coordinate sys- 
temzZ, x is associated with the boundary) are determined by the relations (16) - (18) where the projections of 


the wave vector g in the coordinate systemsZ, ¥ and x’, x* are linked through the relations (15). 


The matrices C describing the wave properties of the layers [9-10] used in [5] and being the key ele- 
ments in the development of concrete solutions to the boundary problems in the layered structures in this case, 


when the axesZ, x form the angle ao with the frame components @”, @* , have the form: 


om i . (19) 


Miroshnichenko I. P., Sizov, V.P. General scalarization method of dynamic elastic fields in transversally isotropic media 





where 


























—— - 2 e*= 0 0 0 
Tt 12 11 12 
Coy -( Co -| igs 0 0 
Cy, Cay Cy) Cay 0 e”™ ies 
C= (L) . Z 
Ci = C3, C39 C= Cy, C39 0 0 e 18:7 0 
3) ~ (4) 7] 
Cy, Cay Cy Can 0 0 7) 
0 -ig.z 
(T) (L) 
£72,785 (20) 
4) © cosOcosa ©) 
C,, =ig,(D,+ 2) > 
cos(a + @) 
(T) () cosa 
C,, = g,(—cos 8D, + ————_) ; 
cos(a + @) 
© cosOsina “ © sinOcosa “ 
C,, =ig,.—\——— (D+ D,) + ———D]; 
cos(a + 8) cos(a + @) 
C 7) sin(a+8)cos @ sina |. 
22 Zz 1 > 
cos(a + @) cos(a + @) 
(ZL) (ZL) (LZ) (LZ) cos” 0 (L) : 14 ©_ eos@sind (L) : (L) (L) sin? a) (L) 
C,, = cos’ afd, + d,—(d,+d,) ———— g? ]+ sin 20(—d,+ d,) — g.++sin’ a(d,—d, —.—— g°) 
cos’ (a +6) 2 cos’ (a +6) cos’ (a +6) 
. cos8 3) ) © OT © Eos*Q sn6 © 1 ©  ¢os’@ 
Cy, = Cos” o ig {dj+d,+d,—d, 2 2.) lg. 37 44 2 &:)+ 
cos(a + 8) cos’ (a+ 0) cos(a + 8) 2 cos’ (a+ 8) 
oD ve 
a coss6 © © sin? % 
+sin” a ig (d, g)s 


cos(a+@) ~* * cos*(a+0)°" 


1 Oo © © ese © © ginr?fof © 1% © cosOsind “ 
C,, =—sin2o[d,—(d,+ d,) ———"—- *]-(1-2sin’ a)(—d,+d a 
0 Sim tala ota) (a+) costar) GO os (a4 0) 


cosO0 .™ © eos’ @-sin’? 0 ™ — sin@ .1 ©  ¢os’0 
,—d, : g.) +(1-2sin* a) 
cos’ (a+) 








(T) 
g:)- (21) 


U 3 4 9 Zz 
cos(a+0) 2 cos’ (a +8) 


Where a =0, the elements (21) of the matrix C coincide with the corresponding elements of this ma- 











trix for the case of coincidence of the normal to the boundary and the material symmetry axis direction that are 
represented by the formulas (3.82) in [5]. If a= > then we get the case when the symmetry axis of the materi- 


al is tangent to the boundary surface, and the formulas (21) coincide with the expressions (3.89) - (3.91) in [5]. 
Knowing the expressions for matrices C , one can construct solutions to various problems using the research 
and methodological apparatus described in [5]. 
Conclusions. The present research and methodological apparatus were successfully used for determining the 
stress-strain state in the layered structures made of transversely isotropic materials, and for analyzing the diagnostics 
results of the state of flat-layered and layered cylindrical structures in operation. 
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Hlouckoii rocyyqapcTBeHHbI TexHu4ecKHH yHuBepcurTerT, r. Pocros-Ha-Jlony, Poccniickaa Degepauna 


*5 }Ookuptii (beqepabHEI yHuBepcurer, r. Poctos-Ha-J[ony, Poccniicxaa Dexepanua 


Introduction. The investigation of the contact problems for 
cylindrical bodies is urgent due to the engineering contact 
strength analysis on shafts, cores and pipe-lines. In the present 
paper, a new contact problem of elastostatics on the interaction 
between a rigid band and an infinite two-layered cylinder, 
which consists of an internal continuous cylinder and an outer 
hollow one, with a frictionless contact between the cylinders, 
is studied. The outer cylindrical band of finite length is press 
fitted. By using a Fourier integral transformation, the problem 
is reduced to an integral equation with respect to the unknown 
contact pressure. 

Materials and Methods. Different combinations of linearly 
elastic materials of the composite cylinder are considered. 
Asymptotics of the symbol function of the integral equation 
kernel at zero and infinity is analyzed. This plays an important 
role for the application of the analytical solution methods. A 
key dimensionless geometric parameter is introduced, and a 
singular asymptotic technique is employed to solve the inte- 
gral equation. 

Research Results. On the basis of the symbol function proper- 
ties, a special easily factorable approximation being applicable 
in a wide variation range of the problem parameters is sug- 
gested. The Monte-Carlo method is used to determine the 
approximation parameters. The asymptotic formulas are de- 
rived both for the contact pressure, and for its integral charac- 
teristic. Calculations are made for different materials and for 
various relative thickness of the cylindrical layer including 
thin-walled layers. 


Discussion and Conclusions. The asymptotic solutions are 
effective for relatively wide bands when the contact zone 
length is bigger than the diameter of the composite cylinder. It 
is significant that the method is applicable also for those cases 


* The research is supported by RFFI Grant no.18-01-00017. 


Beedenue. AkTyalIbHOCTb MCCJI€¢OBaHHA KOHTAKTHBIX 3a]la4u 
WIA WMIMHIpHyeckux Tel obycNoBlIeHa HeOOXOAHMOCTBIO 
IIpOBeeHHA WHXKCHCPHBIX PacueTOB Ha KOHTAKTHY!O IIpou- 
HOCTb BaJIOB, CTepxHei u TpyOonpoBozoB. B HacTosulel pa- 
OoTe w3y4aeTCA HOBad KOHTAKTHad 3afjaya CTaTH4eCKOH Teo- 
PHM ympyrocTu o B3avMMOjelCTBUM %KecTKOTO OaHaxKa Cc Oec- 
KOHCYHBIM JIBYXCJIOHHBIM IMJIMHIPOM, COCTOALIHM H3 BHYT- 
PeHHero CIUIOWHOLO HW BHELMHerO MOOT WHIMHAPOB, MeXxKIy 
KOTOPBIMH BBITIOJIHAIOTCA yCIOBHA aqKOrO KOHTAKTA. 
HapyxHblit WuIHMHApHyeckHit OaHyax%K MocaxkeH Cc HaTATOM U 
wMeeT KoHeuHy!0 WWiMHYy. Tipu nomMouiM uHTerpasibHoro pe- 
oOpa3oBanua Dypbe 3aqaya CBOWMTCA K HHTerpayibHOMy 
YPaBHeCHHIO OTHOCHTCJIbHO HeH3BeCTHOrO KOHTAKTHOFO ]jaB- 
JIeHHA. 

Mamepuaavi u memoooi. PaccmatpuBarotca pa3Hbie KOMOHHa- 
WH JIMHeMHO-ypyrux MaTepHasioB COCTaBHOTO WWIMH pa. 
Uccnenyetca acuMmToTuka (yHKUHH-cMMBOsIa Apa HHTe- 
TpasIbHOrO ypaBHeHHA B Hye M OeCKOHeYHOCTH, urparoulad 
B@KHYIO POJIb VIA HCHOML30BAHHA AHAIHTHYeCCKUX MeTOJOB 
pemenua. Jia petieHua HHTerpasIbHOrO ypaBHeHHA BBOJMTCA 
OCHOBHOM Oe3pa3MepHbIi reoMeTpHuecKHH MapaMeTp Hu TIpH- 
MeHACTCA CHHTYJIAPHbIM ACHMNTOTHYCCKHM MeTO. 
Pesyasmamei ucciedoeanun. B cooTBeTCTBHH CO CBOHCTBaMH 
(PYHKUHH-CMUMBONIa TpesyoweHa cieWMabHat JIerKo (akTOo- 
pu3yemMas allmpokcuMalua 3TOM PyHKUMH, IpHroquad B WIH- 
pOKoM JManla30He H3MeHeHHA Hapamerpos 3ayqauH. IIpu no- 
mMoujH Metowa Moute-Kapmo paccuutaHbl lapamMeTpbl 9TO 
anmpoxcumayun. Tomryaenbl acuMnToTuueckHe (OpMysIbI Kak 
WI KOHTAKTHBIX JaBIeHHH, TaK HM [It HX MHTerpasibHOH xa- 
paktepuctuku. Pacuerbl cyemaHbi JJId pa3sHbIX MaTepHayioB U 
OTHOCHTCJIBHBIX TOJIMHH WMWJIMHpH4ecKoro COA, B TOM 4HC- 
Jie JJL1 TOHKOCTCHHBIX CJIOCB. 

O6cyarcdenue u 3akmo4uenua. WomyaeHuble acHMnTOTHYecKHe 
pelwenua 9pdekTHBHBI [Jit OTHOCHTeILHO WIMpoKHXx OaHya- 
%eli, Kora pa3sMep OONacTH KOHTaKTa TIpeBbilaeT THamMerp 
cocTaBHOoro WwWInHpa. BaxkHo, 4YTO HCHOUb3yeMbIl MeTOT 
ocTaeTca IIpHMeHHMBIM M JIA cilyyaeB, Kora BHeEWHHM WH- 


“ E-mail: pozharda@rambler.ru, zolotov.nikita.borisovich@gmail.com, ivan.sk6@gmail.com, pozharskaya.elizaveta@rambler.ru, 
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when the outer cylindrical layer is treated as a cylindrical 
shell. The asymptotic solutions can be recommended to engi- 
neers for the contact strength analysis of the elastic barrels 
with a flexible coating of another material. 


Keywords: elasticity theory, contact problems, composite 
cylinder, approximation, asymptotics. 


For citation: D.A. Pozharskii, et al. Contact problem for a 
two-layered cylinder. Vestnik of DSTU, 2018, vol. 18, no.3, 
pp. 265-270. https://doi.org/10.23947/1992-5980-2018-18-3- 





JIMHApH4ecKH COM MO%KHO paccCMaTpHBaTb Kak WMWIMHApH- 
yeckyio oOon0uKy. ACHMNTOTHYeCKHe pellieHHA MO2%KHO pe- 
KOMCHJOBaTb HH)KeHepaM JIA aHasiv3a KOHTaKTHOM Mpo4Ho- 
CTH yIIpyrux Jetaet WHIMHIpHyecko opMbI c yupyruM 
MOKpbITHeM 43 Wpyroro MaTepuaia. 


Ku1rouesBbie c10Ba: TeOpHaA yUPyrocTH, KOHTAaKTHBbIe 3aa4un, 
COCTaBHOM WWIMHp, alliipOKCHMallMa, AaCHMMTOTHKA. 


O6pa3ey dna yumupoeanua: Moxapckui, A. J. Kontaxtuaa 
3aqjaya Jia WByXxcioMHoro yHuHApa / J. A. Woxapcxnit [un 
ap]. — Bectaux Jjonckoro roc. TexH. yH-Ta. — 2018. — 


265-270 T.18, Ne3. — C. 265-270. https://doi.org/10.23947/1992- 


5980-2018-18-3-265-270 


Introduction. A dynamic contact problem for a prestressed elastic liquid-filled cylinder is studied in the paper 
[1]. Static contact problems for homogeneous elastic cylindrical bodies were considered in the papers [2-6] using regu- 
lar and singular asymptotic methods. It was found [4] that for cylindrical bodies, kernel symbols of integral equations of 
the contact problems are characterized by a more complex asymptotic behavior at zero and infinity than, for example, in 
contact problems for resilient strip. This required applying complicated approximations of the symbols easily factored 
by functions when using the singular asymptotic method. The approximation for hollow thin-walled cylinders is particu- 
larly complicated [6]. The suggested approximation [6] is advantageous even for the cases when a thin-walled elastic 
cylinder can be considered as a cylindrical shell [7]. The contact problem on the interaction of an elastic ring and a flex- 
ible cylinder was studied [8]. The elastic cylinder wear was analyzed in the paper [9]. This work objective is to obtain a 
solution to the contact problem for a composite two-layer elastic cylinder on the basis of a singular asymptotic technique 
and an effective approximation of the kernel symbol of the integral equation. 

Materials and Methods. In the cylindrical coordinates r, z (with axisymmetry), we consider an infinite elastic 
composite cylinder of outer radius R which consists of an inner solid cylinder of radius R/ <R with elastic parameters 
v,, G, (Poisson's ratio and shear modulus) and an outer cylindrical layer with elastic parameters v, G. Between the layer 
and the inner cylinder, the sliding binding conditions are met. We consider the contact problem on the interaction of the 
described composite cylinder and a rigid band in the domain |z|<a. For the given band preload 4, it is required to esti- 
mate contact pressures o. =—g(z) (|z|< a). Using the integral Fourier transform to solve the properly mixed (contact) 


boundary-value problem for Lame’s elastic equilibrium equations and introducing non-dimensional notations (the 
primes are further omitted) 








R 5 Zz q(6)d—-v) G, R 
A=-, F=-, C=-, WC)= , €=—, k= <1, 1 
. G o (C) G G 2 (1) 
we obtain the following integral equation for g(C): 
1 «oO 
J er( 53 Jae - m5 (<1), k()=|L(u)costut)du, (2) 
-1 0 
where the kernel symbol takes the form 
d d 
re eer (3) 


d= (A; Ags ~ As, As MA, As, am A, Ayy (Aj; Aug ~ A,, Ay) a 
—A,,4,, (Ay Aj, 7 A,,Ay) +r AAs, (Ay Ay; A,;Ay) + 

+A, As) (Ay Ay, = A,,A,)) ~ AAs, (A, Aj; = Ay, Ay, ) + 
+Ass Ase (434g, i Ay, Ag; A, Ag, mi A, Ag ) + 


+(A,) Ay, a Ay Ay (An Aes a AA.) + (AA, iz Ay Ay Ayer = A, Ag) + 
+(4,, Ay, = A,, Ay, M(Ai34o0 ~ A,,A;) + (Ay Ags iz A, Ay, MWA Aga ~ A,Ag) + 
+(A,, Ay G Ay, Ay, MWA Ags ~ AA, I, 

d, = (A; Ay ~ A, Ay, As, (As.Ags = A,; Age) mS A, Ass Ag) ] a 


—A,; As, [A, (Ay Ag, 7 A,,A,;) = Ay (AA, ~ AAs, )1. 
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d, = -(AAy, — 4,4 Ay; Ag) (Ag 65 — A55-Ase ) - 435 Ase | + 
+Ay5 Ase [Ay (Ay3-Ag4 — Aaa o3) — Aa (A364 — A)4-o3 ) I> 
A, =ul,-2(1-V)I,, Ay =-uK,-201-v)K,, 4, =ul,, A,=-uK,, 
A, =G—2v)ul, —(u2 + 4(1-v)h,, Ay =-B—2vuK, — 2 +4(1-v)K,, 
A,, =u(,—-ul,), A,, =—u(K, +uK,), 
A, =2(0l-v)I", Ay =201-v)Kt, 4,,=—40-v, I", 
A, =ukI{ -201-v)I", Ag =—ukKi -20-v)K*, Ay =ul?, Ay, =—uK*, 
Aog =ukEs -2(01-v, IP, Ag <ul’ 
A, =B—2w)ul* —Wk +40 WEYL, 
Ag =—3—2v)uKs —(W?k + 4(L—vyk) KY, 
Ag =uk Tf -w Ij, Ay =—uk'K; —u°K;, 
Ag = [3 —2v, ul? —<WkK+40—V EOI], Ag =—eu(I*h —ull), 
I,=1,@), K,=K,@, [, =1,uk), K;,=K,(uk), n=0, 1. 


Here, /,(u), K,(u) are modified Bessel functions. The dimensionless parameter 4 characterizes a relative width 
of the contact region. 
Function Z(u) at zero and infinity behaves as follows: 
(v-l(l+e-v, +ev,)+k’e, 
Av—-D[(v+ Dd +e-v, +ev,)+h'e,]’ 





lim L(w) = L(0)= &, =(v+ lv, —D-eWv—-Diy, +. (4) 


D 
Hj 4 26) (u—+0), D=1-2v. 
uu 


At k=0, value L(0) is the same as the known value for the homogeneous solid cylinder [4]. 

To solve the equation (2), we apply the singular asymptotic method [3,4] effective for sufficiently small values 
of 2. To apply the Wiener-Hopf technique [10], we used an easily factorable approximation of the function L (u) (3) by 
the expression 








L(u)= 


2 2 2: 2792 AB D 
Vu? +B of D Ss a s( } 6) 


UL 
— ; e 
u>+C V2 +10°) w+G L(0) 10° 


At calculations, two cases were taken: 1) iron inside, zinc outside (e=2.126, v=0.27, v\=0.28); 2) aluminum 





inside, zinc outside (g=0.779, v=0.27, v;=0.34). Table 1 gives the approximation parameters values (5), its relative error 





on the real axis 8 (%) calculated using the Monte Carlo method for different relative thicknesses of the outer layer k. 
Table 1 


Approximation parameters 





A B G 0 A B G 0 
k Iron inside zinc Aluminum inside zinc 
0.01 1.230 1.549 4.638 2.5 0.451 6.113 3.394 
0.03 1.291 1.563 4.042 2.5 0.888 5.977 3.029 
0.05 1.296 1.851 5.006 2.399 1.467 2.524 
0.07 1.258 2.057 6.376 0.720 7.945 2.512 
0.09 1.024 4.360 92.381 1.481 9.890 9.369 
0.099 3.540 24.178 5.347 1.178 5.535 235.946 10 
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Considering, as follows from (4), 
l-v, 


260) 


(6) 


then, in the case of small ¢, the approximation (5) for thin outer layers is to be complicated by increase in 
number of the parameters included into it. 
Research Results. As a result of applying the Wiener-Hopf technique, the principal term of the asymptotic solution to 


the integral equation (2) for small 4 can be constructed in the form of 





elf ie | 
q(G) | of 4 eof 7 i (|x|sD, (7) 
_ W(s) +13) oP es. 2 
cs) = TO” I(s) =-— J W(s—1)K,(10°2)dt, 
W(s) = OR eet +(4-1oca6,s) 


F-C 
VB-F 


Here, erf(x) is error function integral. 





O(F,s) = exp(—Fs)erf(,/(B — F)s)+ F erf (VBs ). 


For the integral characteristic of the solution 





1 
P= [a@ac (8) 
=i 
on the basis of the formulas (7), we obtain the expression 
Cee [z(2}+.(2)| =. ee [20 -9K,(10?2)dr, (9) 
5 'L(0) r Xr AL(0) Tt 








210)= (+ ate) erf(/Bs HS exp( ayl(4 1]T(4G,s), 


r1F.9)= (5 rare =) et) | ex(-80)] 
C )exp(—Fs) ——- 
[1 =) qs es aes 


Calculations show that the asymptotics (7), (9) error for A<1 does not exceed (5+0) %, where 0 is the ap- 
proximation (5) error. 








Table 2 shows the values of the integral characteristic PS” calculated from the formulas (9) at different values 


of k and 2. 
Table 2 


Values of PS" 


























A= 2 1 0.5 0.25 2 1 0.5 0.25 

k Iron inside zinc Aluminum inside zinc 
O.1[,]) 3.25 5.77 10.9 21.1 3.22 5.71 10.8 20.9 
0.03 3.37 6.02 11.4 22.2 3.20 5.65 10.7 20.7 
0.05 3.67 6.65 12.7 24.8 3.09 5.51 10.5 20.4 
0.07 4.22 7.78 15.0 29.4 2.96 5.33 10.2 19.9 
0.09 5.60 10.3 19.6 38.2 2.93 5.28 9.99 19.4 
0.099 6.91 12.3 23.2 44.9 2.99 5.31 9.94 19.2 
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Discussion and Conclusions. As Table 2 shows, with decrease in A, the integral characteristic of the contact 
pressures increases, which is associated with the extension of the contact area. For the case of a stronger material inside 
zinc (iron), contact pressures are higher than for aluminum inside zinc. 


With thinning the zinc layer around the iron (with increasing 4), the contact pressures increase essentially. This is 
hardly in evidence during thinning the layer of zinc around aluminum, since the modulus of longitudinal elasticity (and 
also shear modulus) of aluminum is slightly less than that of zinc. The asymptotics found can be recommended to 
engineers for analyzing the contact strength characteristics of the coated cylindrical parts. 
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Introduction. Reliability-Based Design Optimization (RBDO) 
model reduces the structural weight in uncritical regions, does 
not only provide an improved design but also a higher level of 
confidence in the design. 

Materials and Methods. The classical RBDO approach can be 
carried out in two separate spaces: the physical space and the 
normalized space. Since very many repeated researches are 
needed in the above two spaces, the computational time for 
such an optimization is a big problem. An efficient method 
called Optimum Safety Factor (OSF) method is developed and 
successfully put to use in several engineering applications. 
Research Results. A numerical application on a large scale 
problem under fatigue loading shows the efficiency of the 
developed RBDO method relative to the Deterministic Design 
Optimization (DDO). The efficiency of the OSF method is 
also extended to multiple failure modes to control several out- 
put parameters, such as structural volume and damage criteri- 
on. 

Discussion and Conclusions. The simplified implementation 
framework of the OSF strategy consists of a single optimiza- 
tion problem to evaluate the design point, and a direct evalua- 
tion of the optimum solution considering OSF formulations. It 
provides designers with efficient solutions that should be eco- 
nomic, satisfying a required reliability level with a reduced 


computing time. 


Keywords: Reliability-Based Design Optimization, Structural 
Reliability, Safety Factors 
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Beedenue. Mojenb, ocHOBaHHad Ha ONTHMH3alHH HayexKHO- 
ctu, (RBDO) yMenbuiaeT cTpyKTypHbIii Bec B HEKPHTH4eCKHX 
perHonax, oOecieuHBaeT He TOKO yilyYLIeHHy!O KOHCTpyK- 
I{MI0, HO H OosIee BLICOKHH ypOBCHb yBepeHHOCTH B WH3aliHe. 
Mamepuaavi u memoovi. Knaccuyeckuii noyxoq RBDO mo- 
%KeT ObITb BEITIOIHEH B ABYX OTCbHBIX IpOCTpaHctTBax: (pu- 
3H4eCCKOM IIpocTpaHcTBe H HOPMHMpoBaHHOM IIpocTpaHcTse. 
Ilockoubky B BbILIeyNOMAHYTBIX JBYX IIpoctpaHcTBax Tpedy- 
eTca OYeCHb MHOTO MNOBTOPHBIX McceqoBaHuli, pacuéTHoe 
BpeMA JIA TakKOM ONTHMM3AaIHH ABIIAeTCA OONbWOM mpobse- 
MoH. OddexTHBHEIM MeTON, Ha3bIBaemBIii Optimum Safety 
Factor (OSF), paspadotaH H ycnelIHO MmpHMeHseTcA K He- 
CKOJIBKHM HH)KeHEPHBIM MpHy10x7KeHHAM. 

Pesyiesmamel uccaedosanua. UncnenHoe mpysioxKeHHe 110 
KpynHomaciiTaOHolt 3ayaue pH ycTaslocTHOM 3arpy3Ke MOoKa- 
3bIBaeT 9(peKTHBHOCTB pa3paboTaHHoro MeTogwa RBDO ort- 
HOCHTeIbHO JeTePMHHHpOBaHHOM ONTHMM3alHH v3alina 
(DDO). Spdexrusyocts Metoya OSF tTawxKe pacipoctpans- 
eTCH Ha HE€CKOJIbKO peKHMOB OTKa30yYCTOMYMBOCTH IA 
yhipaBsIeHHA HECKOJIBKHMH BbIXOJHBIMH WapaMeTpaMH, TakH- 
MH KaK CTPyKTYpHbIi OObeM HM aTpHOYT NoBpexkAeHuA. 
O6cyorcdenue u 3akmo4uenus. YUpowleHHad cTpaTerua BHe- 
peHusa ctpyktyppi OSF cocTouT 43 eqHHCTBeHHOM 3aa4H 10 
ONTHMH3alHH OWeHKH NpoOeKTHOM TOUKH MH TIpAMOH OWeHKH 
ONTHMaJIbHOTO pellieHHa C y4ueTOM cocTaBoB OSF. Ox mpeyo- 
cTaBlIdeT pa3paOoTYHKaM 9PPeKTHBHbIe pelleHH1, KOTOPbIe 
TOJDKHEI ObITh IKOHOMHUYHBIMH, YHOBICTBOPAIOMIMMH Tpebye- 
MOMY YPOBHIO HajlexKHOCTH C COKpalljéHHBbIM pacuéTHBIM Bpe- 
MeHeM. 


KunoveBbie c10Ba: ONTHMH3al[Ma Ha OCHOBe HajexKHOCTH, 
CTpyKTypHad HajlexKHOCTS, (bakTOpbI Ge3onacHocTH 
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1. Introduction 

When Deterministic Design Optimization (DDO) methods are used, deterministic optimum designs are usually 
pushed to the design constraint boundary, leaving little or no room for tolerances (or uncertainties) in design, manufac- 
ture, and operating processes. So, deterministic optimum designs obtained without consideration of uncertainties may 
lead to unreliable designs, therefore calling for Reliability-Based Design Optimization (RBDO). An RBDO solution 
that reduces the structural weight in uncritical regions does not only provide an improved design but also a higher level 
of confidence in the design. The basic idea is to couple the reliability analysis with optimization problems. This cou- 
pling is a complex task, drawing on a high computing time and convergence stability, which seriously limits its applica- 
bility in real problems. To overcome these difficulties, several methods have been elaborated [1]. These methods can be 
classified into two categories: numerical and semi-numerical methods. An efficient numerical method called Hybrid 
Method is based on simultaneous solution to the reliability and the optimization problem. It has successfully reduced 
the computational time problem. The advantage of the hybrid method allows us to satisfy a required reliability level for 
different cases (static, dynamic, ...), but the vector of variables here contains both deterministic and random variables. 
To overcome both difficulties, an efficient semi-numerical method called Optimum Safety Factor (OSF) method has 
been proposed to solve problems in statics [2], and also an efficient alternative semi-numerical method called Safest 
Point Method (SP) has been proposed to solve problem in dynamics [3]. Recently, a Robust Hybrid Method (RHM) is 
also developed to overcome the HM difficulties in order to solve multiaxial fatigue damage analysis problems [4]. The 
RHM leads to robust solution comparing with the HM, but the computing time is still a big drawback. 


2. Reliability-Based Design Optimization 

2.1 Developments 

The computational cost of sequential RBDO approaches is much higher than the DDO procedure. Several devel- 
opments accelerated the use of the RBDO model. The Reliability Index Approach (RIA) and the Performance Measure 
Approach (PMA) have been proposed [5]. Next, the sequential optimization and reliability assessment (SORA) is de- 
veloped to improve the efficiency of probabilistic optimization [6]. The SORA method employs a single-loop strategy 
with a serial of cycles of deterministic optimization and reliability assessment. The major difficulty lies in the evalua- 
tion of the probabilistic constraints, which is prohibitively expensive and even diverges for many applications. It is clear 
that efforts were directed towards the development of efficient techniques to perform the reliability analysis. Here, the 
reliability index is computed iteratively that leads to an enormous amount of computer time in the whole design process. 


2.2 Basic RBDO formulations 


Traditionally, for the reliability-based design optimization procedure, two spaces are used: the physical space 
and the normalized space [7,8]. Therefore, the reliability-based design optimization is performed by nesting the follow- 
ing two problems: 


1. Optimization problem: 


min f (x) 
subject to g,(x)<0 , k=1...,.K (1) 
and B(x,u) =B, 


where f(x) is the objective function, g,(x) <0 are the associated constraints, B(x,u) is the reliability index of the struc- 
ture, and /;, is the target reliability. 


2. Reliability analysis: the reliability index B(x,u) is the minimum distance between the limit state function H(u) and the 
origin, see Figure 1b. This index is determined by solving the minimization problem: 


min d(u) 
me (2) 
subject to H(u)=0 
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where d(u) is the distance in the normalized random space, given by d = bys , and H(u) is the performance function 
(or limit state function) in the normalized space, defined such that H(u) <0 implies failure, see Figure 1b. In the physi- 
cal space, the image of H(u) is the limit state function G(x,y), see Figure la. Using the classical approach, the RBDO 
process is carried out in two spaces. That leads to a high computational time problem. Therefore, there is a strong need 
to develop efficient methods [9]. 


Physical Space uy, Normalized Space 












Failure 
Region 


Design H(u)=0 


Failure 
Region 


fi x “ G(x, y) =0 
: A 


a) b) 
Fig. 1. Physical and normalized spaces 


In the field of reliability-based design optimization, we distinguish between two types of variables: 


1. The optimization variables, which are deterministic variables to be adjusted with a view to optimizing the 
sizing; they represent the control parameters of the mechanical system (i.e. dimensions, materials, loads, etc.), 
and the probabilistic model (i.e. means and standard deviations of random variables); 

2. The random variables, which represent the uncertainties in the system. Each of these variables is identified 
by the type of distribution law and the associated parameters. These variables may be the geometric dimen- 
sions, the characteristics of the material or the external loads. 


3. Optimum Safety Factor (OSF) 


The Partial Safety Factors (PSF) presented in [10] use the calibration methods that need to propose some con- 
straints during the calibration process to increase the efficiency and the accuracy. The resulting solution when using 
PSF may not represent a global or even local optimum. It may satisfy the required reliability level because of the effi- 
ciency of the used optimization algorithm. An efficient OSF method essentially depends on the satisfaction of the opti- 
mality conditions of the reliability index problem (2). This method provides the designer at least with a local reliability- 
based optimum without additional computing cost. This method has been basically developed for a normal distribution 
case [11]. In this work, it is generalized to be applied to a single and multiple failure cases. 


3.1 Single Failure Mode (SFM) 
The SFM reliability problem can be written as: 


BY’ = min d(u,) = fu; tuy +...+u; st: H(u,,u,...u,)<0 (3) 


where d(u,) is the minimum distance between the design point and the optimal solution. And H (u, ) <0 represent the 


failure mode. The corresponding analytical formulation using OSF can be written as follows: 





i=l,....n (4) 








where the sign of + depends on the sign of the derivative, i.e., 


Spee Si and eae ai, i=l....,n 


i i 
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Formulation (4) provides different optimum values of the normalized variables at the design point and taking into ac- 
count a single failure mode. In [2], a similar formulation can be found for several distributions. 


3.2 Multiple Failure Mode (MFM) 


Using the same OSF developments for a single failure mode [2], the MFM problem can be written as: 


pw = min d(u,) = fu; tu; +..+u; st: A, (u,,u,.u,,) <0 (5) 


H, (u,) <0 represent the different failure modes. The corresponding analytical formulation using OSF can be written 


as follows: 


(6) 





Formulation (6) provides different optimum values of the normalized variables at the design point and taking into ac- 
count several failure modes. 


3.3, OSF algorithm 


The Optimum Safety Factor (OSF) algorithm can be easily implemented in three principal steps (Fig. 2). The 
first step is to determine the design point considering the most active constraint as a limit state function G(y). The opti- 
mization problem is to minimize the objective function subject to the limit state and the deterministic constraints. The 
resulting solution is termed the design point. The second step is to compute the safety factors using the equations (4) 
and (6). The third step is to calculate the optimal solution including the values of the safety factors in the computation 
of the values of the design variables and then determine the optimum design of the structure. 






Setinitialpoint: X)-¥o.- 







Evaluate Objective -function 


a 
z= 
3 
& 
5 
= 
s 
3 
eae | 








Optimization Formulation 


YES 


Compute OSF 


Reevaluate Functions 


Optimum Point 
Steps 2 & 3 





Fig. 2. OSF algorithm 


4. Numerical Application 


The objective of this application is to show the advantage of the RBDO by OSF relative to the DDO when deal- 
ing with SFM and MFM. The corresponding material properties are: Young’s modulus EF = 206.8GPa, Poisson’s ratio 
v= 0.29 and density p = 7820 kg/m’. The endurance limits for the reversed tension stress and torsion stress f-1 and t-1 
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stated after 2.106 cycles are equal to 252MPa and 182MPa, respectively. For more details about the fatigue data and 
methods of this example, the interested reader can refer to [12]. The length and the height of the studied plate are: 
L=0.14 m and H=0.1 m, respectively (Fig. 3a). In finite element analysis, the plate is supposed simply fixed on its four 
edges and is modelled by 32 eight-node square elements which produce no out-of-plane stress (Fig. 3b). Here, the 
thickness of each element 7. is considered as a random variable and its mean as a deterministic variable that leads to 32 


i 
5 §5 6 65 7 76 8 85 


x 10° 


deterministic variables and 32 random variables. 














a b 
Fig. 3: 
a - Dimensions of the perforated plate and 
b - Thickness distribution for resulting optimum design 


The objective of the DDO procedure is to minimize the volume subject to the fatigue damage constraints. Here, we 


consider a global safety factor 5S, =1.25 applied to the damage D,,,, and based on the engineering experience. The 


RBDO procedure cannot control not only the reliability level but several output parameters. In [2], the reliability level 
has been controlled considering a target reliability level, however, here we seek to control the other output parameters 
such as the structural volume and the damage criterion. 

The DDO and RBDO results in table (1) show that the DDO cannot provide the designer with a required reliabil- 
ity level while RBDO by OSF allow controlling the safety levels. According to the problems 3 and 5, the value of the 
global safety factor is applied to the upper damage limit D, =1 to be S, =1.25. This way the allowable damage will 


be: D,, =0.8. The standard deviations are considered as proportional of the mean values: 6, = 0.5m, , i=/,...,32. After 


having optimized the structure, the resulting DDO volume (Table 1) was found to be Vppo=/05.64cm’. The correspond- 


ing reliability index was found to be: B,,, =2.73. This resulting value does not belong to the standard structural engi- 
neering norm B,,5 =2.73¢ [3-4.25] . However, for the same optimum volume, the RBDO for multiple failure modes 


(using 5 and 6) provides the designer with a more reliable optimum structure with B = 3.64. Figures 4a and b show 


RBDO” 
the interval of the damage distribution of all structure thickness for DDO and RBDO in the same volume Vppo * Vrgpo’- 
The resulting damage distribution interval of RBDO by OSF [0.52,0.8] is better than the resulting one by DDO 
[0.4,0.8]. While Figures 4b and c show the interval of the damage distribution of all structure thickness for DDO and 
RBDO for the same maximum damage values Dppo=Dpgpo' =0.8. The RBDO by OSF provides the designer with an 


optimum structure with an increase 7% but more reliable by 45% relative to the resulting structure by DDO (f,,,,,,0 = 


3.78>B no = 2.73). Here, the resulting damage distribution interval of RBDO by OSF [0.56,0.84] is also better than the 


resulting one by DDO [0.4,0.8]. Thus, when obtaining a better damage distribution interval, the volume is reduced in 


the non-critical structural regions that leads to economic and reliable designs. 
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Table 1 
DDO and RBDO results for Bp” 






































































































































ee Design oe Solutions : 
Point RBDO DDO RBDO 
Volume 81.9770 105.53 105.64 112.70 
Dy max 0.9945 = 1 0.84 0.8 0.8 
aa 0 3.64 2.73 3.78 
T; 5.5784 8.914 7.7839 9.665 
T, 7.7249 11.477 9.2979 12.434 
T; 8.7796 8.340 10.1057 8.732 
Ty 8.7876 8.328 10.1055 8.718 
Ts 7.7435 11.571 9.2976 12.536 
Ts 5.604 9.001 7.7838 9.759 
T; 6.8081 12.024 8.6421 13.129 
Ts 6.0581 8.675 8.1291 9.265 
T» 5.3822 3.715 7.5599 3.738 
Tio 5.3866 3.745 7.5597 3.770 
Ty 6.0644 8.678 8.1290 9.266 
Tp 6.8267 12.063 8.6422 13.169 
T)3 7.605 10.503 9.2265 11.214 
Ti4 4.7228 4.081 7.1685 4.157 
Ti5 4.7194 4.061 7.1689 4.136 
T16 7.608 10.507 9.2268 11.218 
Ti7 7.6074 10.504 9.2267 11.214 
Tis 4.7191 4.057 7.1688 4.131 
Tio 4.7231 4.082 7.1685 4.159 
Tx 7.6056 10.509 9.2264 11.220 
To 6.8256 12.057 8.6421 13.162 
Tx 6.0637 8.675 8.1290 9.262 
T>3 5.3861 3.738 7.5597 3.762 
T24 5.3825 3.721 7.5599 3.745 
Ts 6.0588 8.680 8.1291 9.270 
Tx 6.8093 12.031 8.6421 13.137 
T27 5.6029 8.992 7.7838 9.749 
Ts 7.7425 11.557 9.2976 12.521 
Tx 8.7871 8.311 10.1055 8.699 
T30 8.7799 8.351 10.1057 8.744 
T31 7.7257 11.491 9.2979 12.450 
T32 5.5794 8.922 7.7839 9.674 
F The structural reliability level is improved because the OSF-based solution essentially depends on the sensitivity 
= study which determines the role of each parameter relative to the failure probability. For the computing time, the DDO 
z procedure needs to solve two sequential optimization problems (1) and (2) while the RBDO by OSF can realize the op- 
5 eration in only one single optimization problem. 
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Figure 4. Resulting damage distribution intervals by DDO and RBDO procedures 


Table 2 presents the different OSF results considering the linear and nonlinear distributions (normal, lognormal, 
uniform, Weilbull and Gumbel distributions), see [2], satisfying a required reliability level B = 3 . The required reliabil- 
ity level can be considered as a given data and the algorithm converges to the optimum design that verifies the require- 
ments (controllable designs: reliability, volume, damage ...). For the computing time consumption, for a single failure 
mode, when using DDO procedure, two optimization processes are used (the first is to find the optimum solution and 
the second is to find the design point). The problem becomes much more complex for multiple failure modes where 
each failure mode needs a separate optimization process to find the corresponding design point. However, the RBDO by 
OSF needs only a single optimization process to find the design point and next the optimum solution is analytically 
computed using OSF-SFM or OSF-MFM formulations. The RBDO by OSF is then carried out without additional com- 
puting time because it has a single variable vector that defines the design point. 


Table 2 
Linear and nonlinear RBDO result for required reliability index B =3 
























































Structural Design Optimum RBDO solutions 
Parameters Point Normal Lognormal Uniform Weibull Gumbel 
Volume 81.97 86.63 86.92 88.22 85.59 83.28 
Dinax 0.99 = 1 0.94 0.93 0.92 0.95 0. 98 
B 0 3.00 3.00 3.00 3.00 3.00 
Py 50% ~0.1% ~0.1% 20.1% =0.1% ~0.1% 
T, 5.5784 5.8679 5.8890 5.9694 5.7954 5.6604 
T> 7.7249 8.1947 8.2203 8.3520 8.0980 7.8541 
T; 8.7796 9.3409 9.3684 9.5255 9.2330 8.9342 
T, 8.7876 9.3488 9.3763 9.5334 9.2407 8.9428 
Ts 7.7435 8.2130 8.2387 8.3704 8.1160 7.8739 
T¢ 5.6040 5.8933 5.9145 5.9949 5.8204 5.6873 
T; 6.8081 7.1965 7.2204 7.3293 7.1097 6.9176 
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Ts 6.0581 6.3822 6.4046 6.4950 6.3039 6.1518 
To 5.3822 5.6633 5.6836 5.7618 5.5934 5.4579 
Tio 5.3866 5.6677 5.6880 5.7662 5.5978 5.4625 
Ty 6.0644 6.3885 6.4109 6.5013 6.3101 6.1584 
T 2 6.8267 7.2148 7.2389 7.3476 7.1277 6.9374 
T 13 7.6050 8.0591 8.0848 8.2121 7.9634 7.7337 
Ti4 4.7228 4.9483 4.9670 5.0288 4.8862 4.7866 
Tis 4.7194 4.9451 4.9638 5.0257 4.8831 4.7830 
T16 7.6080 8.0621 8.0877 8.2151 7.9663 7.7369 
T17 7.6074 8.0615 8.0872 8.2145 7.9657 7.7362 
Tis 4.7191 4.9448 4.9635 5.0254 4.8828 4.7827 
T19 4.7231 4.9486 4.9673 5.0291 4.8864 4.7869 
Tx 7.6056 8.0597 8.0854 8.2127 7.9640 7.7343 
To) 6.8256 7.2137 7.2378 7.3466 7.1267 6.9362 
Tx 6.0637 6.3878 6.4103 6.5007 6.3094 6.1577 
13 5.3861 5.6672 5.6875 5.7657 5.5973 5.4620 
Tx 5.3825 5.6636 5.6838 5.7620 5.5937 5.4583 
Ts 6.0588 6.3828 6.4052 6.4956 6.3045 6.1525 
Tx 6.8093 7.1976 7.2216 7.3304 7.1109 6.9189 
Tx 5.6029 5.8922 5.9135 5.9938 5.8193 5.6861 
Tg 7.7425 8.2121 8.2378 8.3695 8.1151 7.8729 
T9 8.7871 9.3483 9.3759 9.5329 9.2402 8.9422 
T30 8.7799 9.3412 9.3687 9.5257 9.2332 8.9345 
T3; T4251 8.1954 8.2210 8.3527 8.0987 7.8550 
T32 5.5794 5.8689 5.8899 5.9704 5.7964 5.6615 





























The OSF method is shown as a distinctive tool for RBDO problems. It shows the following advantages: 

e The obtained reliability-based optimum solutions should be more reliable than those obtained by DDO proce- 
dure for the same optimum volumes, 

e The OSF procedure needs only a single optimization process for the design point without additional computing 
time because it has a single variable vector that defines the design point while the DDO procedure needs two 
optimization processes. 

e Since the major difficulty lies in the evaluation of the probabilistic constraints, which is prohibitively expen- 
sive and even diverges for many applications, the OSF procedure provides the designer with an analytical 
evaluation with small computing time relative to DDO and other RBDO procedures. 

e = All reliability index evaluations for RBDO-MFM studies can be analytically carried out for different probabil- 
istic distributions. There is no need to optimization processes. 


7. Conclusions 


The RBDO using OSF has several advantages: small number of optimization variables, good convergence stabil- 
ity, small computing time, satisfaction of the required reliability levels and global optima and more economic solution 
for the same reliability index relative to the DDO process. The OSF-MFM formulation can be applied easily to different 
distribution laws. In fatigue analysis, the use of the classical method leads to an extremely high computing time and also 
local optima. 
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Mexannueckne cBoiicTBa CepBOBHTHBIX WJ1eHOK, POPMHPY!OWHXCA pH TpeHHH B BOAHBIX pacTBopax KapOoHo- 


CK 
BBIX KHCJIOT 
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34 HWM «Hanotexnosiorun 4 HaHOMatepnastbp> TamGoBcKoro rocyapcTBeHHOoro yHuBepcnteta um. I. P. Jlep»xasuna, TamOos, 


Poccuiickan Dexepayna 


Introduction. The effect of the organic component nature in the 
systematic series of monocarboxylic acids on the tribological 
characteristics of the brass-steel friction pair in aqueous solu- 
tions is described. Dependence of the mechanical-and-physical 
properties of the antifriction films formed during friction on the 
nature of the lubricating composition is investigated. The work 
objectives are to study the applicability of carboxylic acids as 
an antifriction lubricant component; to assess their effect on the 
mechanical properties of the servovite film formed under the 
brass — steel friction. 

Materials and Methods. Tribological studies of the brass-steel 
friction pair on the AE-5 end-type friction machine are carried 
out. Roughness parameters of the servovite film were deter- 
mined through the optical profilometry. The microgeometry and 
the object structure at the nanoscale were considered using 
atomic force microscopy. The mechanical characteristics of the 
antifriction film were investigated using the instrument 
nanoindentation. 

Research Results. Tribological characteristics of the brass-steel 
tribocoupling and mechanical-and-physical properties of the 
servovite film formed during friction in the “brass — aqueous 
solution of carboxylic acid — steel” system were studied. It is 
established that the friction factor reduces when increasing the 
hydrocarbon radical length. The dimensional effects are found 
in the mechanical and tribological properties of the servovite 
film formed on the surface of the friction interaction in the car- 
boxylic acids. 


Discussion and Conclusions. The study results show that the 





Beedenue. B padote moKa3aHo, KakHM OOpa30M Mpupoya op- 
TaHHYecKOH KOMIOHEHTBI B CHCTeEMaTHYeCKOM pay OAHOOC- 
HOBHBIX KapOOHOBBbIX KHCJIOT BIMAeT Ha TpHOouorMueckne 
XapaKTepHCTHKH Mapbl TpeHuaA «aTyHb — CTaJIb» B BOJHBIX 
pactBopax. V3y4eHa 3aBHCHMOCTb (p)H3HKO-MexaHHuecKHx 
CBOMCTB aHTH@PUKUMOHHBIX WIeHOK, POpMUpyIOWHxcaA Np 
TpeHHH, OT NpHpoyb! cMa304uHOH KoMNo3HuMH. Liem padorsi: 
W3Y4HTb BO3MOXKHOCTH MCMOIb30BaHHA KapOOHOBBIX KHCJIOT 
Kak aHTH(pHKUMOHHBIX KOMMOHEHTOB CMa304HOrO MaTepHa- 
la; OWCHHTb UX BIIMAHWe Ha MeXaHH4eCKHe CBOHMCTBa CepBO- 
BUTHOM IIeHKH, (OpMUpyrollelica IpH TpeHuu MaTyHH 110 
cTasiH. 

Mamepuaavi u memoovi. Wposeyzenbi TpuOomormuyeckue Hc- 
CleqOBaHHA Wapbl TpeHHA «aTyHb — CTaJIb»» Ha MallMHe 
TpeHua TopyeBoro Tuna AE-5. [lapamMetpbi mepoxoBaTocTu 
CepBOBHUTHOHM IWICHKM ONpeeAIMCb C MOMOLIbIO ONTHYeCKOH 
mpodunometpuu. Muxporeometpua u CTpyktTypa oObeKTa Ha 
HaHOYpOBHe MCCJIeqOBaIMCb C TIOMOLI[bIO aTOMHO-CHJIOBOM 
Mukpockonmu. Mexanwueckwe xapakTepHCTuKH aHTH@puK- 
IMOHHOM WICHKH W3Y4aJIM C MOMOLIbIO HHCTPyMeHTasIbHOTO 
HaHOHH]eCHTHpoBaHHa. 

Pezyiemamei uccaedosanua. Visyaenbi tTpudonormyeckue 
XapaKTepHiCTHKH TpHOoconpmKeHHA «aTyYHb — cCTalIb» HU 
(pu3MKO-MexaHHyeckHe XapaKTepHCTHKH CepBOBHTHOM MIeH- 
KH, (opMupytoljelica pH TpeHHu B CHcTeMe wiaTyHb — 
BOJHBIM pacTBOp KapOOHOBOH KHCIIOTbI — cTalIb». YcTaHoB- 
JI€HO, 4TO Ip yBeIMYeCHHM [MHI yrueBoOpoyHoro payH- 
Kala KooPPUUHeHT TpeHua cHwKaeTca. OOHapyxKeHbI pa3- 
MepHbiIe 9oP@ekTbl B MeXaHW4YeCKHX HM TpHOoMOrM4eCKHX 
CBOMCTBaxX CepBOBHTHOM MWIeHKU, (OpMupyroljelica Ha T0- 
BePXHOCTH (PpHKIUHOHHOTO B3aHMOeHCTBUA B BOJHBIX pac- 


TBOpax KapOOHOBBIX KUHCIJIOT. 
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friction interaction on the wearing surface in the aqueous solu- O6cyorcoenue u 3akiiouenua. Pe3ynbTaTbl MCCIeqOBaHHA T10- 


tions of carboxylic acids forms a nanostructured servovite film Ka3bIBalOT, 4YTO IPH (PHKUMOHHOM B3aMMOLelicTBMM Ha T10- 


which drops the friction factor. Its mechanical, physical and BEPXHOCTH TPeHHA B BOJHEIX pacTBOpax KapOOHOBEIX KUCIOT 


cbopmupyetca HaHOCTPYKTypHadA CepBOBHTHaA WICHKAa, Pe3KO 


tribological parameters depend on the composition of the model 
cHIpKaloljat  KooppuuneHt tTpenusa. Ee du3nKo- 


lubricating medium. It is determined that the local mechanical- 
MexaHWyeckve WM TpHOONOrMYecKHe HapaMeTpbI 3aBHCAT OT 


and-physical properties depend on the method of producing the cocTaBa MOmeNBHOi cMa3z04HOH cpempr. OnpeneneHo, 470 


servovite layer, the load and the size of the deformation zone. HOKasbHbIe :pu3UKO-MexaHMyeckue CBOliCTBA 3aBHICAT OT CIIO- 
The results obtained can be used in the development of lubri- co6a TlOJIy4eHHA CePBOBHTHOFO CJIOA, Harpy3KH WM pa3Mepa 
cants. 30HEI ecbopmMuposanusa. [lomyyeHHble pe3yiIbTaTbl MOryT 


ObITbh HCHOJIb30BaHbI Tipu pa3paboTKe CMa304UHBIX MaTepHa- 


JIOB. 
Keywords: friction factor, selective transfer, servovite film, Kooyepbie c0Ba: KosppuyHesT TpeHua, u30upaTesbHBIli 
dimensional effects, surface roughness. TlepeHOC, CepBOBHTHat TIeHKa, pasMepHble 9@exThl, Wepo- 


XOBATOCTb MNOBCpXHOCTH. 
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Introduction. A rapid development of modern machine building increases opportunities of constructing ma- 
chines that operate under the extreme conditions. Thus, the motivation of the life growth of the friction units rises with 
the demanding requirements for the mechanical properties of materials. 

The friction conditions with high energy density impose a number of specific requirements on the machines 
and equipment in terms of the tribological characteristics, reliability and service life. The critical working conditions of 
tribo-mating surfaces affect the flow of force, loss of energy, and, eventually, dynamic behavior of the complete mech- 
anism. Space distributions of topographic irregularities, as a rule, negatively affect the operational capabilities of ma- 
chines and mechanisms. In other words, the load-bearing capacity of friction units hinges on the basic parameters of the 
surface state. They affect the interaction conditions in the contact area, adhesion, and deformation, thereby generating 
vibrations during the friction process [1] and causing uneven wear. 

In this case, the greatest antifriction efficiency is demonstrated by the lubricants [2-6] containing metal addi- 
tives, such as copper, aluminum, silver. These materials form protective metal plaque shells on the tribo-mating surfac- 
es in the process of friction. Thus, a low friction factor and a medium low wear are provided in a wide range of contact 
pressures and slip velocities [2, 3]. The process of metal plaque filming of the surface layers in-service enables to “heal” 
surface defects (cracks, pores, accumulated fatigue damage). This essentially increases the antifriction characteristics 
and improves the mechanical-and-physical properties, which is of prime importance for the reliable operation of the 
friction units, and the lifetime extension. 

It is notable that the formation of transfer films is possible under friction in aqueous-alcoholic media with no 
visible total wear of the mating surfaces [7, 8]. In the engineering practice, the use of selective transfer gives a real op- 
portunity for designing long-life and effective friction units in the movable mating of parts and components of machines 
and mechanisms [8]. 

In this context, it is interesting to consider the effect of the lubricating medium composition on the mechanical 
and tribological parameters of the servovite film forming in the “brass-steel” tribo-pair (aqueous solutions of carboxylic 
acids are used as a lubricant medium). 

Materials and Methods. Evolution of the friction factor of the “brass 59 — aqueous solution of carboxylic ac- 
id — steel 40X” system was investigated in the “Hybrid Functional Graphene-Based Materials” laboratory of the Re- 
search and Education Center (REC) “Materials” on the AE-5 end-type friction machine. The experiments were carried 


out under the following conditions: rotation speed was 180 rpm; axleload was 98 N; the test time was 10 hours. The 
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monocarboxylic acids of the limiting series with the general formula of R— COOH (R = C,A,.;) were used as an or- 
ganic component of the lubricant composition. 

The mechanical characteristics (hardness — H and Young modulus — £) of servovite films were determined by 
the instrument indentation method [9, 10]. Thereat, we used: 
—Nanotest Platform 3 (Micromaterials, UK) nanoindentation tester with the attached function unit that enables to im- 
pose loads from 0.01 mN to 500 mN; 
— Nanotest Platform 3 (Micromaterials, UK) tribological nanoindentation tester, a multifunctional system for mechani- 


cal testing of materials through the dynamic nanoindentation. 


The nanoindentation technique makes it possible to study the package of mechanical properties of light boundary 
layers of solids and films up to a few tens of nm thick [9-15]. The mechanical characteristics were studied using a dia- 
mond Berkovich indenter. Taking into account the specific thickness of the investigated films (from several hundred 
nanometers to units of micrometers), the operational mechanical parameters were measured at the nanoscale at indenta- 
tion depths from 20 nm to the micron units. At the nanoscale, the friction coefficient was investigated on the Tri- 
boindenter TI-950 with simultaneous application of normal and lateral loads to the indenter [11, 16, 17]. The normal 
and lateral components of the forces (Fy and F’,) and displacements (Ay and h;) implemented during the replication of 
the tribo contacts by the indenter were recorded continuously. The recorded data were analyzed, which permitted to 
describe the friction and wear processes. The friction coefficient k,,, was an evaluation item. 

The Contour GT-K1 optical profilometer with Vision 64 analytical software installed in the REC “Materials” 
was used to determine the servovite film thickness and roughness parameters. The measurements were taken using ver- 
tical scanning interferometry (VS/) at the scan rate of 0.1 um/s with the repeatability of RVS 0.01 nm. 

The servovite film surface topography was studied using the PHYWE Compact atomic-force microscope 
(AFM) in a semi-contact mode with a single-crystal silicon aluminum-coated probe. 

Research Results. A servovite film, spontaneously arising under the friction interaction in the glycerin solu- 
tion in the copper-steel alloy, is studied fairly well [8, 18]. It is common knowledge that it is formed on the friction face 
under the topographically unequal conditions. Thus, in different places of frictional contact, the properties of the servo- 
vite film are different. Its thickness is less than 2 um, and the mechanical, physical and tribological properties differ 
from those of base compact copper. 

It should be mentioned that one of the glycerin tribooxidation [19] products is carboxylic acid. In that con- 
text, it is interesting to consider the effect of the lubricating medium composition on the mechanical and tribological 
properties of the “brass-steel” tribo-pair under friction in the systematic series of monocarboxylic acids. 

Long-term evolutionary tribological studies of the “brass-steel” friction pair in aqueous solutions of carboxylic 
acids have revealed the following dependence of the tribological characteristics on the lubricating medium composition: 
the length of the hydrocarbon chain of the acid radical increases from Cy to Cs, and the friction factor decreases to 
0.007. 

At that, the lowest value of the friction factor, which is characteristic for the systems implementing the wear- 
lessness effect [20, 21], is obtained in aqueous solutions of the valeric and caproic acids. At the same time, the “brass- 
steel” tribo-pair wear is reduced by 25 times. As a result of the selective transfer under the tribo-mating surfaces fric- 
tion, a copper film with different roughness and surface coating density is formed. 

The friction surface was scanned under transition conditions in the series of the “formic - acetic - propionic - 
oil - valeric — caproic” acids. At that, optical profilometry has revealed, first, a decrease in the servovite film roughness. 
Secondly, the dependence of roughness on the initial topography of the of the test steel disc surface (its R, is equal to 


118 nm) and on the lubricant composition was determined [22, 23] (Fig. 1). 
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Fig. 1. 2D visualization, 3D visualization and surface profile from steel (a), (c) and brass-59 (5), (d) after friction in aqueous solu- 
tions of acetic (a), (b) and caproic (c), (d) acids 


Roughness R, at the base length (Z = 800 um) of the servovite film decreases, for example, in the aqueous ca- 
proic acid, up to 69 nm. (For comparison: this index comes up to 580 nm in an aqueous solution of formic acid and to 
401 nm — in an acetic acid solution.) The difference is due to the corrosive activity of the medium. 

In the formic and acetic acid solutions, scores of scratches, irregularities and pores occur on the tribo-mating 
surfaces, as well as areas with a copper film formed in the island growth mode. In the cross-section of the friction sur- 
face, deep longitudinal grooves are revealed, which goes to show the abrasive wear pattern of the tribocoupling [23, 
24). 

A more detailed scanning of the surface using the AFM method and image visualization discovers considerable 
damage under friction in the formic and acetic acid solutions. On a 3D image, they are visualized as macroscopic dark 
bands corresponding to the areas of intense frictional stress accompanied by strong abrasive wear of tribo-combinations 


(Fig. 2). 
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Fig. 2. Topography of servovite film surface obtained under friction in aqueous solution of acetic (a) and caproic (5) acid 


A comparative analysis of the friction surface morphology upon transition to the acids with a longer hydrocar- 
bon radical, suggests the wear pattern transition, and the structural modification of a thin subsurface film of tribocou- 
pling. 

It follows from the mass transfer of the components of the contacting bodies, as well as from the adsorption of 
copper nanoclusters from the working medium to the counterbody as a result of tribo-electrochemical processes in the 
friction zone during the implementation of the selective transfer. At this, a copper film of fine-grained copper nanoclus- 
ters is formed on the steel surface under the friction interaction in the aqueous solutions of caproic and valeric acids. 
The layer formed on the surface is fairly dense, with a small spread of particles in size (see Fig. 1, 2). 

By now, a significant amount of experimental data has been accumulated [25-29] implying a major change in 
the mechanical properties with reducing the characteristic dimensions of the elements of the object structure of less than 
100 nm. This permits to assume a change in the mechanical parameters of the servovite films as compared to those for 
copper in the volume [30]. To identify this fact, an instrument indentation was carried out. 

The hardness values of H and Young's modulus £ were determined by the Oliver-Pharr method [10] from the 
characteristic P-h diagrams (Fig. 3) in accordance with the standards [9]. 
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Fig. 3. Representative P-h diagrams for servovite film obtained through friction of “brass-steel” pair in aqueous solution of valeric 
acid (a), caproic acid (b) 


The undertaken studies showed an increase in the values of H and E films compared to the similar characteris- 
tics of the copper sample. 

The increase in the hardness H and the elasticity modulus £ of the film formed on the friction surface can be 
ascribed to the change in the plastic deformation mechanism. A high vacancy concentration at the grain boundaries dur- 
ing friction interaction promotes grain-boundary sliding of building blocks relative to each other. In this case, the non- 
dislocation plasticity mechanism is implemented [31-34], which results in improving strength characteristics of the cop- 
per surface layer, and in increasing ductility during the transition in the series of the “formic - acetic - propionic - oil - 


valeric — caproic” acids. 
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Besides, the analysis of the dependence of the H and E values of the servovite film on the depth of the h, plas- 
tic impression shows that the hardness is affected by the scale factor. The Young's modulus values remain practically 
constant up to the depths associated with the servovite layer thickness. The films formed through friction in valeric acid 


are of much evidence (Fig. 4). 
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Fig. 4. Dependence of hardness (H) and Young's modulus (£) on the depth of plastic impression (4,) (shown in semilog co- 
ordinates) when indenting servovite film obtained in aqueous solution of valeric acid (a), caproic acid (b) 


The following dimensional effect is identified: 

- H value increases with growth of /, in the shallow depth-region; 
- with further growth of h,, H falls dramatically. 

For other samples, the dimensional effects are feebly-marked. Thus, for example, for the film formed in capro- 
ic acid, the values of H and E remain practically constant within the full tested range of h, (from 23 nm to 1.1 pm). 
Then, with increasing h,, the values of H and F decrease. Obviously, this is due to the increase in the deformation zone 
size: the properties of the transition layer and the substrate material start to affect H and E. And, in the region of deep 
depths of the indentation (4, > 1.1 um), the values of H and E decrease (see Fig. 4). 

The dimensional effect [35, 36] is also found under studying the dependence of the friction factor on the ap- 
plied load in the nano- and microscale [11, 16, 17] (Fig. 5). 
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Fig. 5. Dependence of friction factor & of servovite film on load in nanoscale 


First, it is expressed in a sharp drop in the friction coefficient with an increase in the normal load, up to some 
value of Fyonm. Then, the friction coefficient increases smoothly with increasing Fy. So, the dependences k = f (Fy) have 
sharp minimum if some value of Fryonm is reached [11, 16, 17]. 

The experimental data obtained can be compared to the theoretical dependence of the friction factor on the 
normal load. Hence, it is becoming apparent that they are in good qualitative agreement to all the film samples under 
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study. The Fonm value is affected by the type of carboxylic acid used as a lubricant. For example, for caproic acid, 
FNonm = 30 mN, and for valerian acid, Fyonm = 100 mN. 

Discussion and Conclusions. The mechanical properties of servovite films formed in the “brass-steel” friction 
pair are studied. The dimensional effects in the mechanical (hardness) and tribological (friction factor) properties are 
identified in the nano- and microscale. It is shown that the properties studied and the nature of their changes depend on 
the type of carboxylic acid used. The decrease in the average size of crystallites in the composition of film nanoparticles 
leads to an increase in its strength and the development of superplasticity under tension and shear. 

The results obtained in the paper permit to draw the following conclusions. 

1. During the friction interaction of the “brass-steel” pair in aqueous solutions of carboxylic acids, an antifric- 
tion servovite film is formed on the friction surfaces, which contributes to a sharp decrease in the friction factor. 

2. The servovite film, which is formed under the tribo-interaction of the “brass-steel” pair, is nanostructured. 

3. The tribological parameters of the “brass — aqueous carboxylic acid — steel” system depend on the length of 
the hydrocarbon radical of carboxylic acid used as a lubricant. 

4. The lubricating medium composition affects the roughness parameters of the servovite film: in the transition 
from formic to caproic acid, the film surface roughness is reduced to one-quarter. 

5. Local physical and mechanical properties depend on the conditions for obtaining the servovite layer, the 
load, and the size of the deformation zone (the indentation depth). 

6. The detected dimensional effects (dependence of the determined values on the indentation depth) verify the 
conclusion on the nanostructured character of the servovite film formed in the “brass — aqueous carboxylic acid — steel” 
tribosystem. 
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Introduction. Power parameters of the “rolling — equal- 
channel angular pressing” (ECA) combined process are stud- 
ied. The work objective is to determine forces of rolling and 
pressing in the deformation by the combined method. 
Materials and Methods. The die strength calculation and the 
experiment on deformation of AISI 6063 aluminum samples 
were carried out. During the experiment, the force values were 
recorded using a strain-gauge station. 

Research Results. The strength analysis results show that this 
die design is suitable for creating an experimental stand of the 
“rolling — ECA-pressing” combined process, since the calcu- 
lated safety margin is sufficient to implement the pressing 
under extreme conditions. The rolling forces at all stages of 
the deformation exceed the corresponding pressing forces, 
which is a necessary condition for the combined process. 
Discussion and Conclusions. The obtained results can be used 
in the design of experimental stands that implement the inves- 
tigated combined process. Herewith, the given strain-gauge 
method for studying strength characteristics is suitable for the 


case of calibrated rolls. 


Keywords: combined process, rolling, equal channel angular 
pressing, force, strength analysis, tensometry. 
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sek 


PaOota BBIMOIHeHa B paMKax HHHMUMaTHBHOH HUMP. 


Beedenue. CraTba NOcBaleHa HCCeOBaHHI0 IHEPYOCHJIOBbIX 
TlapaMeTpOB COBMeLICHHOTO Tpolecca «NpokaTKa — paBHO- 
KaHaJIbHoe yrsloBoe mpeccoBanue (PKY)». Hemp paborrr — 
onpeyeseHve BO3HHKarOlIMX yCHIMM WpokaTKH HM TpeccoBa- 
HHA IPH JePopMupoBaHHi COBMELICHHBIM CIIOcoOoM. 
Mamepuaavi u memooot. bait WpoBeyeH MpOuHOcTHOH pacer 
MaTpHUbI HW 9KCIepHMeHT 10 FesopMupoBaHHo OOpa3sOB 13 
amzoMuHuaA MapKH AISI 6063. B xoye onbiros 3HaueHHA ycu- 
MM (bUKCHpOBaJIMCh C MOMOLIbIO TeEH3OMCTPH4eCKOH cTaH- 
HH. 

Pesyiomamol ucciedoeanua. Pe3yibTaTbl IpOUHOCTHOrO aHa- 
JIM3a TO3BOJIAIOT YIBep%*KaTb, YTO Jawad KOHCTPyKIMA MaT- 
PHUbI UpHroyHa [Wid CO31aHHA 9KCIeCpHMCHTAaIbHOH yCTaHOB- 
KH COBMeINeHHOrO mporecca 
IIpe€CCOBaHHe>>, NOCKOJIBKY PaCCUMTaHHOrO 3allaca NPOUHOCTH 


«mpokatka — PKY- 


BHOIHE MOCTATOYHO WIA peallM3allM4H MpeccoBaHHa Ha Tipe- 
WesIbHBIX YCIOBHAX. YCHJIMA NPOKaTKH Ha BCeX 9Tallax ]e- 
(opMupoBaHHaA NIpeBbILWIalOT COOTBETCTBYIOWIMe YCHIIMA TIpec- 
COBaHMA, YTO ABJIAeTCA HEOOXOAMMBIM YCJIOBHEM COBMELICH- 
HOrO TIpowecca. 

O6cyarcoenue u 3axmouenua. HomyaeHuble pe3ybTaTbl MOTyT 
OBIT HCHOJL30BaHbI IIpH NpOeKTHPOBaHH 9KCIepHMeHTAasIb- 
HBIX YCTAHOBOK, PeCaJIH3YIOWIHX MCCIIeqyeMbIM COBMECLICHHBIM 
mpouecc. [Ip sTomM TeH30MeTpHuecKH MeTO] MCCIeqOBAaHHA 
IIPO4YHOCTHBIX XapaKTepHCTHK TIpHroyeH MW Mp UCIOsIb30Ba- 
HMM KasIMOpOBaHHBIX BaJIKOB. 


K.unoveBbie C10Ba: COBMELICHHBIM Mpouecc, MpoKaTKa, paB- 
HOKaHaJIbHOe YIIOBOe MpeccoBaHve, ycusIMe, MpoOuHOCTHOH 


aHasIH3, TCH3OMeTpHA. 
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Introduction. Improvement of the deformable metal quality is one of the major tasks when creating a new 
metal forming process (MFP) or improving the known one. No matter which type (rolling or blacksmithing) the process 
belongs to, it is recommended to apply a deformation scheme implementing a pure shear. In this case, a minimum 
amount of energy is expended on the metal deformation; maximum and uniform workpiece processing through the sec- 
tion is performed; and internal defects are closed. 

Over the past two decades, forming techniques have been actively developed. They enable to obtain ultrafine- 
textured blanks. Their analysis shows that they are based on the principles of implementation of shear or alternating 
deformations. Processes that are a combination of these principles are worthy of special mention. These processes allow 
for the implementation of an unconventional forming operation which is called “severe plastic deformation” (SPD). 
SPD implements such methods as torsion under high pressure [1-3], equal-channel angular (ECA) pressing [4-8], screw 
extrusion [9-12] and others. Special attention should be given to an echelon ECA-die (or ECA-die with parallel chan- 
nels) [13-14]. It enables to realize a shear deformation when the workpiece passes through its channels and concurrently 
two alternating deformation zones, provided that the input and output channels are codirectional. In comparison with 
the usual angular die, the echelon ECA-die is energy-saving, since it allows for implementing a greater degree of de- 
formation in one pass with the equal force. 

The MFP methods discussed above belong to the discrete press forming techniques. These technologies for the 
molded articles production have a number of shortcomings. The principle ones are due to: 

- press operation discreteness (discontinuity); 

- occurrence of reactive frictional forces on the metal-to-container contact. 

Thus, due to uneven deformation and high energy intensity of the press operation, the length of the pressed 
products is limited, and their quality is reduced. The application of the continuous pressing scheme enables to eliminate 
these drawbacks. This method differs fundamentally from the previous ones in that active friction forces are used to 
deform and extrude the workpiece through the die hole. The length of the workpiece is thus not limited. 

In recent years, the so-called “combined” MFP processes, which are a combination of two or more defor- 
mation processes, are being developed. 

Crucially, the most important feature of the combined MFP process is that shortcomings of the components of 
individual processes are often reduced or completely eliminated under its implementation. 

Besides, more and more attention has been recently paid to the energy-saving technologies based on the use of 
active friction forces for deformation. 

The authors of the study have developed a combined technique of rolling and pressing in an equal-channel 
echelon die, which, in comparison to the conventional pressing in an equal-channel echelon die, removes the limitations 
on the dimensions of the initial blanks [15]. It seems to be one of the most advanced and understudied on the practical 
level rolling and pressing processes, and it in many ways surpasses the known forming methods. The principle of this 
straining method is as follows. The workpiece, pre-heated up to the temperature of the deformation onset, is fed to the 
cast rolls. The first pair of rolls, due to the contact friction forces, bites the workpiece and pushes it on exit through the 
channels of an equal-channel echelon die. After the workpiece completely leaves the die channels, it is bitten by the 
second pair of rolls, which completely pull the workpiece out of the die also by contact friction forces. Thus, in this 
case, pressing the workpieces is provided by the friction contact forces that arise on the surface of the metal-to-rotating 
rollers contact. Previously, various studies of this combined process have been carried out — both theoretical [16-19] and 
laboratory ones (on a single-roll machine) [20-23]. 

This paper presents the results of the power analysis performed during the development of a pilot plant with 
two pairs of rolls (Fig. 1). 
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Fig. 1. Pilot plant for “rolling - ECA-pressing” combined process implementation 


Strength analysis of the ECA-die. When implementing the combined process of the “rolling — ECA- 
pressing”, a pilot plant created on the basis of the laboratory rolling mills DUO-200 and DUO-250 was used. Their 
power characteristics enable to cold-form aluminum and copper alloys. 

When designing a new technology of deformation at the pilot stage, strength analysis is essential. It enables to 
determine the tool durability and rigidity, and to estimate its strength margin which will give information on its service 
life. Theoretical studies of the “rolling -— ECA-pressing” combined process were performed using computer simulation 
in the DEFORM program. As a result, it was established [17, 18] that the most loaded part is the pressing die. This is 
explained by the fact that the most rigid straining conditions occur in the joint mouth zone. Thuswise, a strength analy- 
sis was carried out to test the selected die design validity. For this purpose, a model of a joint die mouth with 150 degree 
angle represented in [20], was constructed. 

When implementing the combined process, it is appropriate to select the workpiece width in such a way that 
under forming in the die, the metal does not contact with its sidewalls. This will significantly reduce the die total back 
force. However, one cannot completely exclude such a case. Thus, an ECA-pressing model, in which the workpiece 
height and width correspond to the height and width of the die mouth, was created in the DEFORM software. This case 
of deformation is limiting and is almost impracticable. Therefore, if the die is hard enough, it will be suited for a com- 
bined process, the conditions of which are much “softer”. In particular, the workpiece height and width is somewhat 
less than the height and width of the die mouth, which is a necessary condition for the metal passage through the mouth. 

The initial blank dimensions are 10 x 40 x 90 mm. The conditions and assumptions used in modeling: 

- initial material is isotropic, it does not have any stress or strain; 

- workpiece is divided into 24,000 finite elements, average length of the element edge is 1.5 mm; 

- initial temperature of the workpiece is 700 °C, besides, strain warming-up, and heat transfer from workpiece 
to tool and into the environment are considered; 

- instrument is taken as an absolutely rigid body; 

- dummy blank is taken as elastoplastic; 

- workpiece material is AISI 1045, corresponding to steel 45; 

- friction factor at the workpiece-die contact is 0.1; 

- punch speed is 10 mm / sec. 

As a result of the calculation, the deformation force on all die elements was determined to find the most strained struc- 
tural members. The greatest forces occur in the direction of the normal stresses action (Fig. 2). 
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Fig. 2. ECA-pressing force curves for sidewalls (a) and central segments (5) 
292 Fig. 2 shows that the most strained die elements are the central segments. At deformation, the force up to 260 


KN occurs on them, while it does not exceed 75 kN on the sidewalls. Based on these data, the stress analysis of the tool 
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was carried out in the DEFORM program using a dedicated module. This analysis approach implies that absolutely rig- 
id model elements become elastic bodies taking stress. Interpolation of operating forces on the tool is carried out from 
the workpiece. Then, fixation points are defined so that the elastic bodies do not scatter under the action of forces. In 
our case, the fixation points are ring openings for holding elements. In addition, in this calculation, it is required to 
specify the tool material. 5XB2C instrument steel was used to make the die. It is used for manufacturing: 

- complex dies operating with high impact loads, 

- knives for cold metal cutting, 

- thread-rolling dies, 

- punches and swaging dies for cold operation. 

As a result, the following data were obtained (Fig. 3). 
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Fig. 3. Strength analysis results: stress on sidewalls (a); tension in central segments (5); linear displacement (c) 


Fig. 3 shows that the most loaded zones in the central segments are surface areas near the die mouth joints and 
arc zones around the first two pairs of openings. Here, voltage up to 250 MPa occurs. 

In the areas between openings, the voltage is within 125 + 150 MPa (see Fig. 3, a). 

On the sidewalls, peak voltage occurs in the contact zone with the deformable metal, and it is equal to 95 + 100 
MPa (see Fig. 3, 5). 

The results of calculating the linear displacements of structural elements under the action of elastic defor- 
mation were considered [24] in the specialized calculation APM library to study the structure rigidity. The investigation 
of elastic deformation shows its greatest value in the surface region near the die mouth joints. Two facts should be men- 
tioned. The first: operating stresses on all die elements are much lower than the 5XB2C steel yield limit. The second: 
peak value of linear displacement reaches 0.019 mm (Fig. 3, c), which is negligible and allows us to speak about the 
increased rigidity of the proposed structure. In light of this, we can conclude that this matrix design is suitable for creat- 
ing a pilot plant of the “rolling - ECA-pressing” combined process, since the estimated safety margin is sufficient for 
ECA-pressing implementation under extreme conditions. 

Laboratory experiments to study forcing. As a result of the strength analysis, the location of the strain gage 
sensors in the matrix was determined for measuring the deformation force. The mounting spots are at a distance of 75 
mm from the leading edge of the central segment where peak linear displacements occur. 

The following equipment was used to measure the emerging deformation forces: 

- ZET-017-T8 strain-gauge station (“ETSM” CJSC, Russia); 

- force sensors with TKFO1-2-200 strain gages (“ETSM” CJSC, Russia); 

- laptop computer to control the strain gage and to record the signal. 

The strain-gauge station manufactured by “ETSM” CJSC (Moscow) is designed to measure power and other pa- 
rameters on several channels with a time resolution of writeenable up to 20 kHz [25]. 

On investigating the power parameters of the “rolling — ECA-pressing” combined process, it is planned to 
study the deformation forces at the stages of rolling and molding in dies. This calls for installation of strain gauges on 
the die and on the measuring cells of the roll mill. The principle of tensometry is that during the deformation process, 
the rolling and pressing forces affect the rolls and the die. As a result, the rolls and the die are bent. Strain gages at- 
tached to the die and measuring cells are deformed along with them. Their electrical resistance value changes. The 
strain-gauge station records this variation and, according with the results of the predetermined calibration of the sensors, 
converts it automatically into the force pattern. 

For the predetermined calibration of the strain gages, the die central segment and measuring cell were tested by 
compression on the MI-40KU hydraulic torsion-tensile machine. Usage of a torsion-stretching machine enables to de- 
velop considerable forces (35-40 kN). Herewith, their values are fixed accurately and reliably. The change in the on- 


loading stress-strain state causes deformation and a linear resistance variation of the strain gages. The principle of cali- 
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bration is in the construction of the circuit electrical voltage — applied force relationship. This dependence always has a 
strictly linear character. 

To perform the calibration, the elements under study were sequentially loaded in increments of 5 kN in the ze- 
ro — 35 kN range. The corresponding values of the circuit voltage on-load and after unloading were fixed. In order to 
increase the accuracy, three passes were performed over the specified range — up-down-up (42 measurements). As a 
result, for each element, the “force — displacement” relationship was obtained (“displacement” corresponds to the value 
of the upper striker travel). The calibration test data were statistically processed, and the linear regression equations that 
connect the force applied to the tool (P;, H) and the voltage (U;, mV) in the circuit were obtained from them. 


The equation for measuring the rolling force is as follows: 
P,;=—4852,7U; + 23538. 
The coefficient of determination R’ = 0.998. The standard measurement error based on the results of 42 tests 
was less than 0.2%. 
The equation for measuring the pressing force is as follows: 
P;= —4120,2U; + 21609. 
The determination factor R°= 0,997. The standard measurement error based on the results of 42 tests was less 
than 0.3%. 
The data obtained were entered on the program of registration and processing of measurements of the ZET- 
017-T8 strain-gauge station. This enabled to record the signal as a schedule of effort. 
After the calibration operations, the die with the strain gage was placed on the pilot plant (Fig. 4). 





Fig. 4. Force measurement under “rolling - ECA-pressing” combined process” 


Here, on the measuring cells of the first stand, strain gages were also installed to measure the rolling force 
which was performed only on the first stand. This is a key parameter according to the process concept. When imple- 
menting the combined process, the rolling force should exceed the resulting pressing force. 

The basic material used was 6063aluminum alloy (analog of AD31laluminum deformable alloy) at room tem- 
perature. Its chemistry and mechanical behavior are given in Table 1. 















































Table 1 
Chemistry and mechanical characteristics of 6063aluminum alloy 

Weight content of elements, % 

Si Fe Cu Mn Mg Cr Zn Ti Other 

0.2-0.6 0.35 0.1 0.1 0.45-0.9 0.1 0.1 | 0.1 0.15 
Mechanical properties at 20 °C 

Tensile stress, MPa Yield stress, MPa Hardness, HB Stretching, % 

90 48 25 20 
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Blanks of equal size (15 x 25 x 450 mm) were rolled with the thickness of 9 mm in the first stand. In the se- 
cond stand, the roller gap was 7 mm. The blank length was selected so that its lead end was bit by the second pair of 


rolls, whereas the back end was being still rolled in the first stand (Fig. 5). 
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Fig. 5. Force variation under deformation of 6063aluminum alloy: rolling in the first stand (a); pressing in the die (b) 


The force measurement data are presented in Table 2. 

















Table 2 
Force measurement results 
1“ stage 2" stage 3" stage 
Rolling step, kN 205 217 242 
Pressing step, kN 59 104 172 
Difference (friction margin), % 247 108 40 
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The difference between the values of the rolling and pressing forces was also determined. In this case, this dif- 
ference characterizes the level of active friction force margin: 


P 
A = (—2£_-1)-100%, 
TIPECC 


where Prpox is rolling force in the first stand; Pypgcc is back force in the die. 

A comparative analysis of the values of the arising forces has shown that the rolling forces at all three stages 
exceed the values of the corresponding pressing forces. Thus, the necessary condition for implementing the combined 
process is satisfied. Herewith, the active friction forces margin in the first stand with the advancing of the blank along 
the die mouths reduces dramatically. 

Conclusions. The study results indicate that the die strength margin is quite sufficient for the ECA-pressing 
implementation under the limiting conditions. It should also be mentioned that the rolling forces at all stages of defor- 
mation exceed the pressing forces, and thus, the necessary condition for implementing the combined process is satisfied. 
At this, the active friction forces margin in the first stand with the advancing of the blank along the die mouths reduces 
dramatically. 

The results obtained can be used for designing the experimental facilities that implement the combined process 


of interest. The stain-gauge analysis of strength characteristics is also suitable for the calibrated rolls. 
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Introduction. Generation of composite materials with addition 
of sawdust is considered. The wood component of the compo- 
sites is treated with water to enhance the surface roughness. 
This increases the contact area resulting in intensifying the 
sawdust — polymer fibers interaction. The work objective is to 
study the possibility of strengthening composites obtained from 
sawdust. 

Materials and Methods. Samples are made of composite mate- 
rials based on the unsaturated polyester resins, reinforced with 
wood chips. Water treatment was carried out at room tempera- 
ture for 2, 4, 6 and 8 days. Then the samples were tested for 
bending and compression. 

Research Results. As a result of testing the samples, changes in 
their mechanical properties were recorded. It is determined 
how the bending and compression resistance depends on the 
water treatment time. Graphs that reflect these dependences are 
constructed. 

Discussion and Conclusions. After water treatment, the com- 
posites reinforced with sawdust show a higher resistance to 
bending. This is due to the increased roughness of the sawdust 
surface and, as a consequence, to the extension of the surface 
area adhesion with the composite base. Besides, the water 
treatment enhances the specimen resistance under compression. 
The samples created on the basis of large sawdust come into 
particular prominence. This is due to the formation of holes on 
the sawdust surface which also enhances the adhesion between 


them and the composite polymer base. 


Keywords: wood fibers, sawdust, polymer, composite materi- 
al. 


“The research is done within the frame of the independent R&D. 


Beedenue. PaccmMaTpuBaloTca BOMpocb! co3qaHua KOMNO3HIH- 
OHHBIX MaTepHasioB C WoOaBeHHeM [{peBecHbIx onmMoK. J[pe- 
BeCHad COCTaBJIAIOWat KOMMO3HTOB OOpabaTEIBaeTCA BOTOH 
WIA YCHJICHHA WIepOXOBaTOCTH MOBepXHOCTH. OTO yBeM4HBa- 
eT TWIOW ab KOHTAaKTa, BCJIeCTBMe Yero YCHJIMBaeTCA B3aH- 
MoOjelcTBHe MexkKYy OMMWIKaMH HM BOJIOKHaMH MomMMepa. Lemp 
YccieqOBaHHaA — W3y4eHHe BO3MOXKHOCTH YIPOYHeHHA KOMIIO- 
3HIJMOHHBIX MaTepHasion, MOIYYCHHBIX 3 OMHJIOK. 

Mamepuaavi u memooei. OOpa3tbl U3TOTOBJICHEI 13 KOMMO3HT- 
HBIX MaTepHasIOB Ha OCHOBe HeHAaCbINeHHbIX TMOMM3:pHpHBIX 
CMOJI, aPMHPOBaHHBIX JpeBecHoi cTpyxKKoH. OOpadotKa Bool 
TIpOBOW1aCcb IIpH KOMHAaTHOM Temiepatype B TeyeHHe 2, 4, 6 u 8 
quel. 3aTem oOpa3iibI NpOXOTWIM HcnbITaHHa Ha W3rHO V1 CKa- 
THe. 

Pe3yiemamot ucciedoeanua. B pesyiibTaTe McibITaHHit OOpa3ijoB 
3ac:pPHKCHpOBaHbI H3MeCHeCHHA HX MeXaHH4ecKHXx CBOHCTB. Ompe- 
WeleHo, KaKHM 00pa30M CONpOTHBIeHHA H3rHOy HW MpouHocTu 
IIpH CoKaTHH 3aBHCAT OT BpeMeHH OOpaboTKH Boo. Hocrpo- 
eHbI Tpauku, OTparxkaroulve yka3aHHble 3aBMCHMOCTH. 
O6cyarcoenue u 3axmo4uenua. Llocne o6paboTKH BOOM KOMIIO- 
3HTbI, APMHPOBaHHbIe [PeBeCHbIMH OMWJIKaMH, J]e@MOHCTpH- 
pytoT OosIee BBICOKYIO CTOMKOCTb UpH u3rHOe. DTo OObACHA- 
eTCa YCHJICHHeM IepOxXOBaTOCTH MOBepXHOCTH OMMJIOK HW, Kak 
CileCTBHe, yBesIM4eHHeM MOBepXHOCTH CIeMIeHHA C KOMIIO- 
3MTHOM ocHOBOH. Kpome Toro, oOpadoTKa Bozo NoBbmllaer 
CTOHKOCTS OOpa3sl0B Mp cKaTHH. OcobeHHO 3aMeTHO ynpouHsA- 
TOTCA OOpa3l[bl, CO3{aHHbIe Ha OCHOBE KPYHHBIX OMHWIOK. ITO OOb- 
sCHAeTCA OOpa30BaHHeM yrslyOseHHii Ha TOBePXHOCTH peBeCHbIX 
OMMIOK, YTO TakoKe yCHJIMBaeT aflre3HIo MexKy HAMM HM TIOJIH- 
MepHOH OCHOBOH KOMIO3HTA€. 


KurroueBbie CJI0Ba: JIpeBeCHbIe BOJIOKHAa, ONWIKH, MOJIMMep, 
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Introduction. Composite materials based on polymers reinforced with various impurities are widely used in 
modern production practice. Major advantages of wood impurities (sawdust, chips) are their economic cost and availa- 
bility [1, 2]. 

To determine ratios of the sawdust wood fibers, a vibrating unit with a nest of sieves was used. In the end, 
three volumes of sawdust (5 mm*, 25 mm*, 120 mm*) were selected [2, 3, 4]. Through the experiments, the samples 
were subjected to three types of load (pressure, shear, and bending) to determine the effect of the sawdust volume on 
the strength characteristics of the composite materials. 

In 2001, samples with wood filler were subjected to the tests, the purpose of which was to determine the quali- 
ty of adhesion between fibers and binder. The results suggest that reinforcement with wood fibers increases the material 
[3]. 

In 2002, during the experiments, wood and wood fiber was chemically treated using MAPP (maleated poly- 
propylene). It was found that these samples have higher tensile and flexural strength [5]. 

In 2006, tests were conducted on the surface treatment of NaOH and CICH2COOH wood. The experiments 
have shown that such an effect enables to obtain high values of tensile strength [6]. 

Within the framework of this study, it is planned to demonstrate the practicability of the chemical treatment to 
improve the properties of the composite materials. 

Materials and Methods. Wood sawdust (5 mm*, 25 mm’, 120 mm’) [1] was treated by water at room temper- 
ature for 2, 4, 6 and 8 days. 

100 g of sawdust was placed in a plastic container covered with | liter of water and left for 2, 4, 6 and 8 days. 
Then the water was drained, and the mass was covered with distilled water for an hour. After that, the sawdust was re- 
moved and dried in the oven for 24 hours at 110 ° C. Further on, all tests were also conducted at room temperature. 

Research Results. 

Metering of weight loss of processed sawdust 





10 g of sawdust of different volumes (5 mm*, 25 mm*, 120 mm°) treated by the method described above were 
used. After processing, the samples were examined under an optical microscope. 

In this case, the weight loss of sawdust is equal to the weight difference before and after processing. The 
weight loss ratio is calculated as follows: 


Wy —W, 





=W%, 
1 


where W is the weight of the fiber sample before processing (g); W; is the weight of the fiber sample after treatment 
(g); W% is percentage of the weight loss (Fig. 1). 
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Fig. 1. Dependence of sawdust weight loss on time of treatment by water 
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The graph shows that 120 mm? sawdust does not lose much more weight than 25 mm? sawdust (the difference 
is about 4%). If we compare sawdust of 5 mm’* and of 120 mm’, the weight loss of the latter is 8 times higher. At this, in 
all samples, weight loss depends directly on the treatment time (the longer the test, the less weight). On exposure to 
water for long, some solubles go into the liquid, so, after washing and drying, the material weighs less. 

In the succeeding experiments, three composite samples were considered. For that, sawdust of various volumes 
(5 mm’, 25 mm’, 120 mm’) was impregnated with 10% polyurethane ether and mixed well. Then, a hardener was add- 
ed, and the mixture was poured into dies. Polyester was 90% of the total mass of the sample. 

Bending test of composites 

For the tests, samples of the composite material of 8 x 15 x 160 mm according to ASTM D 790 [7, 8] were 





prepared. 
Resistance to bending was measured at three points. Fig. 2 shows samples of the composite material reinforced 


with sawdust before the bending test. 





Fig. 2. Samples of composite material reinforced with sawdust before bending test 


Resistance of the specimen to bending (MPa) is determined by the formula: 
3-P-L 
Gis = 3 
2-Bet 
where P is load force (H), L is the sample length (mm), B is the sample width (mm), ¢ is the sample height (mm). 





Compression test of composites 
For the tests, samples of the composite material of @12.7 x 25.4 mm were prepared. Resistance was measured 





under the compressed condition. The load was applied until a rupture occurred. 
Samples of the composites reinforced with sawdust and prepared for testing are shown in Fig. 3. Resistance to 


compression (MPa) is determined by the formula [9]: 


i . C 2 
where P is compression force (4), F is c/s area (mm’). 





Fig. 3. Samples of composites reinforced with sawdust before compression test 
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Bending test results 


Fig. 4 shows that samples made of 5 mm? sawdust are more resistant to bending. 
g Ip g 
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Fig. 4. Change in resistance to bending against specimen treatment time 


Maximum bending resistance (230.88 MPa) was shown by the sample with 5 mm?’ sawdust after treatment for 
two days. The sample of 25 mm? exhibited maximum bending resistance (98.6 MPa) after the eight-day treatment. The 
sample of 120 mm’ exhibited maximum bending resistance (194.55 MPa) after the four-day treatment. 

The increase in bending resistance in this case is explained as follows. Under the treatment with water, the 
roughness of the sawdust surface increases, and the area of interaction with the base expands. Thus, the adhesion be- 
tween sawdust and polyester is strengthened and, as a consequence, the strength of the samples increases under bend- 
ing. 


Compression test results 





The samples from sawdust of 5 mm? are weaker in compression (compared to the samples of sawdust of 25 
mm? and 120 mm” (Fig. 5). 
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Fig. 5. Change in compression resistance depending on specimen treatment time 

The peak value of compression strength (156.67 MPa) is shown by the samples with 120 mm? sawdust treated 
within four days. This is due to the fact that under the water treatment, holes appear on the surface of larger wood saw- 
dust. This leads to the increased adhesion between the sawdust and the basic polymer, and enhances the compression 
resistance of the composite. 

Discussion and Conclusions. Testing of the composite samples reinforced with wood sawdust and treated 
with water allows for the following conclusions. 

1. Samples of 120 mm? sawdust lose more weight than samples of sawdust of 5 mm? and 25 mm’. 

2. Sample weight loss depends directly on the water treatment time. 


3. Under treating samples with water, the first sharp increase in bending resistance is observed for two days. 
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4, Maximum bending resistance (230.88 MPa) was shown by the sample with 5 mm? sawdust. 


5. Maximum compressive strength (156.67 MPa) was performed by the sample with 120 mm? sawdust. 
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Ultrasonic effect on electric spark forming and development in electroacoustic spraying * 


S. B. Kudryashev’, A. A. Zakalyuzhny”” 


Don State Technical University, Rostov-on-Don, Russian Federation 


Bumsanne y1bTpa3Byka Ha Wpolecchl POPMHpOBaHHA HU Pa3sBUTHA IIEKTPHU4eCKOH HCKPbI Ipu 
aK 


IIeKTPOaAKyCTHYeCKOM HallbLieHnn 


C. B. Kyapsames’, A. A. 3akamoxKHbtii? ™ 


'? Tlouckoii rocyzapcTBeHHbIii TeExHH4ecKHi yHuBepcute+r, r. PocroB-Ha-Jlouy, Poccuiicxaa Denepanua 


Introduction. The electroacoustic application of wear- 
resistance coatings is studied. The work objective is to obtain a 
mathematical model of the ultrasonic effect on the formation 
and development of an electric spark occurring in the process 
of the electroacoustic sputtering. 

Materials and Methods. The effect of ultrasonic vibrations on 
the processes occurring during the formation and development 
of a spark discharge is analyzed; the equations of continuity, 
energy motion and transfer, with the ultrasonic field 
contribution are considered. Factors affecting the thermal 
conductivity and electrical conductivity of strongly ionized gas 
are studied. 

Research Results. When obtaining the model, it was assumed 
that the heat removal from the channel is carried out by a 
“clear emitter”. Then, for the channel region, a self-similar 
solution is made: pressure, temperature and density are 
constant over the cross-section, and velocity is proportional to 
the radius. A mathematical model that describes the processes 
occurring in the spark channel with the ultrasonic field energy 
effect is obtained. 

Discussion and Conclusions. On the basis of the developed 
model, it is specified that under the ultrasonic radiation effect, 
the radius and temperature of the spark channel increase, and 
conditions of the double ionization under high ultrasonic 
energy are created. 


Keywords: electroacoustic spraying, hardening, hardening 
coatings, highly concentrated energy flows, ultrasound, 
electric spark, mathematical model, conductive channel, 
thermodynamics, ionization. 
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Beedenue. CtTaTba MOcBalleHa W3yYeHHIO mpollecca 
93IEKTPOaKYCTHYeCKOTO HaHeceHHa VM3HOCOCTOHKHX 
mokppirH. Llembro padoThl saBiaeTca nomyyeHne 


MaTeMaTH4eCKOi MOJeIM BIMAHUA yIbTpa3sByKa Ha MIpOrecchl 
(opmupopaHua uM  pa3sBHTHA  9seKTpH4eckKoii 
Tiporecce 


YCKpBI, 
mpoucxoysmmeli = iB 9JIEKTpOaKyCTH4YeECKOTO 
HaribuieHHA. 

Mamepuansi u memoovi. B ocHope axalM3a BIMAHHA 
YIIbTpa3BYKOBbIX KOJIeOaHHi Ha MIpoOuecchl, mpoTeKaromHe 
pH dopmMupopaHvuH UM pa3BHTHH MCKpoBoro pa3psaya, 
PaccMOTpeHbI ypaBHeCHHA HelIpepbIBHOCTH, JBY2KeHHA HU 
llepeHoca 9HeprHv C yYeTOM BKJlajja yIbTpa3BYKOBOTO TOIL. 
YutTeHbI (akTOpbl, BIMAOUIMe Ha TeMIOMpOBOAHOCTh U 
93JIEKTPONPOBOTHOCTh CHJIbHO HOHH30BaHHOLO ra3a. 
Pesyasmamei uccaedoeanua. Iipu nomysaenun Moyen OBuiM 
CieaHbl TpeqWOOKeHHA, 4YTO OTBO Tella M3 KaHasla 
OCYIJeCTBIAeTCA «pO3payHbIM U3yyaTesem». Torga ia 
oOsacTH KaHasIa ObWIO MPHHATO aBTOMOJelIbHOe pelleHue: 
aBlenue, TeMileparypa WW IWJIOTHOCTb TOCTOAHHbI 10 
CeYeHHHO, a CKOPOCTS NponopuMoHabya payuycy. onyaena 
MaTeMaTH4eckas ONMCbIBaIOlad —s HIPOIeCCBI, 
lmpoTekaroujve B HCKPOBOM KaHasie Cc y4eTOM BJIMAHHA 


MOJIEIIb, 


3HeprHu YIbTpa3sByKOBOPO TOJIA. 
O6cyarcdenue u 3akmo4uenus. Ha ocHoBaHuu WocTpoeHHoli 
MOJeIM yCTaHOBJIeHO, 4YTO NOX WelicTBueM yiIbTpa3Byka 
YBCJIM4YMBACTCA payMyc HW TeMIlepaTypa MCKpOBOrO KaHalla, a 
TakWKe CO3JaKOTCA yCJIOBHA JByKpaTHOH MOHH3allMH pH 
BbICOKHX 9HEPrHAax YIbTpasBykKa. 


Kosro4ueBbie 
ylpouHeHue, 


COBa: + 9JIeEKTpoakycTHyecKoe 
ynpouHsrouye 


BbICOKOKOHIUCHTPHPOBaHHble MOTOKH 9Heprun, yiIbTpa3ByK, 


HaliblIeCHHe, 
TIOKPbITHA, 


SIEKTpHAecKaAr MCKpa, MaTeMaTHuecKar MOJICJIb, 


TOKONpOBO TALIM KaHaJI, TCDMOJMHaMiKa, HOHHM3alnA. 
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Introduction. When applying wear-resistant coatings by electroacoustic spraying, a relatively narrow 
conductive channel with high temperature and ionization is formed, in which Joule heat is released, which leads to an 
increase in pressure and expansion of the channel. The expanding channel acts like a piston on the rest of the gas, and, 
since the expansion occurs at supersonic speed, it causes a shock wave, which propagates in front of this kind of piston. 
The temperature in the shock wave region is much higher than in the undisturbed gas, and the temperature in the 
channel itself is many times higher than in the shock wave. Accordingly, the density of the gas in the channel is very 
small; the vast majority of the mass of the moving gas is displaced from it, which makes it possible to consider the 
channel boundary as a piston [1]. 

The very fact of the formation of a narrow channel can be understood as follows: under the ultrasound action, 
as well as after the gas breakdown and the appearance of conductivity in the places of the current flow, Joule heat is 
released. The electrical conductivity of the gas is known to increase strongly with temperature. Thus, at a high degree of 
ionization, when the collision of electrons with ions is significant, the electrical conductivity is proportional to 7%”. 
With small ionization, this dependence is sharper, since with the growth of 7, the degree of ionization rapidly increases, 
and, consequently, there is a tendency to the current concentration in a relatively narrow channel. In places where the 
temperature is higher, the conductivity of the current is greater, so there is more current flowing, and more heat is 
released. This leads to even more heating and so on [2]. 

Physical processes determining the width of the channel and the limit of current concentration is heat transfer 
from the channel and the extension of the heated area under the action of pressure. The channel can be considered to be 
the area to the point where the temperature and the degree of ionization decreases significantly. In the channel, the 
inertia of the gas can be neglected, but it is necessary to take into account the release and transfer of heat. In the shock 
wave region, inertia must be taken into account, but electrical and thermal conductivity can be neglected. These two 
areas are separated by a transition layer — the “shell” of the channel. Heating and ionization of the gas entering the 
channel occurs in the shell [3, 4]. 

Materials and Methods. The main equations are the equations of continuity, motion and energy transfer 
taking into account the action of the ultrasonic field. They have the form of [5, 6] 

dp_y OP, a(rV) _ (1) 
ot or ror 





(2) 


OV _OV)\ Op 
+V + =0 
(e _ or (3) 


2, 2 fa) 
per Pe be rpV pa! + (0) _ jesu 
ot 2 ror p 2 ror 


where ( is density, V is velocity, p is pressure, € is internal energy per unit mass of gas, q is heat flux, / is 











current density, E is electric field, U is contribution to the ultrasonic field. 
The equation of state has the form 


p=(n,+n,)T =pT(z+1)/m, (4) 
where m, is the average mass of an atom; 7,,n; is the number of electrons and ions per unit volume; z is the 


average charge of an ion; n.=zn,. Temperature is expressed in energy units. 
We assume that the ionization in the channel can be calculated by Saha ionization equation. The internal 


energy of the gas is 

3p I pls. 7 (5) 

e= i +—= + 
2p m, pl2 (zt+I)T 


where / is the total ionization energy plus the dissociation energy per atom. The formula (5) is convenient to 





use in the case of full ionization; with incomplete ionization with increasing 7, the value of / also increases. At the same 
time, as follows from Saha ionization equation, //T~const, so for €, the following formula is more convenient: 
1 
e=— iz 
y-1P 
(6) 
where /is the effective adiabatic index; for air y=1.22. 


Electrical conductivity O and thermal conductivity v of highly ionized gas are equal to 
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o=0,(T)T*? =30'(z)T%” /(4e? l2nm.) 
(7) 
aie” 
Here, e, m is the charge and mass of an electron, a= In(37 ue / (ze? 4nn,)), Oz) is the dimensionless 


coefficient. For z=1, the value is 0/= 1.95. The value y e*/o7 , under the Wiedemann—Franz law, of the Ist order. For 


A=5 O1(1)=3,4.1013sec—leV-3/2, 

X1(1)=3,9.1020 cm-1 sec-leV-5/2. 

Research Results. Let us suppose that the time dependence of the channel radius, the boundary of which plays 
the role of the piston displacing the gas, has the form of a(t)=Atk; the motion in the shock wave region is determined by 
two dimensional parameters A, /. 

In the equations (4)-(6), we introduce dimensionless notation (a, — the radius of the wave front) 





x=r/a.(t), p'(x)=p/Po> V'(x)=V/4.. P(x) = P/ Poa (8) 
Neglecting heat release and transfer, let us write the system (1)—(3) in the form 
; dp’ _,dxV' 
V —_ ————i —<=<$<—=. 

( *) dx ms xdx 
i 9 
=z)" aad a er 0) 

k dx ' ax 





1 dp' dxV' 
2) 1-—|p'+(V' +y p’ =0 
( i} OSA) Pe 


with boundary conditions at x=1 
pe’ =(v+1)/(y-1), V' =2/(y-1), p’=2(y-1) (10) 
The position of the piston is determined by the point where V =x. The pressure on p, piston can be expressed 


through the piston speed. 
Px =K,,Poa°, (11) 


where “resistance coefficient” K,*0,9 (K,=p (a)a2/a2,) is deduced from the numerical solution of the system (9). 


We neglect the radiation and assume that 


at (12) 
pe dr 
Let the temperature 7 in the channel be much greater than that required for the complete ionization, therefore, 


on the edge, 7 is much less than in the center. Let us suppose that 7=0 at r=a. We introduce dimensionless notations 


f° T lr(r V r}q 5S(V +r (13) 
= 2 2 0(s)= 2 u= Y= + , 
a’(t) is 0@2a\la a 2al|pa 2\a a 


where 7) is the temperature on the axis. We consider the pressure to be constant along the channel cross-section. 

If the heat is removed from the channel by a “clear emitter”, a simple self-similar solution can be specified for 
the channel area: pressure, temperature and density are constant over the cross section, and the speed is proportional to 
the radius. The temperature drop is concentrated in the shell. In the same place, the radiation is absorbed and the 
ionization of the gas entering the channel takes place. Considering the shell to be thin, it is possible to obtain a system 
of equations for the main parameters of the channel. In general, you can use these equations to estimate them as a 
mathematical model that describes, although roughly, the main processes in the channel. This takes into account the 
approximate action of ultrasound and thermal conductivity [7]. 

The energy balance equations for the channel and the shell have the form 











dW dna (14) 
—= =0.+0,, 
a a 2; +2 
(15 
[er2 | 949, ) 
p) dt 


where M, W is the mass and energy of the gas in the channel. The equation (15) is obtained by integrating (3) 
over the cross section of the channel (including the shell) without assuming the form of distribution of values over the 
cross section. For a homogeneous model, let us suppose 


Kudryashey S. B., Zakalyuzhny, A.A. Ultrasonic effect on electric spark forming and development in electroacoustic spraying 





W=M-s,M=n-a’p. 
The equations (14), (15) are the consequence of the law of energy conservation. The expressions for the release 
of heat Q; due to the electric field and heat Oy due to the ultrasonic field, as well as for the heat removal by radiation Op 
and thermal conductivity QO; , can be taken as 


O,=;'/na’s, O,=no (16) 
0,=na¢O(pT), O,=1,3-20y T (17) 


where @is the frequency of ultrasonic vibrations, 77 is the dimensional coefficient. 
Comparing (14) and (15), we obtain that 


O, +O, =» (Q,+Q,) . (18) 
where J is a coefficient of the Ist order. If T is independent of ¢, then 
aa aul i (19) 
=y|1+(y-1)2a’ 
Bayt) - ) 


Let us consider the channel in air with the conductivity of oO =2'10'4sec"; Kp=0.9; v=1.2; j~t; hence & =4,5. 
For the radius of the channel, we get the expression 





a =0,93(1+0) p76 i412 eo) 


Here, a is measured in mm, j — in kA, ¢— in pts; Pp =1.29'10° g / cm’ under the atmospheric pressure. If the 
shock wave is weak, the radius is similarly deduced from (19). 
Table 1 gives values of the radius calculated by the formula (20), with different values of @ and ¢ (1s) at the 


discharge voltage V=30 V, and the battery capacity c=0.15 pF, and inductance of circuit L=4 nH (corresponding to 
j=V/L=7.5 10° Alsec). 




















Table 1 
The channel radius for different values of 6 t 

T 

6 0.3 0.5 1 

0 0.65 1.00 1.62 

1 0.73 1.12 1.82 

2 0.78 1.20 1.95 

3 0.82 1.26 2.04 














Let us estimate the temperature in the channel. We believe that [t ~1, for the same discharge as above at the 
time f=1 ps at L=4 nH we have that O,+Q, = (1+0)1,7-10" erg / cm sec. If we assume that all heat is transferred by 
the electronic thermal conductivity, and the radiation is neglected, we obtain that 
T4(1+O0)eV. Taking T=4eV we find that the number of ions per unit volume in this case is n; =9.1017, which in order 
corresponds to the experimental values. 

Discussion and Conclusions. Based on the constructed approximate model, the following conclusions can be 
made about the effect of ultrasound on the development of the spark channel. 

1. The radius of the channel increases by (1+0)1/6 times in comparison with the case when there is no 


ultrasound, where 0 is the ratio of the energies of the electric and ultrasonic fields. 


2. The temperature in the channel increases proportionally (1+9) under the assumption that the outflow of heat 
is carried out by the electronic thermal conductivity. 

3. Already at the moment of the shock wave formation, almost complete ionization occurs in the channel, and 
conditions for double ionization at high ultrasound energies can be created. 
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'? Tlouckoii rocyapcrBeHHblii TeXHH4ecKHii yHuBepcuter, r. Pocros-Ha-Jlony, Poccniickaa Denepannsa 


Introduction. The sources of damage and wear of the main shaft 
of the drier toaster are analyzed. The repair know-how and 
welding operations execution limitations which must be consid- 
ered when developing the technique providing the restoration of 
the structure performance features are studied. The work objec- 
tive is to develop a technique of repair without dismantling for 
the main toaster shaft. To solve the task, a design repair struc- 
ture was installed, and postwelding operations that meet the 
engineering and regulatory requirements for this structure were 
performed. 

Materials and Methods. In “Kompas 3D” software, the follow- 
ing models were developed: integral shaft (project shaft design); 
damaged shaft as a result of long-term operation (more than 15 
years); and damaged shaft with a welded repair structure. Nu- 
merical simulation of the stress-strain state (SSS) was carried 
out. 

Research Results. Software for the computational modeling of 
the repair structure SSS is developed. The repair shaft structure 
in which the maximum stresses do not exceed the shaft stresses 
in the project design is obtained using the model. To eliminate 
the aggressive medium effect on the corrosion fatigue strength 
of the shaft, an insulating method is used. A technique for 
mounting the repair structure to the shaft allowing for the out- 
run limitation 0.12 mm is developed. 

Discussion and Conclusions. Torsion shafts damaged deeply by 
wear and corrosion are considered. To restore their structural 
integrity, it is worthwhile using the following complex of tech- 
niques: 

— constructive (consists in the installation of optional parts that 
compensate for insufficient strength, and provides a reduction in 
stress concentration in the most loaded zones); 

— processing (reduces residual welding stresses due to the rea- 
sonable sequence of deformation that contributes to generating 
favorable residual compressive stresses); 


— isolation (is based on the application of anticorrosion coat- 


ings). 


* The research is done within the frame of Contract No. 130 of 19.04.2014. 


* E-mail: lyudmirskiy40@mail.ru, assaulenko_s@mail.ru 
“ PaGota BbIOmHeHa 10 ZoroBopy Ne 130 or 19.04.2014 r. 


Beedenue. Crarba MocBaljeHa avalM3y IpH4ynh NoBpexqeHuit 
WM wW3HOCa TylaBHOrO Basia TOCTepa cylmMiKH. PaccMoTpeHbl 
cieljHuka PEMOHTHOTO MpOH3BOACTBa H OrTpaHH4eHHA B BBI- 
TIOJHCHHM CBapOuHBIX paOoT, KOTOpbie HeOOXOZHMO y4HTBI- 
BaTb IIpH CO3qaHHH TeXHOJOrMH, OOecreyMBaloleli BOCCTa- 
HOBJICHHE SKCIUIyaTallMOHHBIX CBOMCTB KOHCTpyKUMH. 

Lemp padotpl — co3qaHHve TexHOJIOrHH PeMOHTa TaBHOTO 
BaJla TOCTepa 6e3 ero pa300pxu. Jia pearu3alun 9TOM 3aqauH 
Oblia ycTaHOBsIeHa pacueTHad PeCMOHTHat KOHCTpyKUMA U 
BBINOJIHEHBI MOCHecBapoyHble onepalwu, oOecnewnBaromjue 
TpeOoBaHHa TexHHYeCKOM H HOPMaTHBHOM JOKYMeHTalHH K 
WaHHOW KOHCTpyKUHH. 

Mamepuaaet u memodvi. B mporpamMMHoli cpeye «Kommac 
3D» pa3paOoTaHbI Moyes: Wemoro Basia (KOHCTpyKUHA Balla B 
IIPOCKTHOM COCTOAHHH); Balla, MOBpexKCHHOTO B pesyIbTaTe 
TIMTeIbHOM 9KcIIyaTauMH (Oomee 15 mer); MoBpexqeHHOTO 
BaJla C IIpHBapeHHOH peMOHTHOM KOHCTpyKuHeli. [IpopeqeHo 
YHMCIeHHOe MOeMpOBaHHe HalipsKeHHO-ePpOpMUpOBaHHOTO 
coctosHua (HC) mMogenei. 
Pezyiemamet ucciedoeanua. Pa3spadoTaHo mporpaMMHoe 
oOecneyeHue AJId YHceHHOrO MozemMpoBanHua HJIC pemont- 
HOM KOHCTpykKHH. C MOMOM[bIO MOJeIM MOyueHa KOHCTpyK- 
WMA PCMOHTHOrO Balla, B KOTOPOM MaKCHMaJIbHbIe HalipsKe- 
HHA He TIpeBbUalOT HallpsKeHHH Basia B IIPOCKTHOM COCTOS- 
Huu. Jd HCKIOUCHHA BIIMAHHA AarpeCCHBHOM CpeybI Ha KOp- 
PO3HOHHY!0 YCTAJIOCTHYIO IPO4HOCTh Basia HCHOJIb30BaH H30- 
JIAWHOHHEIM cocoé. Jia mpHcoeqHHeHHA PeMOHTHOHM KOH- 
CTpyKUMH K Bally pa3paOoTaHa TeXHOJIOrHA BBbINOMHeHHA 
UIBOB, MO3BOJAIOMat OrTpanvu4ntTb Onenue 0,12 mM. 
O6cyarcoenue u 3akmouenua. PaccMoTpeHs! padotaromve mpu 
KPY4eHHH BaJIbl, MOBPeCxKCHHbIe B pe3yJIbTAaTe H3HOCAa MU KOp- 
po3uu Ha Oombutyio rryOuny. Jt BOCCTaHOBICHHA UX KOH- 
CTPYKTHBHOM MpouHOcTH WerecooOpa3HO HUCIONb30BaTb KOM- 
TWIeKC TepeyHCJIeCHHbIX Jalee MeToqoB. KoHcTpyKTHBHbIi 
(gakKmOuaeTCA B YCTaHOBKe J{OMOMHHTEIbHBIX JeTamel, KOM- 
TICHCHPYIOWMX HeOCTAaTOYHYy!IO MpOyHOcT, HW ObecneuHMBaeT 
YMeHBbUIeHHe KOHICHTpallui Hallps2KeHHH B HanOoslee Harpy- 
2KCHHBIX 30HaX). TexHOJIOrHyecKHM (yMeHbUIaeT OCTATOUHBIe 
CBapO¥Hble Hallps2xKeHHA 3a CUeT PallMOHaJIbBHOM MOcIeqOBa- 
TeIbHOCTH JesopMupoBaHHsA, clocoOcTByIOMIero HaBeeHHIO 
OnaronpHATHBIX OCTATOUHBIX HallpsKeHHH cxKaTua). V30nn- 
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The economic expediency of the developed repair technique is pyloljui (OCHOBaH Ha HaHeCeHHM aHTHKOppO3HOHHBIX TI0- 


obvious. The repairing of the shaft without dismantling costs kppituit). 


180,000 rubles, while a new shaft costs 3.8 million rubles. SKOHOMMYECKAA HeMeCooOpasHOcTh paspaboraHHOH TexHO.0- 


THH PeMOHTa HecoMHeHHa. PeMoHT Bala 6e3 ero pa300pKH 


ctoun 180 Trica pyOnel, a WeHa HOBoro Bala — 3,8 MIH 
pyOumei. 
Keywords: welding repair under factory conditions, shaft re- Kusouesbie C10Ba: PEMOHT CBapKO B MpOH3BOJ{CTBCHHBIX 
ry P Pp. Pp’ 
pair, repair type selection, 3D modeling, finite-element method YCHOBHAX, PeMOHT Basia, BbIOOp Tuma pemouta, 3D- 


(FEM), optional parts, stress factor, repair technique, plastic MOR RBD O EaH BGs CTO hos sae ee ce eae 


TIOJIHHTCJIBHBIC WeTasiM, KOHWCHTpauyHa Halips»KeHHH, TeCXHO- 


deformation, economic expediency. 
JIOTHA PeMOHTa, IylacTHyeckoe edpopMupoBaHHe, IKOHOMH- 


uecKkaa leslecooOpa3Hocts. 
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2018-18-3-311-317 — T. 18, Ne 3. — C. 311-317. https://doi.org/10.23947/1992- 
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Introduction. Failure of machines is often due to the negative impact of structural, technological or operation- 
al factors. To choose the type of repair, you need to know: 

— operational aspect features , 

— conditions of its operation; 

— possible failure causes of the main elements. 

The specificity of the repair is that you have to deal with already created structures during work performance. 
This requires the selection of certain welding methods, the location of repair welds, restricts access to welding sites. In 
some cases, the possibility of heating the metal during welding and subsequent heat treatment of the welded joint is 
reduced [1]. All this leads to the need for developing special technical solutions in each case. 

Oil extraction plants for the distillation of solvent (gasoline) from the meal (cake) use an evaporator (toaster), 
equipped with heating jackets and steam supply pipelines. The toaster, shown in Fig. 1, is a column apparatus assem- 
bled from vats (1) mounted one above the other. 
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Fig 1. Scheme of dryer toaster 
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Through all the vats, shaft 2 passes, on which pair agitators of meal 3, installed above the bottoms of the vats 
with a small gap, are fixed. The shaft drive is provided through the reduction gear 4 of the engine capacity of 160 kW 
and rotates with a frequency of 10 Rev/min. Torque on the shaft of the toaster is 14500 kgem. 

Dryer toaster, investigated in the framework of this work, was in operation at the oil extraction plant of the 
“MEZ Yug Rusi” branch for 15 years, which led to significant local wear of the shaft. 

The design dimensions of the shaft are as follows: diameter is 285 mm, length is 6860 mm. As the diagnostics 
had shown, the shaft wear was localized in the lower part, where it was in contact with the meal (acting as an abrasive). 
The temperature of the meal was 105 °C. The contact point was blown with moist air. In this part of the shaft at the site 
of 700 mm, its diameter decreased to 245 mm as a result of corrosion and mechanical wear. 

The worn part of the shaft is shown in Fig. 2. 





Fig. 2. Shaft wear after 15 years of operation 


The calculations performed according to [2] showed that the polar moment of inertia of the shaft in the design 
was 6.6*10* cm’, and due to the operation, it decreased to 3.8*10* cm*. Further operation of the shaft without restoring 
its bearing capacity could lead to the destruction and explosion of gasoline vapor in the toaster. In this regard, various 
options for repair of the main shaft of the toaster were evaluated. Replacement was not considered due to the high cost 
of the individual shaft production. In addition, its dismantling required disassembly of the toaster and the room where it 
was installed, which was almost impossible. 

Currently, the most widely used methods of repairing worn shafts are the following: 

— welding metal over the damaged area [3] to restore the size and the carrying capacity of parts; 

— replacement of the damaged section of the structure with a new one; 

— installation of additional parts or structure repair [4]. 

The first repair option was rejected, because in this case it was required to weld a large amount of metal (more 
than 140 kg) on a vertically located shaft, which would inevitably lead to defects and large deformations. 

The second method is difficult. In its implementation, there would be difficulties with cutting the damaged part 
off the shaft, cutting edges for welding, assembly and welding. This would result in unacceptable welding deformations. 
In this case, even with a successful replacement of the damaged area, the operational properties of the toaster would 
most probably be violated. 


The third method is usually recommended for repair if the surface damage is large [4, 5]. Thus, it was decided 
to repair the main shaft of the toaster by installing a repair structure (Fig. 3). 
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Fig 3. Scheme of repair structure 


The design of the shaft 1 consists of two half-rings 2 and two half-couplings 3. The sizes of the coupling and 
the ring were determined according to strength and technological considerations. The length of the couplings equal to 
800 mm, exceeded the length of the damaged section of the shaft by 100 mm. The inner couplings diameter is 285 mm, 
and the outer one is 325 mm. To restore the static strength of the shaft in its worn parts, it was needed to weld the ele- 
ment, the polar moment of inertia of which would make up for its wear loss (6.6°10* cm* — 3.8*10* cm* = 2.8¢10* cm’). 
Then, the above-mentioned diameters of the polar moment of inertia of a coupling half is equal to 3.1*10* cm*. Thus, 
the dimensions of the repair structure welded to the shaft not only restore the calculated moment of inertia of the shaft, 
but also increase it by 15 %, which provides static strength of the shaft. 

The papers [6, 7, 8] show that the durability of the welded designs working at repeatedly static loadings in the 
damp environment, mainly depends on the stress-strain state (SSS) where weld metal comes to the basic metal. The 
finite element method (FEM) was used to estimate the effect of weld geometry on the SSS. When designing the repair 
structure, we sought to ensure that the stresses in the weld area did not exceed the stresses in the area of keyways on the 
shaft to transfer torque to the meal agitators. 

The method of SSS calculation is presented in [9]. In the “Compass 3D” software environment, 3D-models of 
the complete shaft with the original keyway, damaged by the shaft corrosion, and the damaged shaft, were developed 
with account of the installation of the repair structure. The Ansys software product [10, 11] was used for the calculation. 
The 3D model was imported in Ansys Workbench. With the help of this software shell, pre-processing of geometry and 
selection of the material corresponding to the material of the shaft-17G1S was carried out (Table. 1). 














Table 1 
Parameters of the selected material 
Property Unit of measurement Value 
Density kg/m? 7850 
Yield stress mPa 250 
Ultimate tensile strength mPa 460 

















Further, the calculation was carried out, which includes the creation of a finite element mesh, as well as the de- 
termination of boundary conditions by applying torque from the reducer and rigid fastening of the shaft on the reverse 
side. 

Fig. 4 provides the results of the calculated structures, and Table 2 - maximum stresses on the shaft: in the are- 
as of keyways (Fig. 4, a), in the place of shaft wear (Fig. 4, b), in the most loaded area of the repair design (Fig. 4, c). 
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a) b) c) 


Fig. 4. Models of stress-strain state calculation in areas of: keyways (a); shaft corrosion (6); repair structure(c) 














Table 2 
Maximum calculated stresses 
Area of stress calculation on the shaft Maximum voltage, mPa 
Keyways 74 
Shaft corrosion 124 
Repair structure 75 














It can be seen that the stresses in the developed repair structure do not practically exceed the design stresses in 
the keyways areas. 

SSS depends not only on the operating voltage, but also on the magnitude and the sign of residual stresses in 
the area of weld metal connection with the main metal. Therefore, the method of plastic deformation was used in the 
repair technology [12]. It allows you to reduce residual tensile stresses or even change them into favorable compressive 
stresses. 

Before starting, degassing of the toaster and the shop where it was located was performed. The damaged part 
of the shaft and the welding site were cleaned with rotating brushes to the metallic sheen. 

The assembly of the repair structure started with the installation of two half-rings 2 (see Fig. 3), their tighten- 
ing up and tack welding. To weld the half-rings, the UONI electrodes 13/55 [13] were used in the modes specified in 
Table 3. 




















Table 3 
Modes of manual arc welding with coated electrodes 
Electrode brand Recommended current strength when AWM * A, in the position: 
Diameter, mm 
bottom vertical horizontal and ceiling 

UONT 13/45, 3.0 80-100 60-80 70-90 
PONE D> 4.0 130-160 100-130 120-140 

5.0 170-200 140-160 150-170 

















* AWM-arc welding mode. 











After that, the welds were cleaned flush with the base metal. 


Machine building and machine science 


315 


http://vestnik.donstu.ru 


pa 


Vestnik of Don State Technical University. 2018. Vol. 18, no. 3, pp. 311-317. ISSN 1992-5980 eISSN 1992-6006 





Then the ring was set perpendicular to the shaft, as shown in Fig. 3, and welded; circular welds were welded 
connecting the ring 2 to the shaft, according to the scheme shown in Fig. 5, a. 


been | mt 





. 


' \ : 


























a) 
b) 


Fig. 5 Diagram of implementation of: circular welds (a), longitudinal welds (b) 


The welding was performed simultaneously by two welders. The order of bead welding (1-18) is shown in Fig. 3. 

Two coupling halves were assembled on the welded ring, so that the distance from the longitudinal welds of 
the ring to the longitudinal welds of the coupling was 120-130 mm. Half couplings were tightened around the shaft by 
the means of two centralizers; at each junction there were three tack welds with the length of 50-60 mm. To reduce the 
probability of occurrence of defects during welding, the beginning and the end of each tack weld was rubbed down with 
an abrasive tool. After this, the upper centralizer was removed, and the welding of longitudinal welds of the coupling 
started. The welding was performed by the step-back method [14] from the bottom up — in the order shown in Fig. 5, 
b. After performing 70 % of the length of the longitudinal weld, the lower centralizer was removed, and the root weld was 
completed. Filler and cap beads were done in the same order. Each subsequent layer was shifted relative to the previous one 
by 25-30 mm. The current strength recommended in Table 2 for vertical joints welding, was reduced by 10 %. 

The circular welds connecting the coupling to the shaft and the ring were made simultaneously by two welders 
in diametrically opposite places relative to the shaft (see Fig. 5, a). First, the beads 19 to 27 and then, 28 to 37 (see Fig. 
3) were welded. After circular welds were welded, their deseaming and plastic deformation at the edges of the weld 
metal with the base metal of the shaft was made. The quality of the welds was checked by the visual inspection in lay- 
ers. 

The radial-motion variation of the shaft was measured by hour-type indicators above the repair structure by 
100 mm in two mutually perpendicular directions. Its value did not exceed 0.12 mm. 

In order to protect the coupling and the shaft from corrosion, the space between them was filled with sunflower 
oil [15], and the outer surface of the shaft was insulated with epoxy glue of the EDP brand according to TS (technical 
specification) 07510508.90-94. 

The integrity of the repaired shaft is monitored quarterly. It has been working reliably for four years with full 
compliance of functional indicators with the requirements of technical and regulatory documentation. 

The economic feasibility of the developed technology is significant. The repair cost 180 thousand rubles, 
whereas the price of a new shaft is 3.8 million rubles. 

Conclusions. To restore the structural strength of the shafts working in torsion, it is advisable to use a set of 
the following methods. The constructive method involves the installation of additional parts that compensates for the 
lack of strength, and provides a decrease in stress concentration in the most loaded areas. The technological method 
reduces residual welding stresses due to the rational sequence of deformation, which helps to induce favorable residual 
compression stresses. The insulating method is the application of anti-corrosion coatings. 
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Bred 
Harpy3KH Ha padouem oprane MOOMJIBHBIX Mall 


B. A. Tepmmu’, T. A. XunuKaz3e” 


'? Wucrutyt ciepsi oGcryxuBaHna M MpesnpHHumarenbersa (unuan) r. Uaxtst oucKoro rocyapcTBeHHOro TeXHM4eCKOrO yHUBepcuTeta, 


Poccuiicxaa Dexepayua 


Introduction. Issues on the functional unification of the adaptive 
hydraulic drive module are studied. For the first time, self- 
adapting mechanisms are considered taking into account 
adaptive intercommunication of the load control and agreement 
of motions on the working body of the mobile machines. The 
work objective is to create and analyze the technique of the 
functional unification of the adaptive hydraulic drive module. In 
the furtherance of this goal, a number of tasks are solved. The 
selection of technical equipment — unified adaptive hydraulic 
drive modules of the mobile machines — is validated. The 
methodology and indicators of the module functional 
unification are described. Intercommunications are considered: 
direct positive and back negative ones. Their effect on the 
functional unification property of the adaptive module is shown. 
Materials and Methods. For the synthesis and analysis of the 
functional unification indicators of the adaptive module, a 
similarity method of the technical systems operation is adopted. 
structural-functional 


Research Results. Techniques for 


unification of the self-adapting modules are developed. 
Optional versions of the unified modules modification and 
proper combinations of hydraulic motors, regulating equipment, 
and mathematical models of adaptive communications are 
presented. Criteria and indicators of similarity are proposed. 
The functional unification of the adaptive intercommunications 
of the module and different types of the hydraulic motors and 
fluid throttling elements in the hydraulic system are analyzed. 
Recommendations for implementing the functional unification 
under typing and operation (adjustment) of the adaptive module 
are formulated. 


Discussion and Conclusions. The methodology is recommended 


for the functional unification of the hydraulic self-adapting 


* The research is done within the frame of independent R&D’.Padora BpmomHeHa B paMKax MHuUMaTHBHO! HP. 


* E-mail: vpershin2013@gmail.com, khinikadze@mail.ru 
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Beedenue. VUccneyqopanbt Bompochl dyHKUMOHaIbHOn 


yHuukalu ajanTHBHoro MOsysIA THyzpompuBoga. Bnepsere 
paccMOTpeHbI MpHHIMMbl caMoayalraywH Cc y4eroM 
BHYTPCHHHX ajallITHBHBIX CBA3eH yIpaBeHHA Harpy3Koli u 


COrylacOBaHHOCTH JBMKeHHi Ha padouem opraHe MOOWJIBHBIX 


Mallu. 
Ilemp = padoTbr — co3qaHHe uM aHav3 MeTORHKM 
(YHKUHOHAIbHOH = yHuuKalyMu = aalTHBHOTO = MOys14 


THAPOMpuBoya. J{na TOCTWKCHMA MOCTaBIICHHOM eI penieH 


pia 3aqau. OOocHopaH BEIOOp TexXHHYeCKHX yCTpolicTB — 


YHHUUMpyeMBIX ajalITHBHBIX MOyyNel ruyqpoupuBora 
MOOWJIBHEIX Maui. II[pequoxKeHbI MeTOqMKa HM WoKa3aTesI 
(YHKUHOHAIbHOH = yHuduKayHH Mogyia. Paccmotpexsi 


BHYTPpCHHHe CBA3H: TWpxAMavA TOJIOKUTCIbHAA U oOpatHaa 


oTpuuatembyHas. JloKa3aHo ux BIIMAHHe Ha CBOLMCTBO 
(PYHKUHOHAIBHOM yHupukalHu alalnTHBHOrO MOLysIA. 
Mamepuaavi u memoovi. Jina cuuTe3a HW aHas3a MoKa3aTeset 
(PYHKUHOHAIbHOM yHUpUKalHN aanTHBHOTO MOY IpHHAT 
MeTO] Toq0OuA (PYHKUMOHMpOBaHHA TeEXHHYeCKHX CHCTeM. 
Pe3zyiemamet —succaedoganus.  PazpadoTaHa  MeToyMKa 
TIpOBeeHHA KOHCTPyKTHBHO-(yHKUMOHaIbHOH yHupuKalHn 
Monynei c camoayanrayueit. IloKa3aHbI BO3MO2KHBIe 
BapHaHTbl MOAMPuKaWMH YyHHPMYMpyeMbIxX Moyle u 
COOTBETCTByIOWIMe HM coyeTaHHa ruyqpoyBurareseit, 
perylupyroljei alimaparypbI Hw MaTeMaTHYecKHxX MoOjerIel 
alalITHBHBIX cBA3el. [IpeqoxeHbl KPHTepHH HW MHMKATOPEl 
mogoous. BrmowHeH avarM3 (PyHKUMOHAaIbHOM yHuduKalyHu 
BHYTPeCHHHX ajjallTHBHBIX CBA3eH MOJYIA C pa3HbIMH THNaMH 
THAPOABHTaTeseH HU dIeEMeHTOB ApPOccesMpOBaHHA %KMAKOCTH 
B  rHqpocucreme. CidopmymmpopaHbI pekomeHyalHu 
TIpoBeqeHuaA PyHKUMOHAIbHON yHuuKalMN pu TuMM3alHH 
Wf OKCIIyaTalun (HaslayKe) ayjanTHBHOrO MOAysIA. 
O6cyarcdenue u 3axiiouenun. Metoquka peKoMeHyoBaHa WIA 


(PYHKUHOHAILHOM yHuPUKayHH Tu~paBIMYecKorO MOA Cc 
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module. It can be used for the development of unit sizes and caMoaylantauvei. OHa MOXKeT  MHCMONb30BaTBCA IPH 


under its operation as an independent drive or a hydraulic drive paspaOoTKe TuMopa3mMepoB Moya HM B Mpolecce ero 
. . . . KCIIyaTalwH B KayecTBe aBTOHOMHOTO MIpHBoya HIM 

subsystem of a multifunctional or combined machine. oer Bua eS SOM BOL SEU EBOL 

HOJCHCTeMbI THApONpHBoAa MHOro@yHKUMOHaIbHOH WIM 
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Introduction. Multifunctional and combined mobile machines are abundantly used under difficult conditions: 
in road construction and agricultural work. As a rule, their operation involves manual or hard-programmed control. 
Thus, it is important to create simple functionally unified drives (modules) of mobile machines capable of adapting to 
the variable properties of the processed environment [1]. 

It should be noted that the issues of unification are related to the solution of a number of scientific and practical 
problems. One of them is the problem of standard sizes of functionally unified products. For its solution, the 
corresponding means of fundamental [2] and applied [3, 4, 5] mathematics is used. Studies on the adaptive systems of 
different structures and purposes [6], hydraulic drives of devices with adaptively coordinated movements of the 
working body taking up stochastic fluctuating loads are carried out [7, 8, 9, 10]. In the framework of this paper, the 
principle of functional unification in the modular hydraulic drive capable of load self-adaptation is considered. The 
research and practical testing of such a module is conducted for the first time. 

Materials and Methods. A device for drilling rocks with variable properties [11] and a device for processing 
curved surfaces [12] are accepted as the basic objects of research a 1, 2). 


8 





Fig. 1. Schematic diagram of device for drilling rocks with variable properties [3] 


Machine building and machine science 


319 


http://vestnik.donstu.ru 


1eS) 
N 


0 


Vestnik of Don State Technical University. 2018. Vol. 18, no. 3, pp. 318-325. ISSN 1992-5980 eISSN 1992-6006 





Fig. 2. Schematic diagram of device for processing curved surfaces [2] 


The following elements are numbered in the diagrams (for Fig. 2, designations are given in brackets): 1 (10) — 
fixed displacement hydraulic pump; 2 (11) — safety valve; 3 (12) — filter; 4, 7, 8 — control valves; 5 (4) — choke 
(mode control); 6 (1) — main motion hydraulic motor; 9 — flow control; 10 — retaining adjustable choke; 11 — feed 
hydraulic cylinder; 12, 13, 14 — pressure gauges; 15 — tank; (2) — working body; (3) — hydraulic system; (5) — 
reversible feed motor; (6) — feed motor shaft; (7) — feed mechanism; (8, 9) — adjustable chokes. 

The elements of the hydraulic system of these devices (including their modifications [13]) are structurally 
combined according to the differential scheme. 

A hydraulic cylinder (see Fig. 1) or hydraulic motor (see Fig. 2) are additionally included between the return 
lines of the mode control and the main motion hydraulic motor. Thus, two communications are established: the back 
negative and the direct positive. 

A flow control (see Fig. 1) is on the return line of the hydraulic motor of the main motion of the drilling device 
under consideration. An adjustable choke (see Fig. 2) is on the return line of the hydraulic motor of the main motion of 
the device for processing curved surfaces. This is the fundamental difference between these devices. 

Such variations of the unified modules allow providing a combination of movements of the working body of 
the mobile machine: rotation-feed; rotation-rotation; feed-feed [14]. 

Design differences affect the characteristics of direct positive and back negative communications of the 


devices, both unified and modular (Table 1, 2). 
Table | 


Functional dependences of adaptive communications between the module and the flow control 

















. Functions of module adaptive communications 
Types of hydraulic motors : : — 
Back negative Direct positive 
Hydraulic motor of the main motion 
— feed hydraulic motor Wr = AMg — (BMerg + CMegn) Wry = ky Wry 
Hydraulic motor of the main motion 
Wp, = AMy — (BM,,, + CM Vin = k2w 
—feed hydraulic cylinder o ” 0 — Mer, can) 
Hydraulic cylinder of the main 
Vin = AMy — (BM,,, + CM Vun = k3 Vv 
motion — feed hydraulic cylinder a 0 ~ (BMerg cu) - less 
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Table 2 
Functional dependencies of the adaptive communications of the module with an adjustable throttle 
. Functions of module adaptive communications 
Types of hydraulic motors - - = 
Back negative Direct positive 
Hydraulic motor of the main motion fing |2APing 
Wry = AMy — (BM.,, + CM = — yy Aes [SO aps 
— feed hydraulic motor = 0 ~ BMcrg can) a a igi p 
Hydraulic motor of the main motion fing |2APing 
= AMy — (BM,,, + CF = = RPE [ee 
—feed hydraulic cylinder is 0 ~ (BMerg can) Van = Korg — fun p 
Hydraulic cylinder of the main fing |2APing 
Vun = AMo — (BF erg + CF = — Ae8 Jo aps 
motion — feed hydraulic cylinder me 0 — (BF erg cut) Yun = KeVurg ~ Lae p 














Let us note that the ratio of the velocities of the output elements of hydraulic motors (see Table 2) does not 
depend on the flow rate through the flow control, since this indicator does not depend on the pressure drop on the 
device, i.e. Opn = const. 

Adaptive communications play a special role at action of working bodies of machines in modes with the 
required indicators of speed of performance of operations and stability [14]. The discussed above design and functional 
features of the devices (presence of an adjustable throttle or flow control) allow implementing such communications. 

Tables 1 and 2 provide the legend: v,,, — rod speed; fi, — cross sectional area of the feed hydraulic cylinder; 
W;, — angular velocity of the hydraulic motor main movement; A,B,C — constant coefficients depending on the 
moments of inertia, kinematic parameters and efficiency of the units of the main drive and feed drive module; My — the 
nominal and the calculated total moment of resistance corresponding to the processing without adaptation; fj, — area 
of flow cross-section of the throttle Ips; Men, Meu — the actual moments of resistance taken up respectively by the 
main drive shaft and the feed shaft; Fo.) , Fey, — the actual resistance forces, which are taken up respectively by the rod 
of hydraulic cylinders of main movement and the feed cylinder; AP,,3 — the pressure differential on the throttle J[ps; p 
— density of the working medium in the hydraulic system (oil); k,,k2,k3 — the coefficients of conversion of design 
features and dimensions of the speeds of movements. 

Back negative communication allows stabilizing the actual total moment (force) of the resistance on the 
working body in the process of the module work. In this case, we automatically compare: 

- the actual moment of resistance with the value set by the mode regulator; 

- changes in feed rates and the main motion with the sign, reverse to the error signal. 

The following are the relevant equations of back communication on accelerations of the shaft of the feed 
hydraulic motor. 

If there is a flow regulator: 








d 
oo = 1 | +7p(BAMErg + CAMEyn) (1) 
Kern (BAMEr,+ CAME an) 
If there is an adjustable throttle: 
d 
. ; +a (BAME + CAME) 
ner Kern IC (BAM, oF CAM@,,)°° 


d |{2(APo 44+ 6(AP®,3) 
+Hapsfap3 a J — 5 2)I. (2) 


The methodology of functional unification is determined by its goals and features, as well as the stage of the 
module life cycle. However, we should note the overall steps of its implementation: 

- clarification of the functional purpose and conditions of the unified module operation; 

- preliminary data collection within the framework of constructive and functional unification of hydraulic 
devices of standard and unified module (including single-valuedness conditions; nominal and boundary values of 
operation parameters, especially pressure drops on hydraulic devices under the conditions of dynamic equilibrium of 
hydraulic system); 

- clarification of the adaptive communication type according to the design features of the modules (in 
accordance with Tables 1 and 2). 

In compliance with the mentioned general provisions, the functional unification of adaptive communications is 
performed in a certain sequence. 
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- According to the calculation or by the specification, the passport, we specify the values of working (or 
maximum) loadings on the working body of the mobile machine with the unified module. 

- In relation to the standard module, we determine the functional and design values of the parameters of the 
hydraulic elements of the unified module. Here we use nomograms, graphs or calculations (for example, by the method 
of similarity of functioning) [15]. 

Research Results. When performing unification, the principle of superposition of functions is used. 

Let us briefly consider an example of the method of functional unification of hydraulic motors of the module. 

The output characteristics of the hydraulic motor and the hydraulic cylinder are angular velocity (rotation 
frequency) and torque on the shaft of the hydraulic motor, linear speed and force on the rod of the hydraulic cylinder. 
Functional dependencies that define these characteristics have the following form [16]: 

a) for the hydraulic motor 


M.,, = ff ta ss q, B, J, Mn» ? AY (Oe. cd Tu ? p)), 


Nyy = LO 04,0 (AP n> Sy» P))5 (3) 
b) for the hydraulic cylinder 


V,, = f(O.d,,8,p,AP,,E,V,), 


ul 


i = FAAP AF, ji F ios BE spe BD (4) 


Here n.,, — shaft speed; g — volume of oil; O — flow (supply) of oil; E — elastic modulus; 1 — dynamic oil 
viscosity p — density of oil; P,, — inlet pressure; AP,, — differential pressure on leakage element; O,, — leakage 
value; f;, — the area of the slit through which passes a leak; J — moment of inertia; 8 — oil compression coefficient; ; 
Frp, gr — equivalent force of breakaway and internal friction in the cylinder; r — radius of rotation of the mass; €, v 
— structural coefficients; M, — shaft torque; M,, — friction torque; AP,, — differential pressure on hydraulic cylinder; 
F,, —stress on the rod; AFn — the difference between the forces on the piston. 

When carrying out constructive and functional unification, special criteria of similarity of unification [17] 
obtained by the dimensional analysis [15] are used. 

For a hydraulic motor, the similarity criteria of unification have the form: 


Nons, qd pis. E 








aes y ae uy Mae pte Phe 
Mee AP TP, 
Tvy,, PF Toys a 
(5) 
Ocw’® p ge 
Qe POS Ape P vagl hae : 
Similarity criteria of unification have certain functional or structural sense, namely: 7, — the criterion of 


flow kinematic parameters; 1, — the criterion that characterizes the ratio of constructive and consumable (functional) 
characteristics; ™, — the criterion that characterizes the nominal stiffness of the working fluid; 7, — the theoretical 


(indicator) criterion of the shaft torque of the pump; =, — the criterion of elastic deformation of the working fluid; 


Tam —— the criterion of friction losses in mates; 2, — the criterion of the loss of inertia hydraulic-mechanical 


im 


resistance; ut Q — the criterion of the loss of hydrodynamic resistance of the working fluid; ™,, — the criterion of 
c ’ 


the loss to overcome the leakage of the working fluid; ~,,, — the criterion of the loss on the inner (viscous) resistance 


of the working environment. 

Thus, each criterion reflects the essence of one of the properties of the functional unification of the module, 
and any system of criteria reflects this property generally. 

The generalized similarity criteria for unification obtained by combinations of separate criteria are presented in 
the form of: 


1 wtp Ty Pq?"'rJN,, OP —1 Of,,..f, EAP, 
M,= —_* ;n, =1,0,T, ES 
T,T,T, My MoT, P.ukfy 
Particular similarity indicators of the functional unification of C;, are obtained from the equations (5) for the 
unified model and standard modules [15]. In this case, the respective similarity number zy; and z,,; for the standard and 


the unified modules should be equal: 


(6) 
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Ty i = Tn idem. (7) 
In this example, separate indicators of similarity of functional unification have the following form: 
C C CE Copp C3 
1= C15; 1=— Cris; l= 5 1 =CeCp_ 5 1= : 8 
Cq Su Cq Tu Cap BY Pox PpxJ? : ( ) 
qm Liem a Pg = 50% 083 _ CQ CaP 
CPgy C3” Pox’ nee , Cq? Cp 


As well as similarity criteria of functional unification, separate indicators can be combined into complexes, 
obtaining the effect of superposition. 

The change of scale of parameters C; of the unified module in indicators (8) in relation to similar parameters of 
the indicators of the standard module shall correspond to the condition of equality to 1 [15]. For example, for the output 
characteristic of AP.,, on the indicator 


Ca C. 
ae (9) 
Ca” Cp 


you can analyze the functional unification on the pressure drop, the characteristic (working) volume and flow 
of the hydraulic motor, oil density. 

The required standard size of one of the parameters of the unified module is determined by the substitution in 
(5) and (8) of the known values of the parameters of the standard module hydraulic motor and a priori known 
(established by technical specifications, etc.) values of the parameters of the unified module. 

The features of the procedure of functional unification research depend on the task. When carrying out 
unification of the module of a new standard size, first, you carry out constructively-functional unification, the choice of 
standard sizes of elements of the hydraulic system, and then you check their functional compliance to the internal 
adaptive intercommunications. To solve the first problem, we use particular criteria and indicators of functional 
unification similarity. To solve the second problem, generalized criteria, indicators, as well as the equations of direct 
and back adaptive intercommunications are used. 

Other objectives of the unification of the adaptive module are: analysis of the causes of violations of functional 
unification in the process of operation of a single adaptive module, and the study of the effectiveness of functional 
unification in the modernization of the adaptive module of a hydraulic drive. The method of similarity of functioning of 
technical systems can also be used [18], [19]. 

Discussion and Conclusions. The urgency of development and application of the unified hydraulic drive of 
modular type in mobile machines is proved. This design allows you to adapt the power and kinematic parameters when 
the working body is exposed to fluctuating loads. The method of carrying out constructive and functional unification of 
modules with self-adaptation involves the use of mathematical models of adaptive intercommunications. Optional 
versions of modification of the unified modules with self-adaptation and the corresponding combinations of hydraulic 
motors, regulating equipment and mathematical models of adaptive intercommunications are shown. An example of the 
functional unification technique using the method of similarity of technical systems operation is given. The technique is 
recommended for the study of functional unification in the development of standard sizes and in the operation of the 


proposed type of module. 
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IIporno3suposanue YPOBHA TOT3eCMHBIX BOX MeCTOPOKICHHA WeMeHTHOTO ChIPbH Ha OCHOBe AHHAMHM4eCCKHX 
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' JInmeuxuii rocyzapcrsenueiii Texnwueckuit yHupepcurer, r. JIumenx, Poccuiicxaa Deneparma 


Introduction. The development of a mathematical model for 
the groundwater level of a deposit of cement raw materials 
located in the Zadonian-Yelets aquifer, which is the principal 
domestic water supply source for the city of Lipetsk, is 
considered. Therefore, it is necessary to provide ongoing 
monitoring and to have the possibility to predict the water 
level under the field development. The work objectives are the 
identification and study of a dynamic neighborhood model 
with variable hierarchical neighborhoods of the groundwater 
level that enables to adequately predict value of the water level 
in the examined wells. 

Materials and Methods. The definition of a dynamic 
neighborhood model with variable hierarchical neighborhoods 
is given, differing by time-varying double-level neighborhood 
communications between the first- and second-level nodes. At 
each next discrete instant of time, the neighborhood model 
nodes change their state under the influence of the online 
parameters and node states included in their neighborhood. As 
a subcase, we consider a model with line state recalculation 
functions. Parametric identification of the dynamic 
neighborhood model consists in finding the system parameters 
for each second-level node, and is based on the ordinary least 
squares. 

Research Results. A linear dynamic neighborhood model with 
variable hierarchical neighborhoods for predicting the 
groundwater level in a cement raw material deposit located in 
the Zadonian-Yelets aquifer is developed. The software using 
C++ is developed for the parametric identification and 
simulation of the functioning of the dynamic neighborhood 
model under consideration. It enables to determine parameters 


of the node state recalculation functions for a given structure, 
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Beedenue. CTatba NocBaljeHa paspadoTKe MaTeMaTHYecKoli 
MOJeuIM YPOBHA MO J3CMHBIX BO, MeCTOPOXKeHHA LeEMeHTHOTO 
CbIPbA, PaCMONO%*KeHHOTO B 3aOHCKO-eJIeEHKOM BOJOHOCHOM 
TOpH30HTe, ABJIATOLHMCA OCHOBHBIM HCTOUHHKOM 
XO3AHCTBCHHO-IMTbeBOrO BOWOCHaOxeHHA Topoya JIumenka. 
Ilostomy Ha cTaqMH pa3pa0oTKH MecTOopoxyeHHa 
HeOOXOJHMO MpOBOZHTb NOCTOAHHbIM MOHHTOPHHT HM MMeTb 
BO3MO2%KHOCTb IIpOrHO3HpOBaHHA YPOBHA MO3eMHBIX BOL. 
Lemp padoTrr — upentTudpukalua u uccleqoBanne 
WMHaMHYeCKOH OKpecTHOCTHOM MoyemM Cc TepeMeHHbIMH 
WepapxXH¥ecKHMM OKPeCTHOCTAMH ypOBHA TMO3CMHBIX BOT, 
TO3BOJMIOWeH C BOCTATOUHOM TOUYHOCTbIO MpOrHo3HpOBaTb 
3HayeHHe ypOBHA BO B OOCJIeEAYeMbIX CKBaxKHHAX. 
Mamepuanst u  memoovi. IlpuseqeHo ompeyemeHue 
TMHaMHYeCcKOH OKpeCcTHOCTHOM MoyemM Cc TepeMeHHbIMH 
OTIMMYarIOMleHca 
BYXYPOBHeBLIMH 


MepapxuyecKHMu 
U3MCHAFOWUIMMUCA BO 


OKpecTHOCTAMH, 
BpemMeHH 
OKPeCTHOCTHBIMH CBA3AMH Mey y3laMH MepBoro HM BTOporo 
ypopHa. B kaxkybili clleqyroujMii MCKpeTHbIt MOMeHT 
BpeMeHH y3JIbl OKPeCTHOCTHOM MOeIH MeHAIOT cBOe 
cocTosHHe Mo BO3elCTBHeM TeKYIIMX ympaBieHHuii u 
COCTOAHHM y3I0B, BXOJAMJHX B HX OKpecTHocTH. B kayecTBe 
yacTHOrO cyly4ad paccCMOTpeHa MOjJeIb Cc JIMHeMHEIMH 
dyHKuMAmMH = =epecueta coctosHuli. [lapamerpwueckas 
WeHTHpUuKaLHA WHHAaMHMYeCKOM OKpecTHOCTHOH Moe 
3akIO4aeTCA B HaXOKZCHHH WapaMeTpoB CHCTeMBI IA 
K@kKQOrO y3la BTOporo ypoBHA HM OCHOBaHa Ha MeTOye 
HaMMeHbIUMX KBa[paTos. 
Pesyismamet —succiledoeanua.  Pa3padotaHa _—sIMHeliiHaa 
TMHaMH4ecKkad OKPeCTHOCTHad MOjeyIb Cc TepeMeHHbIMH 
WepapxXHuecKHMH OKPeCTHOCTAMH Id MporHosHpoBaHHua 
YPOBHA TNOZ3eMHBIX BOX MeCTOPOXKACHHA WeMeCHTHOTO CBIpbA, 
pacrosox*KeHHOTO B 3a]OHCKO-eJIeL]KOM BOJ]OHOCHOM TOpH30HTe. 
J\ia mapametpwyueckolt ueHTHpuKaywH UM MOjeIMpoBaHHsa 
(PyHKWMOHMpOBaHHA 
OKpecTHOcTHOH 


TWMHaMvyeckon 
TIporpaMMHoe 


paccmatTprBaemoli 
MOjeuIH pa3pa6oTaHo 
oOeceyeHve Ha sa3bIkKe C++, MO3BOUIAIOMIee JIA 3aaHHOl 


PaOota BEIMOIHeHa Ip PuHAHCOBOL NoAWepxKke POO. (mpoext Nel 6-07-00-854). 
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and also to predict the model behavior in the operation 
process. A hierarchical structure is given, and a parametric 
identification of the linear dynamic neighborhood model of the 
groundwater level is carried out. After the parametric 
identification on the teaching data selection, the mathematical 
model is checked on the test sample. 

Discussion and Conclusions. The obtained average ratio errors 
of the identification and forecast suggest the developed model 
validity and enable to recommend it for predicting the 


underground water level of a cement raw materials deposit. 


Keywords: groundwater level, deposit of cement raw 
materials, dynamic neighborhood model with variable 
hierarchical neighborhoods, parametric identification. 
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CTIpyKTypbl HaxOQHTb apaMeTpbI (yHKMH mepecueta 
COCTOAHHH y3JI0B, a TakoKe TIPOTHO3HPOBaTb MOBeeHHe Moe 
B Tpoyecce dbyHkyHoHupoBanusa. Ipupeyena uepapxwueckad 
CTpykTypa HM UpoBeyeHa MapaMeTpueckad MeHTHpUuKalHA 
JmMHeHHOM = MHaMWYueCKOM OKpeCTHOCTHOH MOjleuIH ypOBHA 
TIOWBeMHBIx Boy. T[locne BbINONHeHHA TNapamMeTpHueckolt 
weHTupukaywH =69Ha =e OGyyalolyeli = BLIOOpKe =——-JaHHbIX 
MaTeMaTH4ecKasd MOJelIb IpoBepeHa Ha KOHTPOJIbHOM BEIOOpKe. 
O6cyacdenue  u Tlomyyenupie cpeqHHe 
OTHOCHTeIBHEIe § ONIMOKH HyeHTHpuKaWMH MM MporHno3a 


3AKINOYeEHUA. 
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TIO3BOJIAIOT PCKOMCHIOBATb Ce JIA TIpOrHOSHpOBaHHA YpOBHA 
TION3CMHBIX BO], MCCTOPO?KCHHA WEMCHTHOYLO CBIPbA. 


Koro4uesbie 
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OKPeCTHOCTHad MOJeIb C MepeMeHHbIMH WepapxuyecKHMu 
OKpecTHOCTAMH, HapaMeTpHyeckad HAeHTHuKallMaA. 


c1OBa: YPOBeHb 
IeMeHTHOTO 


TIO3CMHBIX BO, 


CBIpbA, AMHaMidecKkaAa 
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Introduction. Zadonian-Yelets horizon is the main productive complex of limestone suitable for 
cement and metallurgical industries [1]. Limestone is light and light gray, with a yellowish tinge, medium-hard, 
fractured, fine-grained. You can often find hard limestone represented by silicified varieties. Limestone is often porous 
and cavernous. Clay rocks suitable for cement production lie among the quaternary deposits covering the limestone. 

Overburden rocks are represented by a fertile layer of soil, loam, substandard clay, sand and a layer of crushed 
stone in the upper part of the destroyed limestone. 

Hydrogeological conditions of the deposit are simple. Sokolsco-Sitovsky deposit is in in geomorphological 
terms in the slope part of the valley of the Voronezh River, which determines the hydrogeological situation in the field 
area [2]. 

The aquifer is fed by infiltration of atmospheric precipitation and by absorption of flood flow through the 
valleys of gullies and ravines. Groundwater discharge is in the Voronezh River. 

The current position of the groundwater level was studied by monitoring observations in 7 wells evenly located 
on the field area. The level of groundwater in the low-water period varies from 102.5 m to 109.7 m. 

Currently, the Zadonian-Yelets aquifer is the main source of the domestic and drinking water supply of the 
Lipetsk city. In this regard, the development of limestone can be made only in the water-free part of them with high 
pillar of at least 2.0 m, which is recommended by the sanitary service of the Lipetsk region. 

Water inflow into the developed pit of the deposit occurs only due to atmospheric precipitation. In this regard, 
at the stage of development of the field, it is necessary to constantly monitor the state of groundwater of the Zadonian- 
Yelets aquifer and to be able to predict the groundwater level. 

Dynamic neighborhood models with variable hierarchical neighborhoods are used in this paper to model the 
groundwater level of the field under consideration. These models allow modeling complex spatially distributed 
processes and objects that change their state in time [3-8]. 

Neighborhood models were first proposed in the late 90-ies of XX century [9]. The basic definitions and 
algorithms of the theory of neighborhood modeling are given in [9-12]. The concepts of “neighborhood”, 
“neighborhood communications” are considered in [13-14]. They use agents that move around the neighborhood and 
interact with each other according to certain rules. Today, the theory of neighborhood modeling is being actively 
developed. There appeared dynamic [6-8], non-deterministic models [11], models with variable neighborhoods [15]. 

Materials and Methods. In this paper, we use dynamic neighborhood models to predict the state of spatially 
distributed systems. In them, each node is an independent object, functioning in time and related by neighborhood 
communications to the other objects of the system. This distinguishes the models under consideration from the widely 
used today neural networks, which can be used to simulate the operation of each object or node separately. The method 
of neighborhood modeling is intended for simultaneous joint modeling and prediction of the behavior of all elements of 
a distributed system. 
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Next, we consider the dynamic neighborhood models “input-state” with variable hierarchical neighborhoods. 
They differ in time-varying double-level neighborhood communications between the first- and second-level nodes, and 
in line state recalculation functions. They enable, in comparison with the known one-level neighborhood models, to 
perform forecasting with higher accuracy. 

The dynamic neighborhood “input-state” model with variable [15] hierarchical neighborhoods can be specified 
by the set NS, =(N,X,V,G, X[0],t) , where: 

1) N =(4,0,,0,,0,,) is a two-level structure of the neighborhood model; A= {a,,a,,...,4,} is a set of nodes 
of the first level; O, and O, are neighborhoods of communications of nodes by states and by controls, respectively; 


O,,,. are hierarchical neighborhood communications between nodes. Each node a, ¢A has its own neighborhood 


defined by states O,[a,] CA and controls O [a,]c A; O, = UO, (a), O.= YO, [a,]. 


Some nodes of the first level a, € A are assigned a set of nodes of the second level O,,,[a,] = {a;,...,a°}. All 


nodes a, €O,[a;], a, €O,[a,] have an impact on the second level nodes a €0.fe), 
At each point in time, only a single active node is defined b =1,...,c such that a? € O,,,[a,,¢]. 
The node a? is active at the time f¢ if it meets the specified activation condition Baa tal = true. For all nodes of 


the second level a? O,[a?]=O,[a,]; O,[a?]=O,[a, ]. 


Yr 
2) XeR" is a block vector of states of the neighborhood model in real time, each block of which 


X[a,|= X[i]e R” is a vector of states in the node a, of the system i =1,...,n. 


ym 
3) VeR™ isa block vector of controls in real time, each block of V[a,]=V[i]¢ R” is a vector of controls 


in the node a, ofthe system i=],...,.n. 


4) G:X, xV, +X isa vector function of states recalculating of the neighborhood model, where X, is a 
set of states of the nodes of the first level, included in the neighborhood O, ; V,_ is a set of controls of the nodes of the 


first level, included in the neighborhood O, . 


For the nodes of the first level a, € A, the function G, will be: 
X[t+1,=Gld =>) X[t+1,77]= > GA, (1) 
b=1 b=1 


where G?:X,, OF ee [i’] is a state recalculation function for the second-level node a?. 


For each node of the second level a? € O,,[a,,¢], the function G? in real time in the linear case has the form: 


teri 


Xe+h’]= >) etl? JAB /+ >) ell? AV A+ 220"), (2) 


a;€Oy [ay ] 4 €0, [ay] 
where a,,a, € A (ik = Length) are nodes of the first level of the model; X[t,i’]¢R” is the state in the node 


a’ at time ¢; V[t,i’]€ R” is the control in the node a? at time ¢; g°[i’,jJeR”””’, g?[i’,k]eR™™, gf’ ]eR’™ are 
model matrix-parameters. 
For each second-level node a? ¢ O.,,[a,,t], the function G? is currently zero, that is X[¢+1,7’]=G?[¢]=0. 


jeri 


va 
5) X[0]<¢R“ 1s the initial state of the model. 


6) ¢ is the current discrete time of the model. 

The structure of the neighborhood model can be represented as a single two-level graph of the structure of the 
neighborhood model-oriented graph with two types of arcs: state and control actions — or two-oriented graphs — 
external and internal structures. 

Parametric identification of the dynamic neighborhood model [16-17] is in finding the system parameters for 
each node of the second level, and it is based on the least squares method: 


E=>)|Xt+1i)-GJs]] > min. 
i=l 
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Research Results. Let us consider the linear dynamic neighborhood model of the level of cement raw 
materials underground waters deposits. Fig. 1 shows the graph of the external structure of the model. Node a, is 
external environment. Currently, 7 wells are used for groundwater level monitoring purposes, which in Fig. 1 
correspond to the nodes a,—a, . The control actions of the nodes V[t,7]<¢ R'° consist of the amount of precipitation in 


mm and of the average daily air temperature in °C for the last 5 days before the measurement of the water level in the 
wells at time ¢. The states of the nodes X[t,i] € R are equal to the water level in the well i at time? , 7 =2,...,8. 





Fig. 1. Graph of the outer structure of the neighborhood model 


[4,]={a',a?} and 


correspond to the positive and negative average daily air temperature. Their structure is shown in Fig. 2. 


V{i| 


The first level nodes of the neighborhood model a, (i=2,...,8) are hierarchical: O 


ier 





Fig. 2. Graph of the inner structure of the first level nodes 


The system (2) for each node of the second level a? of the surrounding groundwater level model will have the 
following form: 
XUe+L?]= gh P+ lh Wb?) + 2), (3) 
where g°[i,i]—eR; g°[i,i]eR”"; ge? [i]eR; i=2,...,8; b=1,2. 
The program in C++is developed for the parametric identification and performance simulation of the 
considered dynamic neighborhood model. The initial data for identification are the structure and the training sample. 


The program allows you to find parameters of the recalculation functions of node states, as well as to predict the 
behavior of the mathematical model in the process of operation. 
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After the parametric identification on the training data sample, the obtained model was tested on the control 
sample. The normalized initial and model values of the groundwater level for the node a, for training and control 


sampling are shown, respectively, in Fig. 3 and Fig. 4. 


=—@=X[2] —=X"(2 


Fig. 3. Normalized source and model values of the groundwater level for the node a, 


=—X[2] = x*[2 


Fig. 4. Normalized source and predicted values of the groundwater level for the node a, 


The average relative error of identification (prediction) of the neighborhood model is calculated by the 
formula: 
MX [t+Li]-X,[t+Li 
ae 1 >» mf + Li] mie 4! 100%, 
Mn X,[t+Li] 


k=1 i=l 





where X,,[f+1,i] is the state of the node in the m-th training (control) sample tuple; Xx [iri i] are model 
values of the node state a,; M is the volume of the training (control) sample. 


The average relative error of prediction is shown in Fig. 5. 
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Fig. 5. Average relative prediction error for each tuple of the control sampling 


The average relative error of identification was 0.19%, the average relative error of prediction was 0.23%, 
which indicates the adequacy of the developed model and allows us to recommend it for predicting the level of 
underground waters of the cement raw material field. 

Conclusion. A linear dynamic neighborhood model with variable hierarchical surroundings to predict 
groundwater level deposits of cement raw materials, located in Zadonian-Yelets aquifer, is developed. 

The software in C++programming language is developed for the parametric identification and performance 
simulation of the considered dynamic neighborhood model. 

The implemented model predicts the level of groundwater in the analyzed wells with sufficient accuracy and 
can be effectively used to predict the level of groundwater deposits of cement raw materials. 
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Introduction. Production processes quality depends largely on 
the management infrastructure, in particular, on the information 
system (IS) effectiveness. Company management pays increas- 
ingly greater attention to the safety protection of this sphere. 
Financial, material and other resources are regularly channeled 
to its support. In the presented paper, some issues on the devel- 
opment of a safety enterprise information system are consid- 
ered. 

Materials and Methods. Protection of the enterprise IS consid- 
ers some specific aspects of the object, and immediate threats to 
IT security. Within the framework of this study, it is accepted 
that IS are a complex of data resources. A special analysis is 
resulted in determining categories of threats to the enterprise 
information security: hacking; leakage; distortion; loss; block- 
ing; abuse. The connection of these threats, IS components and 
elements of the protection system is identified. The require- 
ments of normative legal acts of the Russian Federation and 
international standards regulating this sphere are considered. It 
is shown how the analysis results enable to validate the selec- 
tion of the elements of the IS protection system. 

Research Results. A comparative analysis of the regulatory 
literature pertinent to this issue highlights the following. Differ- 
ent documents offer a different set of elements (subsystems) of 
the enterprise IS protection system. To develop an IS protection 
program, you should be guided by the FSTEC Order No. 239 
and 800-82 Revision 2 Guide to ICS Security. 

Discussion and Conclusions. The presented research results are 
the basis for the formation of the software package of intellec- 
tual support for decision-making under designing an enterprise 
information security system. In particular, it is possible to de- 
velop flexible systems that allow expanding the composition of 
the components (subsystems). 


“The research is done with the financial support from the Russian Federation President Council on Grants within the frame of R&D “Building a model of 


Beedenue. KauecTBo MIpOu3BOACTBeHHBIX IIpOleccOB BO MHO- 
TOM 3aBHCHT OT HH(PpacTpyKTypbl yupaBseHHA — B 4acTHOCTH, 
OT 9exTHBHOCTH WHPopMalMoHHOH cuctempr (UC). Me- 
He]~KMeHT KOMHaHHi yelieT Bce OobUIee BHUMaHHe OObectie- 
yeHHIO Oe30MacHOCTH 9TOM ccbepbl, Ha ee NOAAepxKKy perysap- 
HO HalipaBJIAioTcd (MHaHCOBbIe, MaTepHasIbHble WM ~Apyrue 
pecypcpt. B nmpeyctaBieHHol paOoTe paccMOTpeHbI BOTIpPOCch! 
TIOCTpOeCHHA KOMIMIeKCa 3alMTbl HHPOPMalMOHHOM CHCTeMBI 
TpeqnpuaTus. 

Mamepuaaei u memoooi. Oxpana UC npegnpusTua yauTbIBaeT 
ocoOeHHOCTH OObeKTa 3aLIMTbI HW akTyaIbHble yrpo3bl HHop- 
MaljHOHHOM Oe3onmacHocTuH. B pamMKax aHHoro MccreqoBaHHA 
mpHusato, uo MC mpegctapnsaet coOol KommeKc HHPopMaln- 
OHHBIX pecypcos. Ilo pe3ymbTaTaM cleljHabHOrO aHasIM3a 
ompeyeseHbI KaTeropHu yrpo3 HHpopMalHoHHON Oe30nacHo- 
CTH IIpeAMpHATHA: B3JIOM; yTeuKa; HcKaxkeHHe; yTpata; OOKH- 
poBaHue; 30ynoTpeOmeHHe. BaiaBieHa CBA3b JaHHBbIX yIpo3, 
kKomnoHeHToB MC uv sieMeHTOB KOMIIeKca 3alHTbI. PaccmoT- 
peHbI TpeOoBaHHa HOpMaTHBHO-IIpaBOBbIX akTOB PoccuiicKoli 
@Meyepaywu HU Mex] yHAapOAHBIX CTaHapToB, perysIMpyloulux 
qanHyto cibepy. Moxa3zaHo, kakuM o6pa30M pe3ysIbTaTbl JaHHO- 
TO aHasI3a MO3BOJIAIOT OOOCHOBaTb BbIOOp 3ICMCHTOB KOM- 
Tiekca 3amuTbI MC. 

Pe3yiemamel ucciedoeanus. CpaBHUTeIbHbIM aHasiM3 persa- 
MeHTHpyIolel IMTepatTypbl, OTHOCAINelica K aHHOMY BOTIpO- 
CY, HO3BOJIMJI BBIABHT cyleqyroulee. PasHbie JOKYMeHTBI ped - 
WaraioT pa3HbIii HaOop 9dJIeMeHTOB (HOACHCTeM) KOMMIeKca 
3amjuTE MC npegnpuatua. PaspaOaTpiBas WporpaMMy 3allHTbI 
UC, cnenzyer pykosoycrBosatsca IIpuxka3om WCTOK Ne 239 u 
cTaHyapToM 800-82 Revision 2 Guide to ICS Security. 
O6cyacdenue u 3axmouenua. Pe3ynbTaTbl MmpecTaBsIeHHoro 
HCCHeNOBaHHA ABIAIOTCA OCHOBOM Id (OpMupoBaHHA Tmpo- 
TpaMMHOTO KOMIVJIeKca HHTeIICKTYAaIbHON NOAWepxKKH NIpHHsA- 
THA pelleHHi pH MpoeKTHpOBaHH CHCTeMBI 3all[HTbI HHop- 
Mallu Ha pezmpuaTuu. B uacTHocTH, MoxHO pa3paOaTEIBaTb 
THOKHe KOMIVICKChI, MO3BOJIAIOMIHe paclIMpATb COCTaB 9JIEMeH- 


TOB (11oJIcHcTeM). 
oe | 


intellectual support for decision-making when designing an enterprise information security system”. 
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Introduction. Information systems (IS) are increasingly used in the production and management processes. Accord- 
ingly, the problem of cyber security (CS) of IS gets worse. In particular, IS weak isolation simplifies an unauthorized ac- 
cess to them [1, 2, 3]. The consequences of fraudulent attacks on IS may be production downtime, financial losses, and on 
the worst-case scenario — even man-made disasters [4]. Thus, the crucial task is to establish a protective system for the in- 
dustrial IS which could effectively militate against malicious acts. 

Materials and Methods. The development of an information security system is based on the results of a pre-project 
study during which the setup of the asset to be protected and the immediate threats are determined. 

The asset of protection is represented as a set of information resources: 


Object... = {NE,CC, IS, Sts, WS, PE, OS, SS, AS, IP,Sn, RSM, SM, IA}. 

Here, NE is a set of network hardware; CC is a number of communications lines; JS is an array of infrastructure serv- 
ers; Sts a set of data storage systems; WS is number of localhosts; PE is a number of external equipment; OS is an 
array of operating systems; SS' is a set of system software; AS is a set of application software; JP is number of infor- 
mation processes running in the tech companies; Sm is subnets; RSM is a number of removable media; SM _ is electronic 
data storage devices; JA is information assets. 

Number of immediate CS threats Threat is determined [5]: 

Threat = { Breaking, Leak, Distortion, Loss, Blocking, Abuse} . 

Here, Breaking is hacking threats; Leak is information leakage threatening; Distortion is threats of distortion; Loss 
is threats of loss; Blocking is threats of blocking the information resources of a tech company; Abuse is risks of abuse. 

The system which meets these threats is an enterprise IP protection system (ISPS) (Fig. 1). 


ad RSM 


SM 1A 


enterprise ISPS 
SPI 





Fig. 1. Security objects — threats relations in ISPS diagram 
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SPI (system of protection of information) is two-tiered, and it includes subsystems (components) [6]: 
— a set of subsystems (Subsystem) of information security; 
— number of information security facilities (MP, means of protection). 

Fig. 2 shows a general form of the ISPS structure. 





Fig. 2. Generalized structure of enterprise ISPS 


When determining the components of the information protection system, experts proceed from the analysis of 
available regulatory documentation and standards operating at the enterprise. International experience should be taken into 
account as well. 800-82 Revision 2 Guide to Industrial Control Systems (ICS) Security [1] is widely used in the foreign and 
domestic practice. This standard was developed by the US National Institute of Standards and Technology. It, in particular, 
contains recommendations on improving safety in the industrial inspection systems, including the supervisory control and 
data acquisition systems. It shows how the organizational processes and business functions are subjected to threats, and it 
describes usual vulnerabilities. Special attention is given to security measures and counters that should be undertaken in a 
hostile situation. 

Research Results. Domestic regulatory legal acts (RLA) handling the enterprise IS protection can be conditional- 
ly divided into two categories [6]: 

—RLA on maintenance of information safety of the automatic process control system (APCS); 
- RLA on the security of critical information infrastructure (CII). 

Crucially, vulnerabilities in the CII protection can cause major material and environmental damage. Inadequate 
CII protection is fraught with social and military-political problems. 

Designing the information security system (in particular, when modeling the intellectual support for decision- 
making) requires a preliminary comparative analysis of the profile regulatory legal acts of the Russian Federation. As for 
example, the following documents should be considered: 

— Order no. 31 of the Federal Service for Technology and Export Control (FSTEC of Russia) of March 14, 2014, “On 
Approving Requirements for Providing Information Protection in Automated Control Systems of Production and Techno- 
logical Processes at Critical Objects, Potentially Hazardous Facilities, and Objects posing high threat to life and health of 
people and the environment” [7]; 

— Order no. 239 of FSTEC of Russia of December 25, 2017, “On Approval of Requirements for Ensuring the Security of 

Significant Facilities of the Information Infrastructure of the Russian Federation” (draft) [8]; 
— International standard 800-82 Revision 2 Guide to Industrial Control Systems (ICS) Security [9]. 
Comparative analysis of the Orders of FSTEC of Russia no. 31 and no. 239 is presented in Table 1. 
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Table 1 


Elements (subsystems) of enterprise IS protection system in 
Orders of FSTEC no. 31 and no. 239 





FSTEC Order no. 31 of 14.03.14 FSTEC Order no. 239 of 25.12 2017 








Identification and authentication of access subjects and access objects (IAF) 





Access control of access subjects and access objects (UPD) 





Restriction of software environment (OPS) 





Protection of machine-readable media (ZNI) 





Secure event logging (RSB) Security audit (AUD) 





Antivirus protection (AVZ) 


Intrusion detection (OV) Intrusion prevention (PV) 


Protection of information (automated) system and its 








Infosecurity control (analysis) (ANZ) 
components (ZIS) 








Integrity control (OTsL) 





Information assurance (ODT) 





ISPS Subsystems 


Event planning to ensure information protection (PLN) 





Protection of technical facilities and systems (ZTS) 





Security assurance of software development (OBR) Information security incident response (INTs) 





Virtualization environment protection (ZSV) Personnel informing and training (IPO) 








Software update control (OPO) 





Security protection of emergency procedures (DNS) 





Analysis of threats to information security and risks from 





their implementation (UBI) 





Configuration control of data processing system and its security system (UKF) 











Note. For greater clarity, variances are not only placed in different cells, but are also highlighted in gray 








So, FSTEC Order no. 239 provides for the following subsystems in the IS protection system: 
— security audit (AUD); 

— protection of information (automated) system and its components (ZIS); 

— information security incident response (INTs); 

— staff informing and training (IPO). 

It should be mentioned that the decision on the ISPS completeness to a certain extent depends on the financial ca- 
pacities of the enterprise. However, if the cost of the protected resources and the potential damage from hazards is higher 
than the cost of the ISPS, then it makes sense to implement AUD and ZIS. 

Comparative analysis of the FSTEC Order no. 239 of December 25, 2017, and 800-82 Revision 2 Guide to Indus- 


trial Control Systems (ICS) Security is given in Table 2. 
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Table 2 
Elements (subsystems) of enterprise IS protection system in 
FSTEC Order no. 239 of December 25, 2017, and in 800-82 Revision 2 
Guide to Industrial Control Systems (ICS) Security) 





FSTEC Order no. 239 of 25.12.17 800-82 Revision 2 Guide to ICS Security 
Identification and authentication (IAF) 
Identification and authentication 
Access control (UPD) 
System and communications protection, 
Security assessment and authorization 
Restriction of software environment (OPS) System and information integrity 
Protection of machine-readable media (ZNI) 
Media protection 
Security audit (AUD) 
Auditing and accountability 
Antivirus protection (AVZ) System and information integrity 





























Intrusion prevention (hacking) (SOV) Bytom audi iomnatominiehy 





Protection of information (automated) system and its System and information integrity 
components (ZIS) 
Integrity control (OTsL) System and information integrity 











Information assurance (ODT) 
System and services acquisition 





Event planning to ensure information protection (PLN) 
Planning, contingency planning 
Protection of technical facilities and systems (ZTS) 
Maintenance 
Information security incident response (INTs) 
Incident response 
Personnel informing and training (IPO) 
Personnel security 
Organization — wide information security program 
management controls 
Security protection of emergency procedures (DNS) 
Physical and environmental protection 
Awareness and training 
Configuration control (UKF) 
Configuration management 
_ Risk assessment 


ISPS Subsystems 














Software update control (OPO) 























= System and communications protection 











Note. For greater clarity, variances are not only placed in different cells, but are also highlighted in gray. 





In this case, the following differences are most obvious: 
— 800-82 Revision 2 Guide to ICS Security combines the functional of the subsystems of OPS, AVZ, SOV, OTsL, ZIS 
defined in the FSTEC Order, in the subsystem of “System and information integrity”; 
— 800-82 Revision 2 Guide to ICS Security provides a subsystem for the protection of communication systems — “System 
and communications protection”; 
— The FSTEC Order combines the functional of the subsystems “System and communications protection” and “Security 
assessment and authorization” in the “Access control subsystem (UPD)”; 
— The FSTEC Order combines the functional of the “Planning” and “Contingency planning” subsystems in the subsystem 
of “Event planning to ensure information protection (PLN)”; 
— The FSTEC Order combines the functional of the subsystems of “Physical and environmental protection” and “Aware- 
ness and training” in the subsystem of “Security protection of emergency procedures (DNS)”. 

Special mention should be made of the subsystems for risk assessment and protection of communication systems 
[10]. These are the critical items of the enterprise IS protection system among those that are not provided for by domestic 
regulatory documents. Their implementation will enable to enhance protection; to respond promptly to incidents that arise 
in the enterprise IS; to counteract attacks timely and accurately. 


Discussion and Conclusions. The analysis results of the IS protection system components will be used to build a 
model of intellectual support for decision-making when designing the ISPS. Particularly, it is planned to foresee the possi- 
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bility of expanding the ISPS subsystem setup. The selection of the system units will depend on the risk assessment, the ex- 
tent of potential damage by injurious actions, and the cost of the ISPS components. 
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AuddepenuupoBanne NOJMHOMOB HeCKOJIBKHX TepeMeHHbIX Ha NoaMH Tasya He4yerHoi MOWHOCTH 
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Introduction. Polynomials in several variables over Galois 
fields provide the basis for the Reed-Muller coding theory, and 
are also used in a number of cryptographic problems. The 
properties of such polynomials specified over the derived 
Galois fields of fuzzy cardinality are studied. For the results 
obtained, two real-world applications are proposed: 
partitioning scheme and Reed-Muller code decoder. 

Materials and Methods. Using linear algebra, theory of Galois 
fields, and general theory of polynomials in several variables, 
we have obtained results related to the differentiation and 
integration of polynomials in several variables over Galois 
fields of fuzzy cardinality. An analog of the differentiation 
operator is constructed and studied for vectors. 

Research Results. On the basis of the obtained results on the 
differentiation and integration of polynomials, a new decoder 
for Reed-Muller codes of the second order is given, and a 
scheme for organizing the partitioned transfer of confidential 
data is proposed. This is a communication system in which the 
source data on the sender is divided into several parts and, 
independently of one another, transmitted through different 
communication channels, and then, on the receiver, the initial 
data is restored of the parts retrieved. The proposed scheme 
feature is that it enables to protect data, both from the 
nonlegitimate access, and from unintentional errors; herewith, 
one and the same mathematical apparatus is used in both 
cases. The developed decoder for the second-order Reed- 
Muller codes prescribed over the derived odd Galois field may 
have a constraint to the recoverable error level; however, its 


use is advisable for a number of the communication channels. 


* The research is done within the frame of independent R&D. 
“ E-mail: vl.deundyak@gmail.com, nadezhda.mogilevskaia@yandex.ru 
PaOota BBINOJIHeHa B paMKax HHHWMaTHBHOH HUMP. 
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Beedenue. TlonwHOMbI HECKOJIBKHX TIepeMeHHBIX Hay NOJLIMH 
Tastya JlexkaT B OCHOBe TeOpHH KoZ0B Puya-Masuiepa, a TakoxKe 
MCIONB3YIOTCA B pale KpHutTorpaduyeckux 3aqay. B padote 
W3y4aloTca CBOHCTBa TaKHX TOJIMHOMOB, 3aaHHbIX Hal 
TIIpOH3BOJbHBIMH MoM Taya HeyeTHoli MouHOcTH. [1a 
TIOMYYCHHBIX Pe3yIbTATOB MIpeANO*KeHbI WBa WpakTH4YecKHX 
IIPHJIOXKeHHA: CXCMAa pa3esIeHHA JaHHbIX HU JeKOWep KOOB 
Puya-Masuiepa. 

Mamepuanot u memoovi. C Uctonb30BaHvem JIMHeHHOl 
anreOpbl, Teopuu nome Tanya u oOuleii Teopuu MOMHOMOB 
HeCKOJIBKUX TepeCMeHHBIX MOJYYCHbI Pe3yJIbTAaTbl, CBAZAHHbIC 
Cc QuddepesunpoBanveM MW HHTerpHpoBaHveM TOMHOMOB 
HeCKOJIBKHX TepeMeCHHbIX Hay ouaMH Taya HeuerHoit 
mMouHocTH. iif BeKTOpOB HWocrpoeH HM W3yyeH aHasor 
omepatopa qu depesunpoBanns. 

Pesyiesmamel uccaedoeanus. Ha ocHoBe MOJsIyY4eHHBIX 
pe3ylbTaToB oO AuddepenyMpopanHu U MHTerpHpoBaHun 
TIOMMHOMOB TipeyWIO2%KeH HOBbIM Wekolep ia KoWOB Puya- 
Mannepa Broporo nopsaka 4 mpeqioxKeHa cxema 
opraHv3alMH pa3qelneHHoH TMepeqayH KOHUAeCHUMAIbHBIX 
WaHHbIX, T.e. TakKOM CHCTeMbI CBA3H, B KOTOpOH MCxXOJHBIe 
aHHble Ha CTOPOHE OTIMpPaBUTeILA paseMEOTCA Ha HECKOJIBKO 
yacteli MW, He3aBHCHMO pyr OT pyra, Wepelarorca 10 
pa3IM4HBIM KaHaJiaM CBA3H, a Ha CTOpoHe MoOyyaTesla v3 
IIPHHATHIX YacTeli BOCCTaHaBJIMBaKOTCA MCXOJHbIC aHHBIe. 
Oco6OeHHocTBIO MpefaraeMOl CXeMBI ABIIACTCA TO, YTO OHA 
TO3BOJIACT 3allHWaTb jjaHHble, KaK OT HeJIerHTHMHOrO 
WOcCTyla, Tak MW OT HelpeyqHaMepeHHBIX OMIHOOK, Ip 9TOM B 
oOonx cilyyaix UcHOMb3yeTCca OMH UM TOT xe 
MaTeMaTHYeckHHi aluapaT. PazpadOoTaHHbli feKoyep aA 
konoB Puya-Mannepa Broporo mMopsyka, 3aaHHbIxX Hay 
TIIPOH3BOJIbHBIM HeYeTHBIM TWoOseM Tastya, MO%xKeT MMeTb 
HekOTOpoe OrpaHHyeHHve NO YHCIy HCHpaBJIAeMbIX OWIHMOOK, 
oWHakO, eFoO HCMONb30BaHHe WeslecooOpa3sHo Wa pana 
KaHasIOB CBA3H. 
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Discussion and Conclusions. The proposed practical 
applications of the results obtained are useful for the 
organization of reliable communication systems. In future, it is 
planned to study the restoration process of the original 
polynomial by its derivatives, in case of their partial distortion, 


and the development of appropriate applications. 


Keywords: polynomials in several variables, Galois fields, 
polynomial derivatives, differentiation of polynomials, Reed- 


O6cyarcdenue u 3axiouenun. UpeqmoxeHHble mpakTHueckue 
TIpHO%KeHHA TOUYYCHHBIX Ppe3YIbTAaTOB [1peCTaBJIAKOTCA 
TIOe3HbIMH JIA OpraHw3allM HajlexKHBIX CHcTeM cBa3H. B 
uccyleqoBanHe = mporecca 
BOCCTAHOBJICHHA HCXOJHOFO NOJIMHOMA HO erO NPOH3BOHBIM, 


TasIbHenIeM TlanHupyetca 


B Cily4ae UX YaCTHYHOrO HCK@KCHHA, U pa3paboTKa 
COOTBETCTBYFOINHX TIpHoKeHHH. 


KunroveBbie CJ10Ba: TOJIMHOMbI HeCKOJIBKHX TepeMeHHBIX, 
noma Tanya, Mpou3BOAHbIe MOMHOMOB, JHddepeHuMpoBaHve 


Muller codes, decoding, partitioned data transmission. TIOMHHOMOB, xkoybI Puya-Mamepa, WekoqupoBaHnue, 
pa3yleleHHas Mepesaya aHHBIX., 
For citation: V. M.Deundyak, N. S.Mogilevskaya. O6opazeu Ona yumupoeanua: Hleyuaak, B. M. 


Jludepenyupowanne TMOUMHOMOB HECKOJIBKHX TepeMeHHBbIxX 
Hay nomamu Tamya HeaeTHOH MOLJHOCTH HM MpHoxKeHHA K 
Koyam Puya-Mamnepa / B. M. JleyHaax, H. C. Morumescxaa // 
Becruuk J[on. roc. TexH. yH-Ta. — 2018. — T. 18, Ne3.—C. 
339-348. https://doi.org/10.23947/1992-5980-2018-18-3-339- 
348 


Differentiation of polynomials in several variables over Galois 
fields of fuzzy cardinality and applications to Reed-Muller 
codes. Vestnik of DSTU, 2018, vol. 18, no.3, pp. 339-348. 
https://doi.org/10.23947/1992-5980-2018-18-3-339-348 


Introduction. Polynomials of several variables over Galois fields and their derivatives are used in several 
information security domains. Some issues on the integration and differentiation of polynomials of several variables are 
considered in a number of papers. Thus, in [1], polynomials defined over simple Galois fields are studied, in [2-4], 
results for Boolean functions are obtained, and in [5-6], results for polynomials given over ternary Galois fields are 
described. 

In this paper, we consider polynomials of several variables specified over the derived Galois fields of fuzzy 
cardinality. For such polynomials, the results related to the calculation of the directional derivatives, as well as to the 
restoration of the polynomial from the set of its derivatives calculated in the basic directions are obtained. For the 
results, two possible practical applications are proposed: a partitioning scheme and a decoder for Reed-Muller codes 
(PM codes). 

The partitioning scheme can be used to organize the partitioned transmission of confidential data, i.e. such a 
communication system in which the initial data on the sender is divided into several parts and, independently of one 
another, transmitted through different communication channels, and then, on the receiver, the original data is restored 
from the collected parts. The proposed scheme feature is that it enables to protect data, both from the nonlegitimate 
access, and from unintentional errors. At this, in both cases, one and the same mathematical apparatus associated with 
PM codes and polynomial differentiation is used. Partitioned transmission can be used both to improve the 
communication speed, and to ensure data security by complicating the task of interception from several communication 
lines. Some issues on data partition are considered in [7-11]. 

For PM-codes of the second order, deterministic decoders are known only for certain values of the Galois field 
cardinality g . Such as, quite a lot of decoders are known for the case ofg =2, for example [12-13], for the case of 


q=3 and the use of a semicontinuous communication channel, a decoder [5] is constructed. In [14], a second-order 


Reed-Muller decoder for the codes preset over Galois fields of cardinality 2, 4, and 8 is proposed. The second-order PM 
code decoder defined above the derived fuzzy Galois field proposed in this paper is based on reduction to first-order 
Reed-Muller codes. Their codewords can be decoded by any suitable decoder. In the case of PM codes specified over 
the fields with cardinality of more than three, the proposed decoder has some limitation on the recoverable error rate. It 
should be mentioned that the use of the proposed decoding scheme in case of field cardinality of more than three, may 
be advisable at a low noise pollution level of the communication channels used despite the limitation. 

Differentiation of polynomials in several variables 


Suppose g =p’, where p is an odd prime, seN, F, is a Galois field of q cardinality. Let us consider a 


polynomial ring in m variables F,[x,,...,x,,] over F, field. Let us denote the polynomial linear space of F.[%,,...,,,] 
degree not more than r as F isco, . Suppose F” is m_-linear space over F, . 
The result of the differentiation operator action is the polynomial derivative f € gene ] alongh € ae 
(D(X) =f, )-f(®), xe", (1) 


where f, (x) = f(x +b). It is easy to see that D.fe ey cee ], and the operator 





Deundyak V. M., Mogilevskaya, N.S. Differentiation of polynomials in several variables over Galois 





RO) (r-1) 
Dy fF [Xys--+s X pq | ET [2+ Xn (2) 
is a linear one. Let us denote the vector coordinate sum @ € F”’, where p is prime, as a sum of natural numbers, by p(@) 


. The polynomials f € Fiat ] may be written in a canonical form 


{MDS He td axe Dae (3) 


GeF; p(G@)=1 p(@)=2 
where, when writing the monomial ¥% = x!"...x°" , the indices a; will be identified with the elements of the field rs 


and the summands in each sum will be arranged in an lexicographic order of magnitude. If the last sum in (3) is zero, 


‘ : ‘ (1) 
we obtain a polynomial in F’’[x,,...,x,, ]. 


Lemma 1. Suppose g = p*, p is an odd prime, f(x) € Lae eee ] is a polynomial in a canonical form (3), 





b =(b,....b,,) €E” . Then 
S(x)= Fov..00 + He Tracysd nines doe y’ +XAx! > (4) 
(D, f(x) = AG ee oe en ee y’ +2x%Ab" +bAb" =2xAb" + f(b) = Fon.b0 ’ (5) 
where 
Fr00..00 Fi10.00 /2 Fior.00 ID aus Fioo..10 /2 Fivo.01 /2 
Fi0.00 /2 Sor0.00 Sort.o0 fDi Soto.10 /2 Soro.01 /2 
Du Fio1.00 /2 Sor.o0 /2 Soor..00 ree Soot..10 /2 Soot.o1 /2 
Fio0..10 /2 Soto.10 /2 Soor..10 ae Fooo.20 Sooo.11 /2 


Fivo.01 /2 Soro.01 /2 Sooi.o1 /2 a Sooo..1 /2 Fooo.02 


Proof. In the case of a prime Galois field, the proof is contained in [1]. Using (1), (4) and matrix symmetry, we 
obtain: 


fe @)= f(E+B) = fiy.o0 +E +b) (Koo..002 Sor..009s Sun,o1) +E +B)AG +b)", 
(D; f\(%) = fg) — Ff &) = 8 (fo.o0»S0r..008fo.o1). + AB" + BAX +bAB" , 
(D; f(x) — Geers eee ee y +2xAb"' +bAb’ =2xAb* + f(b) foo..o0 eo 


Let us prove the theorem that determines the restoring method up to a polynomial constant term in 





F [x,,2;5-.5X,, | by the set of its derivatives calculated in the basic directions. 


the polynomial f € Ee Eee | of (4): 
= = a ——— 
FS (&) = foo..o0 + * ( fo, bast ons alent) +XAX". 


If 
{(D;/) (X) = 04%, +O4x, +...4+0),%p +o) 7 a (6) 
Then, 
a! oa? a™)(b! ob? Bm)’ 
a=4 G Wy ae Oe ‘ . ae a . (7) 
ie ee ON ke, 
fir.w) ( 0 -BAB" \() by) 
fafa [a a) ° 
Foor) \0g — 5,45,” Jb” by Be 


Proof. From (5), (6), we obtain: 
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. . [oy i i iy? L a et 
Vi=l,...,m: 2Ab/ =(a),0),...,.0,) , £(B)—Soo.00 = M+ (9) 
Then, 

1 2 m 1 2 m 

bob oe B tog? a! 

1 2 m 1 2 m 

oA b, by oe B | 2 Oy 

1 2 m 1 2 m 

b,, b,, _ ais Qin On _ Qin 


Hence, the formula (7) comes out right. 


It follows from (4) that for any b € ig 
f(b) —Soo.00 = BU Tin scuachovaoat 
Forb , we take the vectors b ,< 8 and use the equation f(b,)— fio 9 =), from (9). Then, 
Vi=l,...,m: OL, —b Abt =b A Finwsdonaeee dooce.) : 


>Soo..01 y’ +bAb" 5 


Consequently, 
a —5ABT ) (Bi BL BLY frocon 
ay —b, Ab,” |_| bP By BFW for.oo 
a ~ b,, Ab," bi" b _ o Foo..01 


And the formula (8) is proved. 
Reed-Muller q-ary codes RM, (r,m) . Let us consider RM-codes over the finite field F, where g =p’, p is 


an odd prime, s¢N [15-16]. The elements ak eee aml are information polynomials of the code RM,(r,m) ; we 


suppose thatm=>r>0, m2=2. Vector f made up from the information polynomial coefficients is called an 


information vector. 


In the vector space F”, let us fix some ordering 


{Q,, ..., OL, } (a, = (a1, ©. j9.920+ O ip )) ,n=q". (10) 
The arbitrary information polynomial f(x) € Fx, y+++5%,,] is coded by its evaluation at the ordered space F” points: 
CL) = (FG )o-5 H(G,)) » (11) 


and so, the coder-operator is defined 
G gE pene: 


m 


|=E. 
Reed-Muller codes are defined by the natural parameters r and m (r<m) 
RM, (r,m) = {(f(O Joes (,)) |) € ED 55% | deg(/) <r} cE, 


The parameter r is called a code order. They form a family of linear [n, k, d] q-codes whose length and dimension are 





determined from the formulas 
lig) 
n= q", k i > > (-l)’ CZ CU cd 9 
i=0 j=0 
where [ | is rounding up to the smaller whole number, and the minimum code distance d of the code RM,(r,m) is 
convenient to calculate using the dual code RM, (r*,m) parameters where r+ = m(q—1)—r—1. Suppose p is residue of 
division r* +1 byg—1: r*+1=0(q—-1)+p wherep<q-—l, then the parameter d of the code RM,(r,m) is given by 
the expression 
d=(p+lq°. (12) 
Note that the arbitrary [n, k, d] ,-code enables to correct t = | (d -1)/ 2 | errors in one codeword [17]. 
Next, let us consider PM-codes of orders | and 2 given over Galois fields of fuzzy cardinality, write the 
corresponding information polynomials in the form (3), and use ordering (10) for numbering the information vector 


coordinates. 
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Lemma 2. Letr ¢ {1,2}, q 23, then the minimum code distance of the code RM,(r,m) is calculated by the 
formula: 
d,=(q-r)q"", (13) 
and values of the recoverable t, = | (d —1)/ 2 | - by the codes RM,(r,m), r € {1;2}, are related as follows: 
t,/2<t,. (14) 
Proof. Let us make use of the fact that <q and calculate o and p - subquotient and the remainder of the 
division r~ +1 by g—1 respectively: 
o =(m(q—-1)—r)div(qg-l) =m-1; 
=(m(q-1)-r)mod(q—1) =q-1-r. 
Then, from the formula (12) we obtain (13). From the equalities d, =(q—-l)q”', d, =(q-2)q"', we get that the 


] (q-Nq"" = (q- 2)q"" 1 . 
2 2 iin 2 


Note that d, is odd, and d, is even, hence 


desired inequality (14) is as follows 











2 


(g—2)q"™"-1|_ (q-2)q”"-1 
2 | 2 








5) 





(q-lq"™*-1]_ (q-lq"*-2 
5 | 5 


Consequently, the inequality (14) takes the form 
1 [ee =| -(-2)q"'-1 
- < 





2 2 


It is easily seen that it is equal to the inequality g=>3.e 


Consequence. If q>3, then there is a strict and when q=3, then there is an equality in (14). 
Table | lists parameters of some PM-codes. The top three lines contain the parameters q, m, n of the code in 


question RM ,(r,m) . The following three lines contain the values: k, is code dimension, d; is minimum code distance, 
and t; is the number of recoverable errors for the codes RM_,(1,m). And the following three lines show similar values of 


ko, do, tz for the codes RM, (2,m). 


Table 1 


Parameter values of some RM-codes 


Eee 
mcm i 


EP Poe a Do oe ie ane [a [ee a Pa 
Pe Pe fof Ps Te Popa [fe fo fs fa 
2 Fas fm [om [os [oe [aes [aes [as [as [re [as 
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Now, we introduce the analog of the differentiation operator D, operating in the space of polynomials (see 
(2)), for the space F) where n = q”. The vector coordinates from F’ will be numbered by vectors from the ordered set 
F." (see (10)). Let us consider the shift operator t, : F’ — F} following the formula 
T5 (4) = (Gy, poe 45) 


€F”". Note that the shift operator t- is mixing bijection. The linear 
5 p - g OY 


where @=(a,,...a, )EF’, baGe 5b, 
operator A; : F; + E,, that is the analog to D, , is defined by the formula: 
A; (@)=1;(@)-@, @=(a,,...4, )EF. (15) 
We call A, (a) a vector derivative of vector a along b. 
Lemma 3. Let us consider the polynomial f € ‘ie [x1,%)5--.%, |, vector b =(,, ...,b,) € F)’, operators A; , D; 
and C . Then, 
5 (C(f)) = CFE), CCD, A) =A (CP) - (16) 
The proof is carried out by the direct calculations, and for g =3 is available in [6]. 
Note that from (2) and (16), it follows that if C(w) € RM, (2,m), then A; (C(w)) € RM, (1m). 
Below, we consider examples of practical applications of the theoretical results obtained. 
Data partitioning scheme. For partitioning and restoring data in the proposed scheme, we use [n,k,,d,], - 


code, RM,(1,m) and [n,k,,d,],-code RM,(2,m), specified over the derived Galois field, of fuzzy cardinality. 
Values q and m are parameters of this scheme. 
Data partitioning algorithm. 
Input: information vector we Fe of the code RM_,(2,m) and an ordered set of basic vectors 
B= {b =(G1,0),...5,) Eom (17) 


That is a private key of the scheme under consideration. 





Output: vectors S, € F""', i=1,m. 
Step 1. Let us assign the information polynomial w= w(x) to the input vector w and encode it using (11) into 
the vector C(w) € F’ of the code RM ,(2,m). 
Step 2. Let us form m of the vector derivatives (see (15)): 
A, (C(w)) = C(D, (w)) €E’, i=l, b eB. 





Note that C(D, (w)) = RM, (1,m) . 
Step 3. We concatenate each vector C(D, (w)) €F’, i =1,m with the coefficient f,, 9) = w(0) of the code 


vector C(w) : 
5,=CD; 0) Il foo. € Ey” - 


Then, vectors S, e F"", i= lm are transmitted along m different communication lines. Obviously, during the 
i ek, g y g 





transmission, vectors S, , 1=1,m can be distorted. Thus, from the communication channel, vectors S, ‘will be obtained: 


S,'=(C(D, Oo) "IF. 


00..00 


where (C(D, (w)))' is supposedly distorted vector C (D; (w)), scalar f 


00..00 


‘ee 
eF, i=lLm. 


is supposedly distorted value fo, 9). Let us 


denote scalar f. corresponding to S,'by f- 


Data recovery algorithm. _ 

Input: vectors S,', i=1,m and private key B (see (17)). 

Output: vector ' € F . 

Step 1. We isolate two components: (C(D, (w)))'eF’ and f, er I,m from each vector Se Lm. 

Step 2. We direct the vectors (C(D, (w)))' to the decoders of the code RM, (l,m). Note that you can use 


arbitrary decoders, for example, [16], [18]. At the output of the decoders considered, polynomials 
Dz (w) CE) [x,,X25.%y |>4=L...,m are generated. 


ee 


00..00 m 


Step 3. We develop vector (f.__, 


00..00,1 ? 4 00..00,2 7 


of the code RM,(0,m) which actually coincides with the code of m-tuple repetition. The result of this decoder 


) from the coefficients f, ; », and feed it to the decoder input 


performance is the scalar fo, oq - 
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Step 4. We substitute the values of the coefficients of the polynomials D, (w) CF)? [2,,%)5..%,], i= I,m and 
the key B (see (17)) into the formulas (7) and (8), and derive the polynomial f(x) from the results obtained. Then we 
calculate the desired polynomial w'(x) = f(X)+ foo 00 - 

Step 5. The message recipient is given the information vector w' € F corresponding to the polynomial w'(x) . 


Note 1. The correctness of the data recovery algorithm depends on the number of errors that damage the 
vectors S, € i during their transmission along the communication lines, as well as on the correcting capacity of the 


first-order PM-codes used. It is plain that if the decoders used reconstruct the vectors C(D, (w)) and the value fo9_o5 
correctly, then the initial data restoration using the results of Theorem 1 will be correct, hence the vector w'(x) 
obtained at the output of the data recovery algorithm will coincide with the original information vector we F? . Note 
that operation of the decoder for recovery C (D; (w)) is correct, if 

Vi=lm: d,(C(D; (w)),(C, (w))') $(@,-)/2], 
where d,(x,y) is Hamming distance between the vectors x,y. The scalar fo, 9) is restored correctly if the vector 


CF sai, 


00.001 7% 00..00,2 7 * 


ey ; ae formed in step 3 of the algorithm contains fewer than m/2 coordinates differing from the 
value fn 9 =W(0). If, in S, € EY i= l,m, more errors occurred than can be restored by the decoders used, then the 


_ ky . 
recovery we F’? is not guaranteed. 


Note 2. The proposed partition and recovery algorithms employ Reed-Muller codes of both the first and second 
orders, but the decoders are used only for the first-order codes. 

Note 3. The transmission confidentiality is provided not only by the need for knowledge of the key, but also by 
using several communication lines, because in this case, data interception is more challenging for an intruder than 
illegal data capture from a single communication line. 

Decoder of Reed-Muller codes of the second order. First, let us consider the idea of the decoding algorithm 


organization, and then we describe the algorithm step-by-step. Let us fix some basis B = {b, = (;,0},...55;,) ©F"}i4.n 
in the space F.” where q = p*, pis an odd prime. Suppose that the decoder input receives Y =C(w)+é@(e F’) where w 
is the information polynomial, C(w) is the code vector of the [n,k,,d,],-code, RM ,(2,m) ,@ € F;is the error vector for 
which 
wt, (€) St, , (18) 
where wt,(.) is Hamming weight, t, = | (d, —l)/ 2| . Employing vector Y , we construct m vector derivatives calculated 
in the basic directions using the operator A, : 
A; (Y) =A; (C(w) +2) =A; (C(w)) +A; @), i=1L....m 

each of which is the vector A; (C(w)) ¢ RM, (1,m) distorted by the error vector A; (€) <F’, and can be unerringly 
decoded by an arbitrary decoder of the [n,k,,d,], -code RM, (1,m) operating up to half the code distance (see, e.g., [16], 
[18]), if the error level is less than ¢, = | (4, -1)/2| , le. when 

wt, (A; (2) St. (19) 


If the vectors are decoded correctly, then the desired information polynomial of the code can be reconstructed using 
Theorem | up to one coordinate which can be then found, for example, by maximum likelihood decoding. Thus, for 
proper decoding Y , according to the proposed scheme, the fulfillment of the condition (19) is required. 

Algorithm for decoding the code RM ,(2,m) . 


Input: parameters of [n,k,,d,], -code RM,(2,m), Y =(¥,,,Y,,.--Y, )€ EF. 
Output: decoded information vector w . 
Step 1. Let us fix some basis B = {b, =(b;,4,...b,,)€F"},._, Of the space F” and calculate the derived 


m 


vectors along all the directions b, ep: 
A, (Y)=1%,)-Y. 
Step 2. Let us decode A; (Y), i=l,...,m, using the arbitrary decoder of RM, (1,m) -codes that operates up to 
half the code distance, and as a result, we obtain vectors p, and their polynomial representations 


eet _— 
Dj, (%) = OX, FOX +... +O, Xp, + Oly CF Emer ee eee 


mm 
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Step 3. Using the polynomials p,; (x), 7=1,...,m, we find the polynomial f(x) with the intercept zeroth-order 


term from the formulas (7) and (8). 
Step 4. For all ze F,, we calculate 


> 





¥@)=YICU@+2, -%, 
where C( f(x) PE) is o,- th coordinate of the geet C(f(x)+z) (see (11)). Let us denote z, the value z on which 
the function ‘Y¥(z) attained its minimum. 
Step 5. The decoding result is vectorw that one-to-one corresponds to the information polynomial 
W(X) = f(%)+2y. 
Theorem 2. So that the derived algorithm of decoding the code RM, (2,m) can rectify all errors, it is sufficient 


to meet the conditions (18) and 
wt, (€) St, /2, (20) 


where t, =| (d,-1)/2], d, is minimum distance of the code RM _,(1,m). 


Proof. At step 2, the decoder input receives vectors 
A, (Y)=1, (Y)-Y¥ =; (Cl) +1; (€)-C() -€ = A; (C(w)) +A; (2). 





Recall that A; (C(w)) ¢ RM, (1,m) , and the decoders of the code RM,(1,m) operating up to half the code distance 
rectify up to f, errors in the codeword. From the condition (20), it follows that 

wt, (A; (€)) < wt, (t, (€))+ wt, (@) = 2wt,(@) St, 
So, vectors D;, which are developed at Step 2 coincide with A; (C(w)). Hence it follows that at Step 3, owing to 


Theorem 1, the polynomial {(¥) = w(¥)—w(0) is formed. 
From the condition (18), it follows that the value z, calculated at Step 4 is equal to the constant term w(0) of the 
desired information polynomial w(x) . Thus, the desired information vector w is obtained. 

Note that for the correct decoding according to the proposed scheme, it is required to meet the condition (20), 
from which (19) follows, although the condition (18) is more natural. Let us consider the relationship between these 
conditions. 

Lemma 4. For the codes RM, (2,m), in the case of g = 3, the condition (18) and (20) are equivalent, and in the 
case of g >3, when the condition (18) is satisfied, the fulfillment of the condition (20) is not guaranteed. 

Proof. From the consequence of Lemma 2, we obtain that at q =3, the equality ¢,/2 =t, is true; ie. the right- 
hand sides of the inequalities (18) and (20) coincide; hence, the fulfillment of one of them implies the fulfillment of the 
other. At g>3 from the consequence of Lemma 2, we obtain thatz,/2<t,, viz from the fulfillment of (18), the 
fulfillment of (20) does not follow. 

Note I. In [5-6], the decoder of the RM,(2,m)-code is described, where, as in the proposed algorithm, 
derivative vectors are constructed for the noisy codeword. They are decoded by the maximum likelihood algorithm, and 
then the desired information word is restored from the values obtained. However, the derived vectors are constructed in 
all 3” possible directions, not only in basic ones, in the decoder from [5-6]. And another mechanism is used to obtain 
the desired information vector. 

Note 2. For codes RM ,(2,m) ,q =3, the proposed decoder operates up to half the code distance. For codes 


RM _,(2,m) ,q >3, the proposed decoder does not guarantee the correction of all errors, the number of which is less 


thant, , but the decoder will operate well if the weaker condition (20) is satisfied. Note that the use of the proposed 


decoding scheme in case of the fields with cardinality of more than three may be appropriate, despite the determined 
limitation, for the following reasons. First, the theory of the decoders of the second-order PM-codes is ill-defined. But, 
if there is a first-order decoder, then the proposed decoder which is a suspension over it can fill in this gap. Secondly, 
when using the communication channels, the error probability in which is such that (20) is satisfied, the transition from 
the first-order PM-codes to the second-order codes reduces redundancy (see Table 1). 
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Conclusion. Theoretical results associated with restoring polynomials of several variables over Galois fields of 
fuzzy cardinality by their derivatives are obtained. As practical applications to the results obtained, a data partition 
scheme and a second-order PM-code decoder are proposed. In future, it is instructive to study the process of a 
polynomial recovery from distorted derivatives, and to develop appropriate modifications of the real-world applications 
proposed in this paper. 
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